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INTRODUCTION  -  subject.  Improvised  explosive  devices  (lEDs)  produce  head  injuries  in  nearly  a 
majority  of  surviving  soldiers.  Most  brain-injured  soldiers  do  not,  however,  receive  rapid  brain  imaging  studies  as 
would  their  civilian  counterparts.  Instead  they  are  flown  to  rear-echelon  medical  service  centers  such  as  in 
Germany  for  those  studies,  as  well  as  additional  medical  care.  This  is  because  magnetic  resonance  imaging  is 
unavailable  in  or  near  the  theatre  of  war  (and  a  bad  idea  for  many  patients  because  of  the  presence  of  metal 
fragments)  and  because  of  the  lack  reliably  operational  CT  machines  at  medical  centers  in  Iraq.  In  the  mean 
time,  due  to  the  lack  of  adequate  neuro-imaging,  acute  neurosurgical  care  in  response  to  closed  TBI  often 
requires  a  complete  hemicraniectomy  (removal  of  one  side  of  the  patient’s  skull)  simply  to  assay  the  extent  and 
location  of  injured  brain.  There  exists,  therefore,  a  critical  need  for  robust  brain  imaging  systems  at  and  near 
the  battlefield.  Also,  more  subtle  brain  injuries  are  apparent  clinically  but  difficult  to  definitively  diagnose,  hence 
treat.  These  gaps  in  patient  care  reduce  the  quality  of  care  and  potentially,  therefore  the  quality  of  life  of  injured 
soldiers.  This  gap  also  defines  a  critical  need  for  rugged,  field  deployable  systems  capable  of  imaging  injured 
brain. 

INTRODUCTION  -  purpose.  For  a  variety  of  reasons  (summarized  in  a  recent  book  chapter  of  mine  - 
Mourad,  2013)  it  is  reasonable  to  expect  that  changes  in  the  stiffness  of  brain  accompany  TBI,  and  that 
ultrasound-based  ‘sonoelastic’  imaging  modalities  responsive  to  some  measure  of  stiffness  might  offer  a  useful 
means  for  imaging  the  gross  and  subtle  changes  to  brain  due  to  TBI.  Use  of  such  systems  in  and  near  the  field 
should  improve  clinical  outcome  for  patients  suffering  from  TBI.  Our  long-term  goal  is  to  develop  a  field 
deployable  brain  imaging  system  responsive  to  brain  stiffness.  Our  system  would  do  so  through  the  use  of 
ultrasound  applied  transcranially. 

INTRODUCTION  -  scope  of  the  research.  To  met  the  goal  of  this  proposal  we  sought  to  test  the 
following  hypothesis:  Rodent  brain  after  closed  TBI  (blast,  controlled  cortical  impact  (CCI);  ischemic  stroke) 
manifest  spatial  patterns  of  endogenous  and  exogenous  tissue  displacement,  hence  patterns  of  brain-tissue 
stiffness,  sufficient  to  map  and  identify  the  extent  of  each  of  hemorrhage,  edema,  and  TAI,  as  tested  via  the 
following  iterative  aims. 

Specific  Aim  #1 :  Quantify  via  histology  and  MR  the  spatial  position  and  extent  of  each  of  hemorrhage, 
edema  and  TAI  within  each  of  the  blast  and  CCI  models  of  brain  injury. 

Specific  Aim  #  2:  Directly  measure  the  intrinsic  stiffness  of  brain  tissue  of  different  types:  normal  gray  and 
white  matter;  hemorrhagic  brain;  edematous  brain;  traumatized  axons. 

Specific  Aim  #  3:  Use  displacement-based  ultrasound  to  image  intrinsic  brain  tissue  stiffness  of  normal 
versus  traumatized  rodent  brains  and  their  component  parts. 


4 


BODY 


Specific  Aim  #1:  Quantify  via  histology  and  MRI  the  spatial  position  and  extent  of  each  of 
hemorrhage,  edema  and  diffuse  axonal  injury  (DAI)  within  each  of  the  blast  and  CCI  models  of  brain 
injury,  with  comparable  studies  for  ischemic  stroke. 

Task  (1):  Generate  closed  TBI  (mild  and  moderate-to-  severe)  in  rats  using  blast  methodology. 

1)  Thanks  to  administrative  issues  we  have  delayed  this  task  to  start  in  summer/fall  2014 

Task  (2):  Generate  closed  TBI  (mild,  moderate-to-severe)  using  CCI  methodology  in  ‘mito’  mice 
and  in  rats  and  ischemic  stroke  using  an  occlusion  model  in  mice. 

2)  As  reported  last  year,  we  have  achieved  these  goals. 

Task  (3):  MR  image  rodent  brains. 

3)  UVa  has  performed  its  all  of  its  of  MRI  studies  of  CCI  rats  (Figure  1a,b).  They  have  finished  co¬ 
registration  of  those  MR  images  with  immunohistochemical  images,  in  order  to  develop  direct  correlation  of  the 
brain-injury  features  in  each  imaging  modality  (Figure  1c).  To  do  so,  they  have  developed  an  algorithm  to 
spatially  normalize  each  of  brain  histological  slides  and  corresponding  MRI,  achieved  primarily  through 

the  various  tools  found  within  the  Advanced  Normalization  Tools  (ANTs)  repository.  Due  to  administrative 
issues  this  work  is  on  hold  until  Summer  2014. 
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Figure  1.  Images  of  rat  brain  after  (sham)  controlled  cortical  impact  (CCI)  injury.  (A)  shows  various  MRI 
images  of  coronal  views  of  sham  injured  rat  brains.  (B)  shows  various  comparable  MRI  images  of  injured  rat 
brains.  The  injury  manifests  in  the  upper  right-hand  quadrant  of  the  brain,  encompasses  the  cortex  down  to  the 
hippocampus.  ‘T 1  ’  highlights  the  basic  structural  abnormality  generated  by  the  CCI,  ‘T2’  highlights  fluid  buildup 
at  the  site  of  CCI,  ‘Flair’  refers  to  fluid  attenuated  inversion  recovery  that  suppresses  fluid-specific  signals  to 
highlight  changes  in  tissue  structure,  ‘SWI’  refers  to  susceptibility  weighted  imaging  that  is  sensitive  to  venous 
blood  and  hemorrhage.  (C)  shows  co-registration  of  MRI  images  with  a  rat-brain  atlas,  each  of  the  major 
components  of  brain  highlighted  in  different  colors.  Follow-on  co-registration  of  the  atlas  with  histological 
images  will  complete  a  useful  triad  of  MRI/Atlas/histology. 
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Task  (4):  Immunohistochemical  assays. 

4a)  Each  of  UW  and  UVa  have  performed  all  of  their  immunohistochemical  assays  for  CCI  rats.  At  UW 
we  are  using  those  results  to  finish  two  papers  detailing  the  ability  of  different  modalities  of  ultrasound  (shear- 
wave  imaging  (SWE)  and  vibro-acoustography)  to  detect  traumatic  brain  injury  (TBI),  described  below  under 
Specific  Aim  #3. 

4b)  We  at  UW  have  collected  wide  field  fluorescent  microscopy  brain-tissue  images  from  mito-mice 
subjected  to  severe  TBI  via  CCI  (Figure  2a).  (These  mice  have  the  unique  property  that  the  mitochondria 
within  the  neurons  in  their  brain  fluoresce  when  healthy.)  We  have  identified  unanticipated  increases  in 
fluorescence  at  the  site  of  injury  pointing  to  possible  decreases  in  the  ability  of  the  neurons  and  axons  to 
transport  mitochondria,  or,  to  changes  in  mitochondria  morphology  and  function.  Moreover,  we  observed 
changes  in  both  ipsilateral  as  well  as  contralateral  fluorescence  due  to  injury,  with  comparable  levels  of 
fluorescence  both  acutely  as  well  as  24  hours  after  injury  (Figure  2b).  Interestingly,  we  did  not  see  this  kind  of 
damage  for  mild  TBI  (data  not  shown). 
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Figure  2.  mitoMice  brains  -  lighter  gray  has  greater  concentration  of  mitochondria  than  darker  gray. 

(A)  This  figure  shows  a  coronal  slice  of  brain  subjected  to  actual  CCI.  Note  the  enhanced  fluorescence  within 
the  hippocampus  (as  well  as  edema  and  large-scale  necrosis),  suggestive  of  mitochondrial  dysfunction  and/or 
reduced  transport.  (B)  Averaging  over  all  data  (N  =  3)  we  observe  increase  in  fluorescence  in  both 
hemispheres  due  to  a  single-hemisphere  contusion  injury  that  affects  axonal  transport  in  the  upper  quadrants 
of  both  hemispheres. 
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Specific  Aim  #  2:  Directly  measure  the  intrinsic  stiffness  of  brain  tissue  of  different  types: 
normal  gray  and  white  matter;  hemorrhagic  brain;  edematous  brain;  ischemic  brain;  traumatized 
axons. 

Task  (5):  Calibrate  indenter  system  and  apply  it  in  vitro  as  well  as  to  rodent  brain  tissue 
samples  dominated  by  one  of  each  of  normal  gray  and  white  matter;  hemorrhagic  brain;  edematous 
brain;  traumatized  axons. 

5)  This  system  is  calibrated  and  working  after  quite  extensive  development  of  the  technique.  We  have 
shown  statistically  significant  changes  in  brain  tissue  viscoelastic  structure  at  the  site  of  CCI  -  Figure  3.  We 
have  published  this  paper  whose  results  we  described  in  detail  in  the  last  annual  report  -  (Gabler  et  al,  2013  - 
attached). 


A _ 

Figure  3  shows  our  indenter  results  for  CCI  rats  from  the  final,  published  paper  -  Gabler  et  al,  2013. 

Figure  (A)  shows  measurements  of  rat-brain  displacement  (horizontal  axis)  versus  measured  elastic  response 
in  milliNewtons,  all  as  a  function  of  time,  along  with  associated  mathematical  modeling.  Figure  (B)  shows  the 
results  that  demonstrate  changes  in  the  viscoelasticity  of  brain  after  CCI  at  the  various  spots  within  the  brain 
from  which  data  was  collected. 
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Specific  Aim  #  3:  Use  displacement-based  ultrasound  to  image  intrinsic  brain  tissue  stiffness  of 
normal  versus  traumatized  rodent  brains  and  their  component  parts. 

Task  (6):  Optimize  commercial  ultrasound-imaging  device  based  on  measurement  of 
exogenous  and  endogenous  brain  palpation. 

6a)  We  have  established  a  procedure  for  optimizing  the  analysis  of  images  from  the  commercial 
ultrasound  imaging  system,  based  on  exogenous  brain  palpation  by  the  SuperSonic  Imagine  (SSI)  machine  via 
the  acoustic  radiation  force  and  subsequent  monitoring  of  the  propagation  of  the  shear  wave.  We  have 
performed  this  procedure  for  stroke  animals  and  published  the  results,  where  this  grant  supported  the  analysis 
of  that  data.  In  Figure  5  from  the  second  annual  report  we  described  the  methodology.  We  have  now 
published  a  paper  describing  the  application  of  this  methodology  to  TBI  detection  in  each  of  mice  and  rats, 
discussed  below  in  Task  8. 

Task  (7):  Optimize  research  ultrasound-imaging  device  based  on  measurement  of  exogenous 
and  endogenous  brain  palpation. 

7a)  As  reported  earlier,  we  can  track  shear-wave  propagation  using  our  research  ultrasound-imaging 
device  (the  VUE  system)  -  Fig.  4a, b.  We  have  used  this  research  methodology  to  analyze  shear-wave 
propagation  within  mouse  brain  ex  vivo  by  varying  how  the  shear  wave  entered  brain  (Figure  4c, d).  We 
discuss  the  details  of  this  methodology  in  the  attached  senior  thesis  (Yayun,  2013).  We  report  preliminary 
results  below. 
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Figure  4.  Quantification  of  shear  wave  propagation  in  alg  nate.  We  can  create  a  focal  and  vertically 
elongated,  vertical  push  within  alginate  (and  tissue  -  below)  that  (A)  generates  a  horizontally  propagating 
shear  wave,  whose  speed  we  can  quantify  and  translate  (B)  through  custom  software  that  cross-correlates 
temporally  adjacent  ultrasound  images  along  a  given  horizontal  line.  The  local  stiffness  of  tissue  is 
proportional  to  the  square  of  that  local  speed  of  the  shear-wave  propagation.  (C,  D)  We  studied  shear-wave 
propagation  into  damaged  brain  ex  vivo  from  a  variety  of  angles  relative  to  the  injury  in  brain  in  an  attempt  to 
develop  an  image  of  the  contusion  injury  created  by  the  CCI  mechanism. 
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7b)  We  have  tested  on  sham  and  actual  TBI  rats,  both  acute  and  24  hours  after  injury,  an  exogenous 
brain  palpation  system  based  on  vibrating  brain  with  focused  ultrasound,  whose  methodology  we  review  in 
Figure  5a  and  demonstrate  capable  of  measuring  acoustic  emissions  from  brain  (Figure  5b).  We  report  results 
below. 


vibroacoustography:  vibrating 
tissue  with  ultrasound 
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Figure  5.  Vibro-acoustography  as  a  means  of  assay  ng  for  the  presence  of  traumatic  brain  injury. 

(A)  In  vibro-acoustography,  a  dual-transducer  device  (‘transducer  (annular  array)”)  emits  two,  separate 
pulses  of  ultrasound  at  two  different,  but  close  frequencies  -  pulse  train  ‘A’  and  ‘B’.  At  their  mutual  focus, 
high  frequency  oscillations  move  the  tissue  slightly;  their  beat  frequency  generates  a  low-frequency 
emission,  itself  sampled  from  a  separate  hydrophone  placed  either  in  the  water  tank  or  within  the  iFU 
source  itself  for  in  vivo  applications.  (B)  This  figure  shows  our  device  and  real  data  rather  than  the 
conceptualization  of  Figure  6A,  highlighting  a  measured  waveform  that  demonstrates  the  signal 
generated  by  the  two,  spatially  overlapping  high  frequency  iFU  signals,  and  small  amplitude,  low- 
frequency  acoustic  emissions. _ 
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7c)  To  create  endogenous  ultrasound  imaging  results  we  have  demonstrated  that  we  can  measure  the 
pulsation  of  mouse  brains  completely  enclosed  within  their  crania,  pulsations  created  by  both  their  breathing 
and  their  heart  beat  (Figure  6).  Unfortunately,  despite  this  finely  detailed  work  we  do  not  see  a  difference  in 
pulsatility  between  the  TBI  and  contralateral  side  of  mice  and  rats,  neither  acutely  nor  24  hours  after  injury. 
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Figure  6  demonstrates  our  ability  to  collect  tissue  pulsatility  images  of  mouse  brain  in  vivo  within  the 
mouse’s  completely  enclosed  cranium.  We  developed  a  measure  of  the  average  brain-tissue  displacement 
within  a  small  portion  of  mouse  brain,  following  the  method  of  Kucewicz  et  al,  2006,  2007,  2008.  That 
measure,  as  a  function  of  time,  demonstrates  both  respiratory  peaks  (e.g.,  when  the  mouse  breathes)  as  well 
as  cardiac  peaks  (when  the  heart  beats),  when  shown  in  time  series  (A)  or  in  spectral  analysis  as  shown  in  (B). 
Note  that  the  respiratory  peaks  can  approach  1  micron  in  diameter,  while  the  cardiac  peaks  measure  one 
tenth  of  that,  on  average.  (C)  Filtering  out  the  respiratory  peaks  highlights  the  cardiac  cycles  only  -  inset  of  ‘C’ 
as  well  as  (D). 
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Task  (8):  Image  normal  versus  damaged  rodent  brains  with  commercial  and  research 
ultrasound  devices. 

(8a)  We  show  statistically  significant  variations  in  acoustic  emissions  from  within  severe  TBI  rat  brain 
relative  to  sham  TBI  rat  brain,  with  data  collected  a  few  millimeters  away  from  the  TBI  site  (Figure  7,  for  acute  as 
well  as  for  results  measured  24  hours  after  TBI).  We  will  finish  writing  up  these  results  (Suarez  et  al,  2014  -  see 
draft  in  Appendices)  then  begin  design  of  a  system  capable  of  application  to  humans. 
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Figure  7  -  experimental  results  for 
vibroacoustography.  (A)  shows  a  schematic  of  the 
experimental  setup.  The  large  ellipse  represents  the 
portion  of  the  skull  from  which  skin  and  fascia  were 
removed,  exposing  the  suture  lines.  The  circle  represents 
the  placement  of  the  brain  contusion  after  performing  a 
small  craniectomy.  Black  stars  denote  superficial 
projection  of  placement  of  the  vibroacoustography  focus 
at  8mm  deep  within  the  cranium.  (B)  demonstrates  that 
the  acoustic  emissions  associated  with  acute  sham  (dark) 
are  larger  than  that  measured  from  acute  severe  CCI  for 
most  of  the  interrogation  frequencies  considered  (from 
200-270  kHz).  Figure  (C)  demonstrates  comparable 
results  for  brain  24  hours  after  TBI,  without,  however,  a 
marked  frequency  dependence. 
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(8b)  We  have  finished  use  of  the  commercial  SSI  system  without  any  modification  on  TBI  mice  and  rats 
-  Figure  8,  with  a  paper  describing  this  work  in  press  (Xu  et  al,  2014,  see  attachment  in  appendix).  Specifically, 
these  results  demonstrate  changes  in  hemispherically  averaged  shear  modulus  for  each  of  the  hemispheres  due 
to  ipsilateral  TBI.  The  ipsilateral  decrease  in  shear  modulus  relative  to  control  is  consistent  with  the  combination 
of  focal  ischemia  and  wide  spread  edema  formation  after  CCI,  a  well-established  phenomenon  in  this  animal 
model  for  TBI.  Unexpectedly,  we  also  see  a  contralateral  increase  in  shear  modulus  relative  to  control  by  24 
hours  after  CCI  injury,  consistent  with  other’s  observation  of  reduction  of  blood  flow  contralateral  to  injury  as  well 
as  elevated  global  ICP. 


Figure  8  demonstrates  our 
ability  to  detect  changes  in 
intra-hemispheric  brain 
stiffness  after  TBI  using  our 
commercial  imaging  system. 
Average  values  of  ‘shear- 
modulus’  values  (normalized  by 
the  average  value  for  the  entire 
given  brain)  within  a 
hemisphere  of  mouse  or  rat 
brain,  averaged  across  five 
(A)  mice  and  (B)  rats 
demonstrate  that  the  ipsilateral 
hemisphere  has  a  lower  mean, 
normalized  shear  modulus  than 
the  contralateral  hemisphere  by 
24  hours  after  CCI  injury. 

These  values  change  also 
relative  to  sham  TBI  values, 
which  themselves  agree 
between  hemispheres.  We  see 
comparable  values  for  mice,  as 
discussed  in  Xu  et  al  (2014). 


Sham  Acuta  24-hour  Sham  Acuta  24-hour 

contra  contra  contra  tpal  Ipal  ip*l 


A 


Sham  Acuto  24-hour  Sham  Acuto  24-hour 

contra  contra  contra  ipai  Ipat  ipai 


B 
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(8c)  We  have  finished  use  of  the  research  VUE  system  on  TBI  mouse  brain  ex  vivo  -  Figure  9,  as 
described  in  the  senior  thesis  (Yayun,  2013).  These  preliminary  results  suggest  that  we  can  image  small  brain 
contusions  with  shear-wave  based  methods,  although  translating  this  to  transcranial  results  has  thus-far 
proved  quite  difficult. 


Shear  Modulus  Stat 
Parameter  5 


mean  =  5.977 


mean  =  6.073 


]  Control  vs.  TBI  Mouse 


T 

o 

i 

mean  z 

Control  Left  Control  Right  TBI  Normal 

P  =  0-102  A 


Shear  Modulus  Stat 
Parameter  5 


TBI  Hemisphere 


]  TBI  Mouse 


Normal  Hemisphere 

p  =  0.004 


B 


Figure  9  -  results  of  analysis  of  ex  vivo  shear-wave  propagation  in  mouse  brain.  Figure  (A)  shows  that 
the  measured  shear  modulus  of  each  hemisphere  of  control  brains  are  comparable  to  each  other  and  slightly 
less  stiff  than  the  hemisphere  of  TBI  brain  contralateral  to  hemisphere,  consistent  with  our  published  results 
using  the  commercial  system.  Figure  (B)  shows  that  the  TBI  hemisphere  has  reduced  shear  modulus 
compared  to  the  contralateral  (here,  ‘normal’)  hemisphere,  again  consistent  with  our  published  results  using  the 
commercial  system. _ 


Task  (9):  Compare  images  of  strain,  and  of  bulk  and  shear  modulus  of  blasted  rat  brain  with 
histological  and  MR-based  images  of  hemorrhage,  edema,  ischemia  and  TAI. 

9)  We  intend  to  start  this  research  in  Summer  2014,  after  we  deploy  our  VUE  research  device  at  the 
University  of  Virginia. 

Task  (10):  Compare  images  of  strain,  and  of  bulk  and  shear  modulus  of  brain  tissue  with  images 
of  bulk  and  shear  modulus  derived  by  assigning  directly  measured  values  of  bulk  and  shear  modulus  to 
portions  of  histological  images. 

10)  We  intend  to  start  this  research  in  Summer  2014  after  we  deploy  our  VUE  research  device  at  the 
University  of  Virginia. 

Task  (11):  Generate  first-order  specifications  for  a  device  capable  of  imaging  TBI  in  humans. 

11)  We  can  now  specify  one  of  what  we  anticipate  to  be  three  devices  capable  of  imaging  TBI  in 
humans.  That  device  will  emit  diagnostic  ultrasound  with  a  frequency  around  2  MHz,  with  its  data  analyzed  via 
‘tissue  pulsatility  imaging’  [TPI]  algorithms  we  have  developed.  The  other  device  we  anticipate  will  involve 
vibro-acoustography  in  combination  with  TPI.  We  are  less  optimistic  about  using  shear-wave  methods  to 
image  transcranially. 
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KEY  RESEARCH  ACCOMPLISHMENTS 


•  We  have  developed  and  demonstrated  in  vivo  a  new  means  of  detecting  severe  TBI  both  acutely 
and  after  24  hours,  based  on  generation  and  detection  of  focal,  low-frequency  vibration  within  brain 
via  a  process  known  as  vibro-acoustography. 

•  We  have  demonstrated  our  enablement  of  a  novel  ultrasound-based  imaging  system  capable  of 
translating  diagnostic  ultrasound  images,  derived  transcranially,  into  images  of  endogenous  brain- 
tissue  displacement  that  highlight  focal  TBI,  based  on  our  ability  to  image  mouse-  and  rat-brain 
pulsatility.  Unfortunately  these  results  have  not  yet  demonstrated  sensitivity  to  the  presence  of  TBI 
on  length  scales  of  millimeters. 

•  We  have  observed  changes  in  hemispherically  averaged  brain-tissue  stiffness  in  TBI  generated  by 
the  controlled  cortical  impact  (CCI)  method  in  mice  and  rats  consistent  with  physiological  changes 
in  brain  observed  in  rodent  models  of  TBI  as  well  as  in  patients  with  TBI.  These  observations 
include  ipsilateral  edema  and  hemorrhage  formation  (observed  in  rodents  and  in  patients)  and 
short-term  contralateral  reduction  in  blood  flow  (observed  in  patients  and  sometimes  in  rodents) 
along  with  a  reduction  in  cerebral  perfusion  pressure.  We  have  a  paper  that  describes  these  results 
in  press. 

•  We  have  demonstrated  in  preliminary  analysis  that  we  can  produce  hemispheric-averaged 
estimates  of  mouse  brain  stiffness  in  a  using  our  research  device,  ex  vivo. 

•  We  have  developed  a  means  of  co-registering  MR  and  immunohistochemical  images,  in  order  to 
develop  direct  correlation  of  the  brain-injury  features  in  each  imaging  modality.  The  tools  we  are 
creating  to  accomplish  this  task  will  likely  find  clinical  application,  for  bringing  together  different 
imaging  modalities  for  TBI  and  other  etiologies. 

•  We  have  demonstrated  our  direct  means  of  measuring  brain-tissue  stiffness  of  rat  brains,  data  that 
will  give  us  ground-truth  measurements  of  that  stiffness  for  comparison  with  our  ultrasound-based 
means  of  imaging  brain-tissue  stiffness  variations  in  the  context  of  TBI. 

•  We  anticipate  deployment  at  UVa  of  a  research  ultrasound  imaging  system  capable  of  generating 
tissue-pulsatility  images  and  possibly  vibroacoustography  maps  of  the  brains  of  rats  after  blast 
injury.  Successful  collection  and  analysis  of  these  data,  along  with  direct  measures  of  brain-tissue 
stiffness  and  associated  MRI  and  immunohistochemical  images,  will  allow  us  to  demonstrate  in  a 
direct  fashion  that  we  can  use  ultrasound  to  produce  high-contrast  images  of  damaged  versus 
normal  brain  in  the  context  of  TBI. 
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CONCLUSIONS 


We  have  generated  evidence  supportive  of  the  hypothesis  that  ultrasound-derived  brain-tissue  stiffness 
imaging  using  exogeneous  ‘palpation’  of  brain  with  diagnostic  ultrasound  can  successfully  detect  changes  in 
rat  brain  caused  by  both  ischemic  stroke  and  in  TBI.  This  uses  shear-wave  generation  and  tracking  via 
ultrasound  as  well  as  ‘vibrating’  brain  with  focused  ultrasound.  We  have  also  generated  evidence  supportive  of 
the  hypothesis  that  we  can  detect  endogeneous  palpation  of  rodent  brain  by  cerebral  blood  flow.  We  have 
made  direct  measurements  of  brain-tissue  stiffness  in  rats  ex  vivo  that  show  a  change  in  brain  tissue 
viscoelasticity  due  to  TBI.  We  have  completed  a  novel  means  of  co-registering  imaging  and  histological  data 
that  together  highlight  TBI.  We  have  begun  designs  of  two  ultrasound  systems  capable  of  creating  images  of 
TBI  in  humans. 

In  our  last  extension  year  we  will  generate  comparable  data  on  blast  mice  and  blast  rats. 
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44  Abstract 

45  Objectives:  Traumatic  brain  injury  (TBI)  can  cause  adverse  physiological  changes  in  fluid 

46  content  within  brain,  which  may  lead  to  changes  in  tissue  elasticity  (e.g.,  stiffness).  This 

47  study  evaluated  the  ability  of  ultrasonic  shear  wave  elastography  (SWE)  imaging  to  observe 

48  these  changes  in  brain  after  TBI,  in  vivo. 

49  Methods:  Mice  and  rats  received  a  mild  TBI  or  sham  surgery  and  were  imaged  acutely  or  24 

50  hours  after  injury  using  SWE  and  the  hemispheric  stiffness  values  were  compared. 

5 1  Results:  Stiffness  values  were  consistent  across  brain  hemispheres  of  sham  TBI  rodents.  By 

52  24  hours  after  TBI,  relative  brain-tissue  stiffness  values  for  mice  and  rats  each  decreased 

53  ipsilaterally  and  increased  contralaterally,  both  relative  to  each  other  and  as  compared  to 

54  sham  TBI  rodents  (p<0.05).  The  absolute  value  of  tissue  elasticity  increased  for  rats 

55  (p<0.05)  but  not  for  mice. 

56  Conclusion:  Differences  between  intra-hemispheric  stiffness  values  of  rodent  brain  by  24 

57  hours  after  mild  TBI  may  reflect  the  observed  edema  and  hemorrhage  ipsilateral  to  TBI,  and 

58  the  known  reduction  of  cerebral  blood  flow  in  both  brain  hemispheres.  If  these  hypotheses 

59  hold  true,  ultrasonic  shear  wave  elastography  may  offer  a  method  to  detect  adverse  changes 

60  in  fluid  content  within  brain  after  mild  TBI. 

61 
62 

63  Keywords:  Brain;  Edema;  Elastography;  Ultrasound;  Traumatic  Brain  Injury; 

64 

65 
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66  Introduction 

67  Incidence  of  traumatic  brain  injury  (TBI)  has  increased  significantly  in  recent  years 

68  owing  to  the  growing  use  of  improvised  explosive  devices  in  operations  in  the  Middle  East 

69  [1],  TBI  also  occurs  among  civilians,  primarily  caused  by  falls,  motor  vehicle  accidents, 

70  struck  by/against  events,  contributing  up  to  30%  of  injury-related  deaths  [2], 

7 1  Computed  tomography  (CT)  and  magnetic  resonance  imaging  (MRI)  are  the  current 

72  civilian  standards  for  imaging  traumatic  brain  injury,  with  CT  more  commonly  available, 

73  thanks  to  its  speed  and  relatively  low  cost.  It  is  therefore  generally  used  acutely,  if  available, 

74  while  MRI  finds  more  use  a  few  days  after  initial  injury  because  of  its  sensitivity  to  blood 

75  products  [3].  Both  imaging  technologies  are  often  insufficiently  accessible  for  military  TBI 
7  6  (Ken  Curley,  MD,  personal  communication)  with  the  additional  requirement  of  patient 

77  transport  to  the  imaging  site,  relevant  in  either  the  civilian  or  military  context.  Transcranial 

78  ultrasound  imaging  offers  a  rapid,  portable  and  point-of-care  system  for  imaging  brain  but 

7 9  lacks  sufficient  contrast  to  provide  much  diagnostic  utility  [4] . 

80  Here  we  consider  ultrasonic  shear  wave  elastography  (SWE)  as  a  possible 

81  alternative  or  adjunct  to  CT  and  MRI.  SWE  creates  images  of  Young’s  Modulus,  a  measure  of 

82  the  resistance  to  such  that  larger  values  describe  stiffer  tissue.  Therefore,  maps  of  Young’s 

83  Modulus  describe  the  spatial  dependence  of  tissue  shear  elasticity  or  stiffness.  The  SWE 

84  system  works  by  generating  within  tissue  a  rapid,  net  axial  displacement  of  tissue.  This 

85  spatially  focal  impulse  force  then  creates  a  shear  wave  that  propagates  outward  from  the 

86  ‘pushed’  region  [5-7],  SWE  systems  use  diagnostic  ultrasound  to  track  the  propagating 

87  shear  wave,  with  the  existing  commercial  system  from  SuperSonic  Imagine  (SSI)  translating 

88  the  square  of  that  local  speed  into  a  local  value  of  Young’s  Modulus  [5]. 

89  SWE  via  the  SSI  clinical  system  has  found  well-documented  use  on  humans  for  large 

90  and  primarily  homogenous  organs  such  as  breast  [5, 8-9].  Studies  exist  that  have  used  the 

91  clinical  system  in  animal  models  of  disease  such  as  kidney  [10]  and  ischemic  cerebral  stroke 

92  [11].  Relevant  here,  Xu  et  al  [1 1]  resolved  inter-hemispheric  differences  in  the  stiffness  of 

93  mouse  brain  after  complete  occlusion  of  the  middle  cerebral  artery.  They  hypothesized  that 

94  known  formation  of  edema  and  ischemia  due  to  stroke  ipsilateral  to  the  occlusion  and 

95  diaschisis  contralateral  to  the  occlusion  caused  the  observed  inter-hemispheric  differences 

96  in  stiffness. 

97  When  using  small  animal  models  internal  scatter  of  the  shear  wave  can  produce 

98  aliasing  in  stiffness  images  [10-12].  Indeed,  Mace  et  al.  [12]  used  a  research  version  of  the 
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99  SSI  to  remove  this  aliasing  effect  by  tracking  individual  shear  waves  between  scattering 

100  events.  Their  maps  of  brain-tissue  stiffness  displayed  intra-hemispheric  features  in 

101  exquisite  detail. 

102  As  a  first  step  towards  the  goal  of  making  SWE  applicable  to  human  brain,  we  used 

103  SWE  images  from  the  SSI  of  the  brains  of  rats  and  mice  that  received  a  focal  traumatic  brain 

104  injury  to  analyze  the  temporal  evolution  of  the  stiffness  of  brain  associated  with  that  injury. 

105 

106 

107  Materials  and  methods 

108  All  animal  procedures  were  approved  by  the  Institutional  Animal  Care  and  Use 

109  Committee  (IACUC)  of  the  University  of  Washington  and  the  US  Army  Medical  Research  and 

110  Materiel  Command  (USAMRMC)  Animal  Care  and  Use  Review  Office  (ACURO). 

111 

112  Animal  techniques 

113  Surgery 

114  Fifteen  male  Sprague  Dawley  rats  (250-300  grams,  Charles  River)  and  fifteen  male 

115  mitoCFP  mice  (20-30  grams,  bred  at  the  University  of  Washington  [13])  were  used.  These 

116  were  divided  into  three  experimental  groups  for  both  rats  and  mice:  sham  TBI  group  (n=5), 

117  acute  TBI  group  (n=5)  and  24-hour  TBI  group  (n=5).  We  induced  mild  TBI  followingthe 

118  controlled  cortical  impact  (CCI)  protocol  of  Lighthall  [14]  as  refined  by  Dixon  et  al  [15], 

119  which  produced  a  repeatable  and  localizable  cortical  contusion. 

120  Anesthesia  was  induced  via  inhalation  isoflurane  (5%  isoflurane  for  rats;  3%  for 

121  mice)  and  maintained  (2-3%  isoflurane  for  rats;  1-2%  isoflurane  for  mice)  during  surgery 

122  and  cortical  impact.  The  surgical  technique  was  nearly  identical  for  rats  and  mice;  therefore 

123  we  describe  the  procedure  generally  referring  to  the  subjects  as  "rodents." 

124  The  rodent  was  placed  in  a  stereotaxic  device  within  a  pneumatic  controlled  cortical 

125  impact  device  (AmScien  Instruments,  Richmond,  VA,  USA)  to  hold  the  head  steady.  A 

126  heating  pad  was  used  to  maintain  the  core  body  temperature  of  approximately  37°  C.  The 

127  surgical  site  was  prepared  by  shaving  the  top  of  the  head  then  application  of  a  lotion  hair 

128  remover.  Prior  to  surgery,  lidocaine  (1  mg/kg)  was  injected  subcutaneously  to  the  surgery 

129  site.  For  TBI  rodents  allowed  to  survive  24  hours  after  injury,  bupivacaine  (1  mg/kg)  was 

130  also  injected  subcutaneously  to  the  surgery  site  and  buprenorphine  (0.05  mg/kg)  was 

131  delivered  as  an  intraperitoneal  injection. 
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132  After  establishing  an  anesthetic  plane  via  toe  pinch  and  respiratory  observation,  an 

133  incision  was  made  from  the  bregma  to  just  posterior  of  the  lambda  suture  line  and  the  skull 

134  was  exposed.  A  small  (of  order  1  mm)  hole  was  drilled  in  the  right  (ipsilateral)  parietal 

135  portion  of  the  skull  bone  using  a  surgical  drill  (Dremel,  Racine,  WI).  Bone  rongeurs  were 

136  used  to  enlarge  the  drilled  hole  to  approximately  4  mm  in  diameter  between  the  bregma  and 

137  lambda  of  the  right  parietal  bone  of  each  of  mice  and  rats.  (For  mice  this  cranial  window 

138  encompassed  nearly  the  entire  cranial  vault  above  the  ipsilateral  hemisphere  of  brain.  For 

139  rats,  this  hole  took  up  less  than  40%  of  available  ipsilateral  cranial  vault.)  The  cranial 

140  window  was  then  aligned  with  the  tip  of  the  cortical  impact  device.  At  this  point  sham 

141  animals  were  removed  from  the  apparatus.  For  TBI  animals,  the  cortical  impact  device  was 

142  then  fired  to  create  an  impact  at  the  surface  of  the  brain  of  a  velocity  of  3.5  m/s  and  a  depth 

143  of  1.0  mm  for  rats  and  0.5  mm  mice. 

144  For  sham  and  acute  TBI  rats,  the  ipsilateral  cranial  window  was  enlarged  after 

145  impact  to  approximately  the  entire  right  parietal  bone  up  to  the  bregma,  lambda,  and 

146  sagittal  suture  line.  For  mice  we  worked  with  the  existing  ipsilateral  cranial  window.  For 

147  each  of  the  sham  and  acute  TBI  rats  and  mice  an  identical  cranial  window  was  then  created 

148  on  the  contralateral  side  of  the  skull. 

149  For  the  rodents  in  the  24-hour  TBI  groups,  after  impact  the  unilateral  cranial 

150  window  was  filled  with  bone  wax.  The  surgery  site  was  then  closed  with  3-0  nylon  sutures 

151  (Covidien,  Mansfield,  MA,  USA)  and  the  rodents  were  returned  to  their  cages.  After  24 

152  hours,  these  received  bilateral  craniectomies  as  described  above. 

153 

154  SWE  Image  Acquisition 

155  After  the  creation  of  bilateral  cranial  windows  we  imaged  the  brains  of  each  rodent 

156  with  a  linear  SL15-4  MHz  transducer  of  an  Aixplorer  Multiwave  Ultrasound  System 

157  (SuperSonic  Imagine,  Aix-en-Provence,  France).  The  transducer’s  imaging  frequency  ranged 

158  between  4.0-15.0  MHz,  with  a  standard  center  frequency  of  8.5  MHz.  Two  types  of  images 

159  were  collected,  Shear  Wave  Elastography  (SWE)  and  standard  B-mode.  Images  were 

160  collected  with  an  elasticity  range  of  0-180  kPa  to  minimize  gross  imaging  artifacts  associated 

161  with  the  surrounding  tissue  and  skull  while  blinding  the  user  to  intra-  hemispheric  features 

162  at  low  values  of  Young’s  Modulus,  which  generally  he  between  2  and 

163  25  kPa  [11-12,16-18], 

164  A  micro-positioner  was  used  to  maintain  the  SWE  transducer  over  the  rodent’s  head 

165  perpendicular  to  show  coronal  slices.  Multi-purpose  ultrasound  lotion  (Polysonic,  Cone 
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166  Instruments,  Solon,  OH,  USA)  was  used  to  maintain  ultrasound  conduction  between  the 

167  transducer  and  the  head. 

168  The  brain  was  first  imaged  using  traditional  B-mode;  combined  B-mode  and  SWE 

169  imaging  began  when  the  brain  had  been  located  and  centered  in  the  B-mode  image  (Fig.  1 

170  A).  Significant  effort  was  made  to  minimize  creation  of  imaging  artifacts  in  the  SWE  images. 

171  These  artifacts  (e.g.,  Fig.  1  C,D)  are  areas  of  high  Young's  Modulus  that  do  not  correlate  with 

172  the  neuroanatomy  [11-12,16-17].  For  each  rodent,  we  collected  10  -  30  images.  We 

173  selected  3  images  where  both  the  impact  site  was  identifiable  in  the  B-mode  imaging  and 

174  had  minimal  elasticity  artifacts. 

175  After  imaging  the  animals  were  sacrificed  and  perfused  and  the  brains  were 

176  harvested  for  histological  analysis.  The  brains  were  sectioned  through  the  injury  site  at  8- 

177  10  micron  thick  that  were  stained  with  hematoxylin  and  eosin  (H&E)  and,  separately,  cresyl 

178  violet  (CV).  The  slides  were  digitized  and  the  extent  of  TBI  damage  was  quantified  by 

179  manual  selection  of  the  injured  area  using  ImageJ. 

180 

181  Post-processing  of  SWE  Images 

182  Image  post-processing  was  completed  using  MATLAB  (Natick,  MA,  USA).  The  entire 

183  brain  region  was  manually  selected  in  the  SWE-mode  image  obtained  from  the  Aixplorer 

184  (Fig.  1  B).  The  brain  region  in  the  SWE  image  was  selected  based  on  the  structure  of  the  skull 

185  as  seen  in  the  corresponding  B-mode  image. 

186  After  the  selection  of  images,  the  elasticity  map  of  the  selected  brain  region  was 

187  reduced  to  display  an  elasticity  range  of  0-35  kPa  (Fig.  1  C,D),  down  from  0-180kPa  (Fig.  1 

188  A,B).  Regions  of  interest  from  within  each  hemisphere  of  brain  tissue  within  the  skull  were 

189  then  manually  selected  directly  from  the  elasticity  map  (Fig.  1  D).  Regions  of  comparable 

190  size  and  location  but  free  of  imaging  artifacts  were  selected  from  both  hemispheres. 

191  Regional  selection  and  analysis  was  performed  three  times  for  each  image  and  then 

192  averaged  together  to  produce  the  intra-hemispheric  average  of  brain-tissue  stiffness  for 

193  each  of  the  two  hemispheres  of  brain  within  the  image. 

194 

195  Data  Analysis 

196  For  each  image  for  each  rodent,  a  representative  value  of  the  Young’s  Modulus  value 

197  of  each  hemisphere  was  used  directly  for  analysis  as  well  as  normalized  by  dividing  the 

198  average  value  of  Young’s  Modulus  for  a  given  hemisphere  by  the  average  of  the  two 
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199  hemispheres.  The  normalized  Young’s  Modulus  therefore  provides  a  measure  of  the 

200  elasticity  of  a  given  hemisphere  relative  to  the  average  value  across  the  entire  brain. 

201  Kruskal-Wallis  One  Way  Analysis  of  Variance  tests  on  Ranks  (ANOVA]  was  applied 

202  to  the  data  using  SigmaPlot  (Systat  Software,  San  Jose  CA). 

203 

204 

205 

206  Results 

207 

208  Quantification  of  areal  extent  of  damage  in  H&E  and  Creysl  Violet. 

209  We  quantified  the  anatomical  extent  of  TBI  damage  as  shown  histologically  (Fig.  1 

210  E,F}.  There  was  no  difference  between  groups  with  regard  to  the  extent  of  damage  shown 

211  by  H&E  versus  cresyl-violet  stains  nor  in  the  extent  of  damage  between  the  acute  and  24- 

212  hour  groups  for  either  species  (Table  1). 

213 

214  Directly  measured  values  of  Young's  Modulus  for  rats. 

215  Figure  2  shows  directly  measured  average  values  of  Young’s  Modulus  for  each 

216  hemisphere  of  sham,  acute,  and  24-hour  TBI  rats.  There  was  no  significant  difference 

217  between  average  Young’s  Modulus  values  within  any  of  the  sham  or  acute  or  24-hour  TBI 

218  groups  (Fig.  2  A].  We  did,  however,  observe  changes  in  average  Young’s  Modulus  between 

219  groups.  Specifically,  each  of  the  contralateral  and  ipsilateral  hemispheric  values  of  the  acute 

220  and  24-hour  TBI  groups  have  values  of  Young's  Modulus  that  increased  significantly  relative 

221  to  their  respective  values  for  the  sham  TBI  group  (Fig.  2  B). 

222 

223  Normalized  values  of  Young’s  Modulus  for  rats. 

224  Figure  3  shows  normalized  Young’s  Modulus  values  for  each  hemisphere  of  sham, 

225  acute,  and  24-hour  TBI  rats.  There  was  no  significant  difference  between  average  Young’s 

226  Modulus  values  within  the  sham  or  acute  TBI  groups  while  by  24  hours  there  existed  a 

227  significant  difference  between  ipsilateral  and  contralateral  hemispheric  averages  of 

228  normalized  Young’s  Modulus  (Fig.  3  A).  With  regard  to  between-group  changes,  by  24 

229  hours  after  TBI  the  ipsilateral  values  of  normalized  Young’s  Modulus  decreased  with  time 

230  while  the  contralateral  values  increased  with  time  (Fig.  3  B]. 

231 

232 

233 
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234  Directly  measured  values  of  Young's  Modulus  values  for  mice. 

235  Figure  4  shows  directly  measured  values  of  Young's  Modulus  for  each  hemisphere  of 

236  sham,  acute,  and  24-hour  TBI  mice.  There  was  no  significant  difference  between  average 

237  Young’s  Modulus  values  within  the  sham  group  while  there  exist  significant  differences 

238  between  the  average  ipsilateral  and  contralateral  values  of  Young’s  Modulus  both  acutely  and 

239  at  24  hours  after  TBI  (Fig.  4  A).  Relative  to  the  sham  value  of  Young's  Modulus  the  ipsilateral 

240  value  decreased  by  24  hours  after  TBI  while  the  contralateral  value  increased  acutely  before 

241  returning  to  a  value  comparable  to  the  sham  value  of  Young’s  Modulus  (Fig.  4 

242  B). 

243 

244  Normalized  values  of  Young’s  Modulus  for  mice. 

245  Figure  5  shows  normalized  Young’s  Modulus  values  for  each  hemisphere  of  sham, 

246  acute,  and  24-hour  TBI  mice.  There  was  no  significant  difference  between  average  Young’s 

247  Modulus  values  within  the  sham  group.  In  contrast,  the  average  value  of  Young's  Modulus 

248  for  the  ipsilateral  hemisphere  was  significantly  smaller  than  that  of  its  contralateral  value 

249  both  acutely  and  at  24  hours  after  TBI  (Fig.  5  A].  In  addition,  by  24  hours  after  TBI  the 

250  ipsilateral  values  of  normalized  Young’s  Modulus  decreased  with  time  while  the 

251  contralateral  values  increased  with  time  (Fig.  5  B). 

252 

253 

254 

255  Discussion 

256  We  created  mild  traumatic  brain  injury  (TBI]  in  mice  and  rats  using  the  controlled 

257  cortical  impact  (CCI]  paradigm.  We  then  used  a  commercial  diagnostic  ultrasound  machine 

258  to  measure  intra-hemispheric  averages  of  brain-tissue  stiffness  (the  Young’s  modulus,  a 

259  measure  of  shear  resistance  of  tissue)  of  the  brains  both  acutely  and  24  hours  after  TBI 

260  induction.  We  presented  analysis  of  both  the  absolute  value  and  normalized  values  of  brain 

261  stiffness,  the  latter  documenting  changes  in  brain-tissue  stiffness  relative  to  the  average 

262  value  of  the  entire  brain  at  a  given  time  point. 

263  We  note  that  tissue  stiffness  values  inferred  by  ultrasound  shear-wave  elastography 

264  are  sensitive  to  the  tissue’s  net  fluid  content.  For  example,  all  else  being  equal  more  fluid  as 

265  a  percentage  of  tissue  decreases  the  measured  stiffness,  while  less  fluid  increases  the 

266  measured  stiffness  [5].  One  can  appreciate  this  view  by  considering  Derieppe  et  al.  [10], 
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267  whose  study  of  ischemia  in  kidney  showed  a  decrease  in  SWE-derived  stiffness,  and  Xu  et  al. 

268  [11],  whose  study  of  ischemia  in  mouse  brain  showed  a  decrease  in  stiffness  ipsilateral  to 

269  injury  (a  region  of  extensive  edema)  and  an  increase  in  stiffness  contralateral  to  injury  (a 

270  region  of  known  reduced  cerebral  blood  flow).  Paraphrasing  the  literature  analysis  of  Xu  et 

271  al.  [11],  ischemic  brain  within  a  cranium  doesn’t  sag  while  ischemic  kidney  within  a  belly 

272  does. 

273  We  start  with  our  results  for  rats.  The  absolute  value  of  Young’s  modulus  increased 

274  through  time  in  each  of  the  ipsilateral  and  contralateral  hemispheres  compared  to  controls 

275  (Fig.  2).  There  exist  many  physiological  processes  that  occur  after  TBI  that  can  alter  the 

276  fluid  content  of  brain  that  may  explain  our  observations,  subject  to  direct  testing.  For 

277  example,  our  observed  global  increase  in  measured  brain-tissue  stiffness  across 

278  hemispheres  through  time  may  have  arisen  due  the  combined  effect  of  reduced  peripheral 

279  blood  pressure  and  increased  in  ICP,  which  can  work  together  to  reduce  cerebral  perfusion 

280  pressure  (CPP)  hence  cerebral  blood  flow  throughout  the  brain.  This  hypothesis  is 

281  consistent  with  the  observations  of  some  researchers  [19-20]  though  not  all  [21].  Also 

282  contributing  to  this  net  increase  in  brain-tissue  stiffness  may  have  been  a  reduction  in 

283  cerebral  spinal  fluid  (CSF)  within  the  cranium  [20].  We  note  that  these  hypotheses  require 

284  that  global  peripheral  (blood  pressure)  and  intracranial  (ICP;  CPP;  cerebral  blood  flow;  CSF) 

285  processes  more  than  compensate  for  the  focal  increase  in  fluid  content  within  TBI  brain 

286  ipsilateral  to  the  injury  caused  by  the  observed  edema  and  hemorrhage.  That  increase  in 

287  focal  fluid  content  may  itself  balance  in  part  the  focal  ischemia  observed  by  others  after  CCI 

288  [21],  Attention  to  the  normalized  data  highlights  these  potentially  competing,  focal  effects. 

289  Specifically,  we  note  that  by  24  hours  after  TBI  the  normalized  value  of  brain-tissue 

290  stiffness  within  the  ipsilateral  hemisphere  reduced  in  magnitude  relative  to  controls  while 

291  its  counterpart  within  the  contralateral  hemisphere  increased  in  magnitude  relative  to 

292  controls  (Fig.  3).  Moreover,  by  24  hours  after  TBI  the  ipsilateral  value  of  normalized  brain 

293  tissue  stiffness  was  significantly  reduced  relative  to  the  contralateral  value.  Because  of  our 

294  choice  of  normalization  technique,  these  results  document  changes  in  brain-tissue  stiffness 

295  that  occur  relative  to  the  global  rise  in  the  absolute  value  of  brain-tissue  stiffness  reported 

296  above. 

297  With  regard  to  the  observed  reduction  in  relative  ipsilateral  brain-tissue  stiffness, 

298  we  first  note  the  presence  and  time  course  of  edema  and  hemorrhage  demonstrated  by  our 

299  histology  (Fig.  1  E,  F  and  Table  1)  and  consistent  with  other’s  observations  [2  2-2  3].  In 


10 


300  addition  to  edema,  others  have  documented  ischemia  measured  using  laser  Doppler 

301  velocimetry  and  whole-brain  autoradiography  techniques.  Specifically,  Bryan  et  al.  [21] 

302  observed  immediate  ipsilateral  ischemia  after  severe  CCI  that  lasted  at  least  four  hours  after 

303  injury.  We  therefore  hypothesize  that  the  observed  edema  formation  ipsilateral  to  injury  by 

304  24  hours  lowered  the  normalized  value  of  the  measured  intra-hemispheric  average  of  brain 

305  elasticity  compared  to  that  of  sham  TBI  in  a  manner  that  more  than  compensated  for  the 

306  likely  effect  of  ischemia  known  to  occur  in  at  least  severe  CCI,  but  undocumented  within  our 

307  experiment,  all  riding  on  a  wave  of  global  reduction  in  blood  flow  due  to  elevated  ICP.  This 

308  complex  hypothesis,  besides  requiring  verification  in  future  studies,  is  consistent  with  recent 

309  observations  of  reduced  brain-tissue  stiffness  at  the  site  of  CCI  in  mice  relative  to  the  rest  of 

310  the  brain,  based  on  images  created  with  magnetic  resonance  elastography  [24-25]. 

311  With  regard  to  the  observed  increase  in  contralateral  relative  brain-tissue  stiffness, 

312  our  observations  and  those  cited  above  show  no  edema  formation  contralateral  to  the 

313  injury  while  [19],  for  example,  observed  decreased  cerebral  perfusion  pressure  and  reduced 

314  cortical  perfusion  contralateral  to  the  injury,  both  consistent  with  our  observations. 

315  Our  results  for  normalized  values  of  brain-tissue  stiffness  observed  for  mice  after 

316  TBI  (Fig.  5)  are  consistent  with  what  we  observed  for  rats  (Fig.  3).  This  suggests  it  possible 

317  that  the  same  interplay  of  ipsilateral  edema  versus  ischemia  and  global  reduction  in 

318  cerebral  blood  flow  that  can  occur  for  rats  may  also  have  occurred  for  mice.  Consistent  with 

319  this  idea,  Lundblad  et  al.  [26]  observed  after  severe  CCI  significant  edema  in  ipsilateral 

320  cortex,  a  reduction  of  cortical  blood  flow  in  ipsilateral  cortex  on  time  scales  of  hours  to  24 

321  hours  and  in  contralateral  cortex  by  24  hours.  Comparable  studies  are  lacking  for  mild  CCI 

322  and  warrant  future  work. 

323  Our  measurement  of  the  absolute  value  of  intra-hemispheric  brain-tissue  stiffness  of 

324  mouse  brain  differs  significantly,  however,  from  what  we  observed  for  rat  brain.  For 

325  example,  ipsilateral  to  the  injury  site  the  absolute  value  of  Young's  modulus  value  decreased 

326  by  24  hours,  rather  than  increased  as  for  rats.  Also,  the  absolute  value  of  Young’s  modulus 

327  increased  acutely  contralateral  to  the  injury  as  for  rats,  but  returned  to  baseline  values  by 

328  24  hours  after  TBI,  unlike  in  rats,  which  continued  to  increase.  This  difference  may  have 

329  arisen  due  to  the  large  discrepancy  in  areal  extent  between  the  ipsilateral  craniectomy 

330  required  for  mice  versus  that  for  rats.  Specifically,  the  craniectomy  uncovered  essentially 

331  the  entire  hemispheric  cranial  vault  for  mice  while  less  than  half  of  that  vault  for  rats. 

332  Perhaps  this  relatively  large  craniectomy  reduced  the  possibility  of  elevated  ICP,  hence 
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333  reduced  CPP  and  cerebral  blood  flow  in  our  mice.  The  mouse  model  may,  therefore,  only 

334  display  the  effects  of  intra-hemispheric  dynamics  (edema,  local  ischemia)  and  reduced 

335  global  cerebral  blood  flow  due  to  changes  in  systematic  blood  pressure  rather  than  the 

336  additional  effects  of  global  intracranial  dynamics  (elevated  ICP,  hence  a  further  decrease  in 

337  CPP  and  therefore  in  cerebral  blood  flow). 

338 

339  Limitations 

340  This  study  had  several  limitations  due  to  its  preliminary  nature.  First  and  foremost, 

341  we’ve  generated  a  series  of  largely  untested  hypotheses  regarding  the  possible  impact  on 

342  measured  brain-tissue  stiffness  due  to  the  relative  weights  of  different  fluid  sources  within 

343  brain  tissue,  motivated  by  our  direct  observations  and  the  published  literature.  These 

344  hypotheses  require  direct  testing.  Also,  out  of  necessity  we  performed  a  relatively  large 

345  craniectomy  to  facilitate  the  CCI  method  for  mice,  which  may  have  altered  the  effects  of 

346  elevated  ICP  and  reduced  CPP  known  to  occur  in  these  animal  models.  Future  studies 

347  should  target  use  of  larger  animal  models  of  TBI  to  mitigate  this  possible  effect.  In  addition, 

348  we  used  SWE  technology  not  optimized  for  rodent  brain  encased  in  a  skull.  The  artifacts  we 

349  observed  associated  with  the  skull  produced  elasticity  values  that  are  not  consistent  with 

350  other,  published  measurements  of  brain-tissue  stiffness,  forcing  us  to  carefully  select  around 

351  these  artifacts  so  that  they  did  not  confound  our  data. 

352 

353  Conclusion 

354  We  used  ultrasound  shear- wave  elastography  to  measure  absolute  as  well  as 

355  normalized  values  of  intra-hemispheric  stiffness  of  rodent  brain  both  acutely  and  at  24 

356  hours  after  TBI  via  the  controlled  cortical  impact  (CCI)  method.  These  observed  changes  in 

357  brain-tissue  stiffness  after  TBI  may  reflect  multiple  physiological  processes  known  to  occur 

358  within  brain  after  TBI  caused  by  CCI.  These  processes  include  the  initial  formation  of  edema 

359  and  hemorrhage  known  to  develop  ipsilateral  to  TBI  and  observed  by  us,  as  well  as 

360  the  known  but  un-assayed  by  us  (a)  formation  of  ischemia  ipsilateral  to  TBI,  (b)  reduction  in 

361  cerebral  perfusion  pressure  and  (c)  blood  flow  observed  both  ipsilateral  and  contralateral 

362  to  TBI.  Ultrasound-based  elastography  may  therefore  offer  a  method  to  monitor  subtle 

363  changes  in  brain  after  TBI. 

364  Future  work  should  target  testing  the  links  we  have  hypothesized  between 

365  observed  changes  in  brain-tissue  stiffness  and  the  underlying  physiology  of  TBI.  In  addition, 
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366  future  work  should  address  currently  identified  problems  in  the  translation  of  this 

367  technology  to  human  use,  including  delivery  in  a  transcranial  fashion  of  a  shaped  pulse  of 

368  ultrasound  energy  sufficient  to  generate  usefully  observable  shear  waves,  and  the  known 

369  large-scale  anisotropy  of  human  brain  and  its  likely  impact  on  the  formulae  used  to 

370  translate  shear- wave  propagation  into  estimates  of  brain-tissue  stiffness. 

371 
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379  Table  1 


Table  1 

%  extent  of  damage  of  TBI  in  histological  analysis. 

Stain 

Acute  Mice 

24  Hr  Mice 

Acute  Rats 

24  Hr  Rats 

H&E 

8.02%  +/-  2.2% 

14.2%  +/-  7.5% 

4.62%  +/-  2.7% 

7.57%  +/- 1.6% 

CV 

8.33%+/-  0.6% 

11.2%+/- 2.7% 

4.95%+/-  0.1% 

6.98%  +/-  2.1% 
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Figures 


Fig.  1.  Image  Processing.  (A)  Aixplorer  images  in  SWE  mode  (top)  and  B-mode  (bottom)  (B) 
The  brain  region  was  selected  from  the  elastographic  image  (0-180  kPa)  with  reference  to 
the  B-mode  image.  (C)  The  SWE  image  of  the  brain  was  rescaled  to  0-35  kPa  to  enhance  the 
contrast.  (D)  Regions  of  Interest  selected  manually  to  avoid  skull  and  artifacts  (large  values 
of  Young’s  Modulus  greater  than  that  of  brain).  (E)  Hematoxylin  and  Eosin  and  (F)  Cresyl 
Violet  staining  of  corresponding  brain. 


Fig.  3.  Values  of  normalized  Young's  Modulus  for  ipsilateral  and  contralateral  hemispheres 
of  the  brains  of  rats  for  all  test  groups.  *  =  p  <  0.05  between  the  two  indicated  groups. 
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414  Fig.  5.  Values  of  normalized  Young's  Modulus  for  ipsilateral  and  contralateral  hemispheres 

415  of  the  brains  of  mice  for  all  test  groups.  *  =  p  <0.05  bet\veen  the  U.Vo  indicated  groups. 

416 
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Region  Specific  Viscoelastic  Properties  of  the  Adult  Rat  Brain  under  Indentation  following  Traumatic 

Brain  Injury 
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Abstract  Traumatic  Brain  Injury  (TBI)  is  a  serious  health  epidemic  that  places  high  societal  and  economic 
burdens  on  victims  and  their  caregivers.  Further,  the  associated  neuropathological  consequences  that  result 
from  TBI  are  often  complex  and  cause  secondary  injuries  that  are  focal,  diffuse  and  time  dependent.  Current 
computational  models  can  predict  loading  and  deformation  associated  with  TBI;  however,  accurate  knowledge 
of  region  specific  material  properties  from  both  healthy  and  mechanically  damaged  brain  is  needed.  In  this 
study,  the  mechanical  properties  of  both  uninjured  and  traumatically  injured  brain  tissue  are  presented.  Adult 
male  Sprague-Dawley  rats  were  injured  through  a  controlled  cortical  impact  protocol.  Ramp  and  hold 
indentation  tests  were  performed  at  five  locations  on  the  surface  of  tissue  samples  excised  from  whole  brain 
specimens.  Force  displacement  data  were  analyzed  using  quasi-linear  viscoelastic  theory.  An  analysis  revealed 
the  tissue  to  be  viscoelastic  and  spatially  nonlinear  with  mechanical  properties  that  depend  on  both  region  and 
level  of  injury.  After  normalizing  the  data,  the  nonlinear  components  of  the  instantaneous  elastic  force  and 
shear  modulus  were  found  to  be  significantly  lower,  26%,  in  the  region  containing  the  contusion  cavity  on 
severely  injured  samples  compared  to  uninjured  tissue  at  the  same  region  in  controls. 

Keywords  Traumatic  brain  injury,  focal,  diffuse,  controlled  cortical  impact,  mechanical  damage,  quasi-linear 
viscoelasticity 


I.  INTRODUCTION 

Traumatic  Brain  Injury  (TBI)  is  an  important  national  health  concern  in  the  United  States  [1].  Approximately 
1.7  million  incidences  and  52,000  deaths  are  reported  annually  due  to  automobile  collisions,  sports  accidents, 
falls,  and  other  head  impacts  [1].  Further,  severe  TBI  has  a  high  economic  burden,  costing  nearly  $76.5  billion 
each  year  in  medical  and  societal  costs  [1].  On  the  battlefield,  improvised  explosive  devices  have  led  to  TBI  in  as 
many  as  62%  of  soldiers  sustaining  head  injuries,  and  an  estimated  360,000  service  members  have  been 
affected  by  TBI  over  the  past  decade  [2],  [3].  Since  the  1960s  there  has  been  a  dramatic  improvement  in 
understanding  the  complex  pathobiological  behavior  associated  with  TBI  [4].  Brain  injury,  as  a  result  of  an 
impact  or  insult  to  the  head,  leads  to  a  number  of  complex  neuropathological  consequences  that  result  in 
further  tissue  dysfunction  and  eventually  cell  death  [5].  Flowever,  injury  quantification  remains  challenging 
because  the  neurochemical  cascades  that  accompany  TBI  often  involve  complex  secondary  sequelae  that  are 
focal,  diffuse,  and  time  dependent  [6].  A  better  understanding  of  the  mechanical  response  of  the  brain  during 
these  events  would  improve  diagnosis  and  treatment  of  TBI  in  both  clinical  and  battlefield  scenarios. 

The  material  properties  of  brain  during  an  injurious  event  were  studied  as  early  as  the  1940s  [7],  Not  until 
more  recently  has  an  emphasis  been  placed  on  understanding  the  complex  mechanical  behavior  of  the  brain 
during  loadings  that  lead  to  TBI.  Much  of  this  can  be  attributed  to  an  increase  in  awareness  of  the  damaging, 
long-term  effects  following  even  mild  brain  injuries  [l]-[3].  Material  properties  of  brain  are  reported  in  the 
literature  as  viscoelastic  [8]-[13],  both  spatially  and  temporally  non-linear  [9],  [10],  anisotropic  [11],  age- 
dependent  [11]-[14],  inhomogeneous  [12],  [13]  and  nearly  incompressible  [15].  Additionally,  experimental 
factors  such  as  specimen  preparation,  temperature  [16]  and  level  of  hydration  [14],  [17]  have  been  shown  to 
influence  these  properties.  There  is  considerable  variation  in  the  mechanical  response  reported  in  these 
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studies.  Much  of  this  variation  can  be  attributed  to  differences  in  the  type  of  species,  experimental  protocols, 
and  anatomical  regions  tested.  Still,  no  definitive  set  of  material  properties  exist  for  brain. 

Computational  models  of  TBI  are  commonly  used  to  study  the  mechanical  behavior  of  brain  tissue  during  a 
traumatic  event.  Finite  element  modeling  (FEM)  can  predict  loading  of  various  substructures  through 
simulations  of  mild  to  severe  TBI  in  rat  brain  [18].  These  simulations  have  the  added  capability  of  predicting 
injury,  because  model-prediction  of  tissue  strains  were  shown  to  correlate  well  with  experimentally  determined 
strain  and  strain  rate  injury  tolerance  criterion  [19]-[21].  For  more  region-specific  predictions  of  injuries,  local 
material  properties,  including  both  traumatically  injured  and  healthy  brain  tissue,  need  to  be  determined.  Such 
knowledge  would  provide  existing  computational  models  with  the  added  ability  to  predict  the  subsequent 
mechanical  response  of  damaged  tissue  and  allow  for  a  better  understanding  of  brain  excitotoxicity  beyond  that 
of  the  initial  trauma. 

A  survey  of  the  literature  revealed  a  number  of  studies  that  report  thresholds  for  axonal  damage.  The 
reported  mechanical  limits  for  diffuse  axonal  Injury  are  between  0.1-0. 2  Lagrangian  strain  at  strain  rates  greater 
than  10s'1  [9],  [19]-[22].  Flowever,  there  is  little  information  regarding  changes  in  the  mechanical  properties  of 
traumatically  injured  brain  tissue.  Shafieian  et  al.  [23]  used  an  impact  acceleration  model  (1AM)  to  generate 
diffuse  axonal  injury  (DAI)  in  the  brainstem  of  adult  male  Sprague-Dawley  (SD)  rats.  They  reported  a  35% 
reduction  in  the  linear  coefficient  of  the  instantaneous  shear  modulus  in  injured  rats  compared  to  uninjured 
controls.  Saxena  et  al.  [24]  studied  traumatic  spinal  cord  injury  in  adult  female  SD  rats  over  the  course  of  2  and 
8  weeks  post  injury.  They  observed  a  50%  reduction  in  elastic  modulus  for  injured  compared  to  healthy  tissue. 
The  goal  of  this  study  is  to  acquire  accurate,  region-specific  material  properties  for  both  traumatically  injured 
and  uninjured  rat  brain  to  better  understand  the  mechanical  behavior  of  damaged  tissue  and  to  improve 
existing  models  of  TBI. 


II.  METHODS 

Animal  Injury 

All  animal  protocols  were  approved  by  the  University  of  Virginia's  Institutional  Animal  Care  and  Use 
Committee.  Twenty  adult  male  SD  rats  of  average  weight  (mean  ±  SD),  (320±27)gram,  underwent  surgical 
procedures  for  this  study.  Anesthesia  was  induced  with  a  mixture  of  4%  isoflurane  and  100%  medical  grade  02 
for  3-4  minutes  in  an  induction  chamber.  Once  the  animals  were  sedated,  the  level  of  isoflurane  was  reduced 
to  2-2.5%  and  maintained  for  the  duration  of  the  surgery.  The  animals  were  then  placed  in  a  stereotaxic 
reference  device  (MyNeuroLab  Leica  Digital  Stereotaxic  Instrument,  Leica  Biosystems,  Richmond,  IL)  and 
prepared  for  injury.  A  midline  scalp  incision  was  made  along  the  forehead  of  each  animal.  The  skin  and 
underlying  soft  tissue  were  retracted  exposing  the  sagittal,  coronal,  and  lamboid  sutures  on  the  skull.  A  4.5- 
5mm  diameter  hole  was  drilled  from  the  skull  above  the  right  cerebral  cortex  exposing  the  dura  at  the 
coordinates  of  injury:  A  =  -4mm  bregma  and  L  =  2mm  (Ideal  Micro-Drill™,  Flarvard  Apparatus,  Holliston,  MA). 
An  electromagnetically  driven  controlled  cortical  impact  (CCI)  device  (MyNeuroLab  Leica  Impact  One,  Leica 
Biosystems,  Richmond,  IL)  was  used  to  deliver  repeatable,  severe,  open  head  injury  to  nine  of  the  animals;  the 
remaining  animals  were  prepared  for  Sham  injury  and  used  as  controls.  The  2mm  diameter  impact  cylinder  was 
positioned  directly  on  the  dura  at  the  coordinates  of  injury.  Contact  between  the  probe  tip  and  dura  was 
verified  via  an  electric  circuit.  The  probe  tip  was  then  retracted  from  the  dura  and  injury  parameters  were 
inputted  to  the  device.  Severe  traumatic  brain  injury  was  delivered  to  the  right  ipsilateral  cortex  through  a 
rapid  2.5mm  compression  of  the  dura,  sustained  for  200ms,  with  an  initial  impact  velocity  of  3.5m/s.  In  the 
instance  of  Sham  injury,  the  impact  probe  was  retracted  from  the  dura,  but  no  impact  was  performed  on  the 
tissue.  The  animals  were  then  resuscitated  and  monitored  for  a  period  of  twenty-four  hours  after  which  they 
were  sacrificed  and  their  brain  tissue  immediately  collected. 

Sample  Preparation 

Whole  brain  specimens  were  prepared  for  indentation  tests  immediately  following  tissue  collection.  To 
reduce  the  effects  of  temperature  and  level  of  hydration  on  the  results,  specimens  were  submerged  in  a 
physiological  buffer  (Millonig's  Phosphate  Buffer)  for  five  minutes  at  room  temperature  (19-20)°C.  Flydrated 
specimens  were  then  placed  into  a  coronal  slice  matrix  (Braintree  Scientific,  Inc.)  with  incision  planes  spaced 
1mm  apart.  Tissue  cross-sections  were  cut  to  approximately  8mm  in  thickness  from  each  whole  brain  specimen 


using  a  0.23mm  thick  razor  blade  (VWR  Scientific,  Media,  PA).  To  accomplish  this,  two  incisions  were  made  in 
the  coronal  plane  and  parallel  to  each  other.  The  first  incision  was  made  1mm  posteriorly  to  the  injury  plane  at 
-5mm  bregma  and  the  second  8mm  anteriorly  to  the  first  at  3mm  bregma.  In  both  Sham  and  severely  injured 
specimens,  the  injury  plane  was  identified  by  petechial  hemorrhage  on  the  dorsal  surface  of  the  tissue. 
However,  in  the  case  of  a  severe  injury,  the  hemorrhage  was  more  extensive  and  included  a  contusion  cavity. 
Samples  were  then  removed  from  the  slicing  matrix  and  placed  on  an  aluminum  test  stage  with  the  5mm 
bregma  coronal  plane  oriented  upwards.  The  thickness  of  each  sample  deformed  approximately  1mm  under  its 
own  weight.  The  weight  and  dimensions  of  each  sample  were  measured  and  recorded. 

The  coordinates  of  five  regions,  A-E,  on  the  samples  were  determined  using  a  stereotaxic  reference  frame 
[25]  and  are  illustrated  in  Figure  1.  Substructures  of  the  brain  under  the  indenter  at  these  coordinates  included 
both  healthy  and  damaged  tissue  within  the  cerebral  cortex,  corpus  callosum,  hippocampus  and  midbrain.  For 
severely  injured  samples,  region  A  was  located  directly  over  the  contusion  cavity  and  within  the  injured 
hemisphere  of  the  brain.  Samples  were  assumed  to  be  symmetric  about  the  cerebral  fissure,  and  regions  B  and 
D  were  located  contralateral  to  regions  A  and  C,  respectively.  Region  E  was  positioned  at  approximately  the 
center  of  the  sample  on  the  aqueduct.  Coordinates  were  normalized  to  account  for  differences  in  the  cross- 
sectional  dimensions  between  samples  due  to  intra-specimen  variability.  The  normalization  was  performed  by 
making  length  and  width  measurements  on  the  cross-section  of  each  sample;  the  five  indentation  coordinates 
were  then  multiplied  by  the  ratio  of  the  cross-sectional  dimensions  of  the  first  sample  tested  under  this 
protocol,  ID  NIB00287,  to  those  of  the  sample  of  interest;  the  normalized  coordinate  locations  were  then 
dimensioned  from  region  E  using  digital  micrometers.  Evan's  Blue  Dye  was  used  to  mark  each  coordinate  on 
the  tissue  cross-section  for  a  visual  reference  onto  which  the  indenter  could  be  positioned.  The  amount  of  dye 
under  the  indenter  was  assumed  to  have  a  negligible  effect  on  the  tissue  properties.  A  total  of  30  minutes  were 
allotted  for  sample  preparation,  i.e.  from  the  time  of  tissue  collection  to  the  time  indentation  testing  began. 


Figure  1:  Schematic  of  a  tissue  sample  showing  the  coordinate  locations  of  the  five  un-normalized  indentation 
regions  A  =  (2,  2.5),  B  =  (-2,  2.5),  C  =  (3.5,  -1),  D  =  (-3.5,  -1)  and  E  =  (0,  0).  Dashed  circles  represent 

the  cross-sectional  area  of  the  indenter  overlaid  onto  the  tissue  sample. 

Indentation  Testing 

The  aluminum  test  stage  with  sample  was  mounted  atop  a  50gram  load  cell  (Model  31  Low,  Honeywell 
International  Inc.,  Golden  Valley,  MN),  and  beneath  a  3.18mm  diameter  plane-ended  cylindrical  indenter 
mounted  to  a  linear  actuator  equipped  with  an  LVDT  to  measure  displacement  (ElectroForce®  3100  Test 
Instrument,  Bose  Corporation  -  ElectroForce  Systems  Group,  Eden  Prairie,  MN).  Excess  compliance  in  the  test 
frame  due  to  the  motion  of  the  actuator  induced  an  inertial  based  force  response  in  the  load  cell.  A  500g  linear 
accelerometer  (Model#:  7264B-500,  Humanetics  Innovative  Solutions,  Plymouth,  Ml)  was  mounted  to  the  test 
stage  to  subtract  off  this  effect.  Force,  displacement,  and  acceleration  data  were  acquired  at  20kHz  (DEWE- 


2010,  Dewetron  Inc.,  Wakefield,  Rl).  Regions  were  tested  in  a  randomized  order  for  each  experiment.  At  each 
coordinate  the  indenter  was  centered  on  the  Evan's  Blue  Dye.  This  was  accomplished  by  mounting  a  1.3mm 
diameter  spherical  tip  punch  to  the  actuator  and  positing  directly  over  the  dye  via  visual  inspection.  The 
spherical  punch  was  exchanged  with  the  plane-ended  cylindrical  indenter,  which  was  assumed  to  be  centered 
over  the  indentation  coordinate.  The  indenter  was  advanced  slowly  toward  the  tissue  at  a  rate  of  O.Olmm/s 
until  a  tare  load  of  0.3gram  was  achieved.  The  indenter  tip  was  then  pressed  0.6mm  into  the  tissue,  normal 
with  respect  to  the  local  surface,  in  approximately  8ms  and  then  held  for  30s  to  measure  the  tissue's  relaxation. 
A  peak  displacement  velocity  of  120mm/s,  at  an  approximate  strain  rate  of  17s"1,  was  observed  during  the  ramp 
portion  of  the  displacement.  After  each  test  the  indenter  was  carefully  removed  from  the  tissue.  Five  minutes 
were  allotted  between  tests  to  allow  for  tissue  recovery  and  instrumentation  adjustments  [26].  The  tissue  was 
sprayed  with  Millonig's  in  between  each  indentation  test.  The  protocol  was  repeated  for  the  remaining  four 
regions  and  all  testing  was  completed  within  75  minutes  of  animal  sacrifice. 


Mathematical  Modeling 

All  data  were  filtered  in  accordance  to  the  SAE-J211  standard,  CFC  1000,  using  a  zero-phase,  digital  MR  8  pole 
butterworth  filter  at  a  Low  Pass  frequency  of  1650Flz.  The  data  were  resampled  in  a  logarithmically  scaled  time 
step  to  give  equal  weights  to  both  ramp  and  hold  portions  of  the  test.  Samples  were  assumed  to  be 
incompressible  and  isotropic  [9],  [10].  The  force  response,  to  the  displacement  input,  □  □  □ff.y'T, 

was  modeled  using  a  quasi-linear  viscoelastic  (QLV)  mathematical  framework  [27]: 

□  avoftyn-i-t! - (1) 
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where  is  the  instantaneous  elastic  response,  is  the  reduced  relaxation  function,  y  is  the  time, 

and  t  is  a  dummy  variable  over  which  the  convolution  integral  (1)  is  evaluated.  The  instantaneous  elastic 
response  was  modeled  using  the  solution  to  the  Boussinesq  problem  for  a  flat-ended  cylindrical  punch  [28]: 
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where  □  is  the  radius  of  the  indenter,  0  is  Poisson's  ratio,  which  was  assumed  to  be  0.5,  0  is  the  shear  modulus, 
and  0  is  a  constant  used  to  incorporate  the  effect  of  substrate  on  finite  sample  thickness  [29].  Values  of  0  were 
obtained  for  each  sample  and  found  to  be  between  1.26  and  1.3.  The  shear  modulus  was  chosen  to  be  a 
second-order,  even  function  of  □  to  capture  the  spatial  nonlinearity  of  the  tissue  [23]: 

0  □  mnm  00  □  (3) 


where  00  and  0^,  are  the  instantaneous  linear  and  nonlinear  shear  modulus  coefficients,  respectively.  Using  this 
form  for  the  shear  modulus  results  in 
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where  and  are  the  linear  and  nonlinear  coefficients  of  J(tn<Y,  respectively,  described  in  (2). 
The  mathematical  solutions  for  the  values  of  00  and  0^,  in  terms  of  and  ^  are  determined  through  the 
use  of  equations  (2-4). 
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A  six  term  prony  series  with  five  time  constants  was  chosen  to  model  the  relaxation  behavior  of  the  tissue: 

a  a  2 

0  ^  under  the  constraint  that  0  ^  Q  1  (6) 
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where  x>0's  are  the  normalized  relaxation  coefficients  of  the  corresponding  time  decades  and  -o^is  the 
coefficient  of  the  steady-state  response.  Values  for  the  thirteen  coefficients  00,  x>0,  for  D=1  to  5  and 

^injwere  determined  through  a  reduced  gradient  algorithm  (Excel  Solver®,  Microsoft®,  Redmond,  WA)  that 
was  used  to  minimize  the  sum  squared  error  between  the  model-predicted  force,  resulting  from 
numerical  integration  of  (1),  and  the  experimental  data.  An  individual  set  of  optimal  coefficients  was 
determined  for  each  indentation  test.  Preliminary  analysis  of  the  model  fit  to  the  first  few  data  sets;  test  ID 
NIB00288A  through  NIB00288E,  indicated  marginal  variability  in  the  values  for  the  optimized  time  constants, 
EI^O.OOls,  m^o.ois, 

0g=O.ls,  0g=ls,  02=lOs.  To  simplify  the  model  the  time  constants  were  fixed  at  theses  decades  for 
the 

remainder  of  the  analysis  and  only  eight  parameters  needed  to  be  optimized  through  model  fitting. 

Statistics 

The  number  of  terms,  □,  in  equation  (6)  were  determined  via  an  F-test  [30].  Data  from  test  ID  NIB00318E  was 
modeled  with  four,  five,  and  six  time  constants  to  see  if  there  was  a  statistically  significant  improvement  in  the 
model's  fit  to  the  data.  The  model  with  five  time  constants  gave  a  significantly  better  fit  (F=9,  p  <  0.001)  than 
the  model  with  four  time  constants  and  the  model  with  six  time  constants  showed  no  improvement  in  fit  over 
five  time  constants  (F=0,  p=l).  Therefore  the  model  with  five  time  constants  was  chosen  for  the  analysis.  The 
critical  F-statistic  at  the  0=0.05  level  of  significance  for  both  tests  was  Fc  =1.  A  total  of  n=8  indentation  tests 
were  performed  per  region,  A-E,  and  per  injury  treatment,  Sham  and  severe  injury,  for  a  total  of  10  groups.  For 
each  group,  an  average  and  .off-V'T  were  determined  using  least  squares  optimization  between  the 

average  and  the  eight  individual  measurements.  The  coefficients  of  the  shear  modulus,  00  and  0a  were 
calculated  from  and  of  the  eight  individual  curves  using  expression  (5).  An  average  00  and  0a  were  then 
determined  for 

each  group.  Additionally,  the  coefficients  of  and  at  regions  A  and  C  were  normalized  to  the 

values 

at  the  contralateral  regions  B  and  D,  respectively.  Specifically,  normalization  was  performed  by  dividing  the 
value  of  a  particular  absolute  coefficient  00  and  0a  at  region  A  by  its  corresponding  absolute 

contralateral  value  at  region  B.  For  example,  from  indentation  test  ID  NIB00290A  was  divided  by  from 
indentation  test  ID  NIB00290B,  etc.  The  process  was  repeated  for  the  coefficients  at  region  C,  dividing  by  the 
corresponding  contralateral  values  at  region  D.  The  absolute  structural  and  material  properties  ( 00,  0a, 
om,  for  |Ji=1..5  and  were  compared  separately  to  evaluate  for  the  effect  of  region  and  injury  treatment  using 
a  two-way  ANOVA.  Post-hoc  comparisons  were  made  using  a  student's  t-test  with  the  appropriate  Bonferroni 
correction.  Samples  were  assumed  to  be  independent  measurements  of  a  particular  tissue  property,  normally 
distributed  and  homoscedastic.  To  evaluate  the  differences  observed  in  the  normalized  ratios,  student's  t-tests 
with  a  Bonferroni  correction  were  used  to  make  comparisons  across  injury  treatment  at  a  particular  region. 
Specifically,  a  direct  comparison  was  made  between  the  normalized  ratios  of  Sham  and  severely  injured  tissue 
at  region  A  to  evaluate  the  effect  of  injury  at  the  location  of  the  contusion  cavity. 

III.  RESULTS 

Eighty  indentation  tests  were  performed  on  the  tissue  samples,  and  data  from  16  of  the  original  20  animals 
were  used  in  analysis.  The  tissue  from  animal  ID  NIB00311  was  damaged  upon  collection  and  the  data  from  3 
additional  experiments  were  compromised  due  to  either  unpreventable  noise  from  the  surroundings  (animal  ID 
NIB00287),  or  experimental  error  (animal  IDs  NIB00295  and  NIB00298).  Average  and  .off.y'T  for  region 

and  injury  treatment  are  shown  in  Figure  2.  Region  specific,  absolute  structural  and  material  properties, 
normalized  ratios,  and  results  of  the  t-tests  are  reported  in  Table  1  (see  appendix).  ANOVA  indicated  significant 
(p  <  0.05)  main  effects  of  both  region  and  injury  treatment  on  the  absolute  coefficients.  On  average,  ^  and 
00  were  higher  (p  <  0=0.005)  in  region  E  compared  to  region  A  ( +36.3mN/mm,  p=0.002  and  00:  +2.23kPa, 
p=0.002).The  relaxation  coefficient,  sx$,  was  found  to  be  higher  while  and  flirowere  lower  (p  <  0=0.05)  in 
severely  injured  samples  compared  to  Sham  controls  +0.011,  p=0.0164,  x>g:  -0.003,  p=0.008,  and 
-0.003,  p=0.044).  The  student's  t-test,  revealed  a  significant  decrease  (p  <  0=0.025)  in  the  value  of  the 
normalized  ratios  of  and  0a,  (-26%,  p=0.0084  each),  in  severely  injured  tissue  compared  to  Sham  controls. 
Conversely,  the  t-tests  revealed  a  significant  increase  in  the  normalized  ratios  of  and  0a ,  (+38%, 
p=0.0156  each)  on  severely  injured  compared  to  sham  controls  in  region  C.  The  percentages  reported  here 
are  calculated  as  percent  differences  in  the  sample  means  from  Sham  samples. 
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Figure  2:  Average  ?£.□??  and  .ofty7?  for  region  and  injury  treatment.  Error  bars  are  95%  Confidence 
Intervals.  Average  _S?(T.D<T  curves  for  severe  injury  were  on  average  stiffer  than  the  average  Sham  curves  in 
all  regions  except  for  region  A.  The  reduced  relaxation  functions  were  nearly  identical  in  all  cases. 


IV.  DISCUSSION 

This  study  found  the  mechanical  properties  of  the  rat  brain  to  be  viscoelastic,  spatially  nonlinear,  and 
dependent  on  both  region  and  injury  treatment.  A  linear  viscoelastic  model  was  fit  to  the  experimental  data,  in 
addition  to  the  QLV  model,  where  equation  (3)  was  assumed  to  be  a  function  of  only.  With  the  same 
number  of  prony  series  terms,  the  fits  of  both  the  linear  and  QLV  models  were  visually  assessed  using  tissue 
force  time  histories,  and  compared  statistically  using  the  F-test,  Figure  3.  The  linear  model  fit  the  data  well 
during  the  first  of  half  of  the  ramp  and  long  term  relaxation  of  the  tissue;  however,  the  addition  of  the 
nonlinear  term  through  QLV  showed  a  statistically  significant  improvement  in  the  model  fit  (p  <  0.001)  to  the 
experimental  data,  capturing  the  entire  ramp,  peak  force,  and  initial  tissue  relaxation,  Figure  3a.  To  justify  the 
use  of  QLV  over  a  fully  nonlinear  viscoelastic  model  the  ratio  of  the  relaxation  forces  from  two  different 
displacement  steps  were  calculated  and  then  checked  to  be  approximately  constant  in  time  using  linear 
regression.  Two  displacement  steps,  the  first  to  0.6mm  and  the  second  to  1.2mm,  were  applied  to  tissue  from  a 
severely  injured  animal;  ID  NIB00291  at  region  A  and  region  E.  Results  from  the  linear  regression  were  used  to 
evaluate  whether  or  not  the  slope  of  the  force  ratio  was  statistically  significant  from  zero.  Data  up  to  100ms 
after  the  peak  force  were  not  included  in  the  analysis  due  to  transience  of  the  displacement  ramp.  The  value  of 
the  regression  coefficient,  the  slope,  was  found  to  be  statistically  significant  (slope  =  0.035s1,  p  <  0.001)  at 
region  A  and  (slope  =  0.007s'1,  p  <  0.001)  at  region  E.  However,  the  magnitude  of  the  slope  was  not  thought  to 
be  meaningful.  That  is,  the  ability  to  detect  small  changes  in  the  value  of  the  slope  was  due  to  the  large  amount 
of  data  being  used,  and  that  the  resulting  values  were  not  influential,  suggesting  a  relatively  constant  response 
over  time.  These  results  indicated  that  the  relaxation  behavior  of  the  tissue  was  independent  of  displacement 
and  that  no  temporal  nonlinearities  were  observed  up  to  approximately  18%  sample  penetration. 
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Figure  3:  Examples  of  the  QLV  model  fit  to  the  experimental  data  for  both  the  ramp  (a)  and  hold  (b)  portions  of 
an  arbitrarily  chosen  indentation  test  (NIB00308B).  The  QLV  model  followed  the  experimental  data  more 
closely  than  the  linear  model  during  the  ramp,  peak,  and  initial  relaxation  of  the  tissue  (a).  Both  models 
followed  the  data  closely  during  the  hold  portion  of  the  test  (b).  A  similar  result  was  observed  when  modeling 
other  experimental  data.  The  increase  in  force  near  16ms  and  30ms  is  due  to  the  increase  in  displacement  of 
the  indenter  into  the  tissue  at  these  times. 

Validation  of  the  test  methodology  selected  for  the  current  study  involves  a  comparison  to  other  brain 
studies  reported  in  the  literature.  The  absolute  mechanical  properties  and  uncertainties  found  in  the  current 
study  are  consistent  with  those  determined  experimentally  for  shear  [9],  [31],  compression  [32]  and  indentation 
[23].  Further,  the  data  fit  within  the  reported  range  of  shear  moduli  for  brain  tissue  (0.1-22)  kPa  [8]-[16],  [23] 
and  [32]-[34].  Darvish  and  Crandall  characterized  bovine  brain  using  QLV  under  oscillatory  shear  tests  up  to 
200Flz  and  20%  Lagrangian  strain  [9].  They  observed  16kPa  and  2.62kPa  for  the  linear  and  nonlinear 
instantaneous  elastic  shear  moduli,  respectively.  On  the  other  hand,  Takhounts  et  al.  tested  both  bovine  and 
human  brain  up  to  100%  Lagrangian  shear  strain  in  ramp  and  hold  tests  and  determined  the  instantaneous 
linear  shear  moduli  to  be  approximately  2kPa  and  1.5kPa,  respectively  [10].  Shuck  and  Advani  found  large 
variability  in  the  shear  modulus  [31].  They  performed  oscillatory  torsion  tests  on  human  brain  up  to  60Hz  and 
observed  (3-16)kPa.  In  each  of  these  studies  the  brain  was  modeled  under  the  assumptions  of  isotropy  and 
incompressibility. 

In  regards  to  the  results  reported  for  indentation  studies,  Gefen  et  al.  tested  both  young  and  mature  rats  in 
vivo  and  in  situ  under  spherical  indentation,  and  determined  values  for  the  instantaneous  shear  modulus 
between  (1.2-3.3)kPa  [15].  Further,  Gefen  and  Margulies  [33]  compared  the  effects  of  in  vitro,  to  in  vivo  and  in 
situ  on  the  material  properties  of  porcine  brain.  They  used  a  hemispherical  indenter  and  found  the  short  term 
shear  modulus  between  (0.7-2.8)kPa.  In  both  studies  preconditioning  was  found  to  significantly  reduce  the 
value  of  the  shear  modulus  and  in  vitro  results  tended  to  be  lower  than  in  vivo  and  in  situ.  Samples  were  not 
preconditioned  in  the  current  study  as  work  by  Shafieian  et  al.  [23]  suggests  preconditioning  may  have 
damaged  their  samples  and  reduced  the  effect  of  injury.  Various  other  studies  have  used  micro-indentation  to 
determine  region  specific  material  properties  in  rat  brain  [13],  [34].  These  studies  observed  values  for  the  short 
term  shear  modulus  on  the  order  of  (0.1-1.5)kPa,  nearly  one  order  of  magnitude  lower  than  the  values  reported 
in  the  current  study.  The  material  properties  of  brain  have  shown  to  be  rate  sensitive  and  increase  with 
increasing  strain  rate  or  frequency  [9]  [10].  The  load  rates  reported  in  these  studies  were  less  than  lmm/s, 
(0.4-0. 5)s'\  and  much  less  than  that  of  the  current  study  which  could  explain  this  discrepancy.  Additionally, 
these  studies  performed  indentation  tests  up  to  depths  of  40pm  on  the  surface  of  tissue  cross  sections  that  had 
previously  been  blocked  and  then  mechanically  cut  using  a  vibratome.  Consequently,  their  findings  may  have 
been  significantly  altered  from  that  of  normal,  healthy  tissue.  The  tissue  samples  tested  in  the  current  study 
were  hand  cut  and  tested  in  vitro.  As  a  result  of  sample  preparation,  a  thin  layer  on  the  surface  of  the  tissue 
likely  sustained  damage.  Even  though  the  indentation  depths  used  in  the  current  study  were  much  higher  than 
those  used  in  micro-indentation,  it  is  reasonable  to  expect  some  alterations  in  the  mechanical  properties  from 
that  of  normal,  healthy,  living  tissue. 

In  a  similar  study  to  the  current,  Shafieian  et  al.  performed  cylindrical  indentation  at  two  locations,  PDx  and 
PmJ,  on  the  brainstem  of  impacted  rats  and  found  a  statistically  significant  reduction  in 
between  uninjured  and  injured  specimens  [23].  They  reported  values  for  00  and  0^  between  (l-10)kPa 
and  (1- 

25)kPa/mm2,  respectively  with  similar  uncertainties  to  that  of  the  current  study.  The  reported  reduction  in  00 
between  injured  specimens  and  Sham  controls  was  35%.  This  value  was  an  average  taken  from  both 
preconditioned  and  un-preconditioned  samples  and  found  to  be  in  agreement  with  the  work  of  Darvish  and 
Crandall  [12]  who  reported  a  33%  reduction  in  the  linear  shear  modulus  after  non-recoverable  strain 
conditioning.  Flowever,  Shafieian  et  al.  did  not  specify  the  significance  of  change  in  00  and  0^,  with  respect  to 
injury  treatment,  making  it  unclear  which  parameter,  if  any,  was  driving  the  reduction  in  shear  modulus  after 
tissue  damage  [23].  On  average,  they  saw  a  28%  and  47%  reduction  in  00  and  0$,  respectively,  for  un¬ 
preconditioned  samples  at  PDx.  Further,  they  observed  a  14%  and  29%  reduction  in  00  and  0a  for  un¬ 
preconditioned  samples  at  PmJ;  however,  this  result  was  not  statistically  significant.  In  the  current  study  00  and 
0^1  were  reduced  on  average  17%  and  26%,  respectively  at  the  location  of  the  contusion  cavity,  region  A.  This 
was  observed  after  normalizing  the  data  and  the  reduction  in  00  was  not  statistically  significant.  In  the  case  of 


the  absolute  material  properties,  the  observed  differences  in  00  and  between  severely  injured  and  Sham 
samples  were  not  significant.  Additionally,  the  reported  38%  increase  in  0^  at  region  C  on  severely  injured 
compared  to  sham  tissue  contradicts  the  findings  for  the  normalized  properties  at  region  A.  Further 
investigation  into  this  matter  is  needed;  however  there  was  no  evidence  of  mechanical  damage  in  this  region 
when  compared  to  the  visible  damage  of  the  contusion  cavity  in  region  A. 

Possible  reasons  for  these  discrepancies  can  be  attributed  to  a  number  of  factors.  In  the  current  study,  a  CCI 
model  was  used  to  deliver  injury  directly  to  the  cortex  [35].  On  the  other  hand,  Shafieian  et  al.  chose  an  impact 
acceleration  model  to  generate  DAI  in  the  brainstem  due  to  its  predictable  pattern  of  injury  there  [23],  [36]. 
Further,  regions  tested  in  the  current  study  incorporate  both  white  and  grey  matter  and  are  heterogeneous 
compared  to  the  brainstem  which  is  comprised  of  predominantly  white  matter  and  is  relatively  homogenous 
and  stiffer  than  the  cortex  [37],  Another  explanation  for  these  discrepancies  is  the  rate  at  which  the  tissue 
samples  were  loaded  under  indentation.  Peak  Loading  rates  determined  from  the  tests  conducted  by  Shafieian 
et  al.  [23]  were  between  (16-33)mm/s,  (8-16)s"\  while  that  of  the  current  study  was  approximately  120mm/s, 
17s1.  High  rate  inputs  are  necessary  to  characterize  the  mechanical  response  of  the  brain  during  TBI  [38].  The 
peak  loading  rates  reported  in  the  current  study  are  within  the  range  of  those  reported  for  impact  traumas  [13], 
[38],  and  [39].  An  understanding  of  the  brain's  mechanical  response  at  high  load  rates  would  be  useful  for 
computational  models  of  these  events.  However,  testing  at  these  rates  may  cause  further  damage  to  the  tissue 
and  reduce  the  effect  of  injury  on  the  mechanical  properties.  Strain  and  strain  rate  tissue  tolerance  thresholds 
have  previously  been  studied  and  the  reported  values  associated  with  axonal  injury  are  between  0. 1-0.2 
Lagrangian  strain  with  strain  rates  greater  than  10s1  [9],  [19]-[22]. 

The  hypothesis  that  the  tissue  samples  were  damaged  upon  loading  was  examined  in  a  separate  analysis.  A 
QLV  model  with  only  1  time  constant  @9=  0.001s  was  fit  to  the  ramp  portion  the  experimental  data  (n=16)  for 
region  A.  The  form  of  was  kept  the  same  as  expression  (4).  The  model  was  fit  in  three  ways.  The 

sum 

squared  error  was  minimized  between  the  model  and  experimental  forces  from:  (1)  the  toe  region  to  the  force 
corresponding  to  the  peak  loading  rate;  (2)  the  force  corresponding  to  the  peak  loading  rate  up  to  the  peak 
force;  and  (3)  over  the  entire  ramp  from  the  toe  region  to  the  peak  force.  Average  Gl^and  were  determined 
from  the  (n=16)  individual  fits.  ANOVA  revealed  significant  (p  <  0.05)  differences  between  the  three  model  fits. 
Post-hoc  Bonferroni  comparisons  were  made  between  individual  samples  and  revealed  the  following 
information.  The  value  of  00  was  approximately  constant  across  all  three  models,  while  the  value  of  0^ 
decreased  by  34%  (p  <  0.001)  after  the  peak  displacement  rate  of  120mm/s.  This  observation  was  made 
independent  of  injury  treatment.  The  stability  of  the  model  was  evaluated  during  each  of  the  three  fits  due  to 
the  relatively  small  amount  of  data  being  fitted.  Regardless  of  initial  inputs,  model  parameters  optimized  to  the 
same  values,  indicating  a  stable  solution.  These  results  suggest  that  the  tissue  may  have  been  damaged  during 
loading  up  to  5%  tissue  penetration  and  at  a  peak  rate  of  17s"1. 

V.  CONCLUSIONS 

This  study  presents  an  experimental  methodology  and  analytical  framework  for  modeling  the  region  specific 
structural  and  material  properties  of  mechanically  injured  and  uninjured  brain  tissue.  Force  data  were  acquired 
under  high  rate  loading  inputs  in  the  range  of  those  related  to  impact  traumas.  Quasi-linear  viscoelasticity,  a 
popular,  constitutive  model  that  is  commonly  used  to  model  soft  biological  materials  was  chosen  to  fit  the 
experimental  data  over  both  linear  and  fully  nonlinear  viscoelastic  models.  The  tissue  was  assumed  to  be 
isotropic  and  incompressible  for  model  simplification.  The  mechanical  properties  were  found  to  be  viscoelastic, 
nonlinear,  and  regionally  dependent.  Additionally,  these  properties  were  quantified  under  a  well  characterized 
injury  model  capable  of  delivering  repeatable  levels  of  mechanical  damage  directly  to  the  dura.  To  the  authors' 
knowledge,  this  is  the  first  study  to  examine  the  mechanical  properties  of  the  brain  after  CCI.  The  results  for 
the  shear  modulus  were  within  the  range  reported  in  the  literature  for  shear,  compression,  and  indentation 
tests.  The  material  and  structural  properties  were  found  to  be  roughly  constant  across  region  and  injury 
treatment.  However,  after  normalizing  the  region  specific  mechanical  properties  in  regions  A  and  C,  there  were 
statistically  significant  differences  in  the  values  of  the  nonlinear  coefficients  of  the  instantaneous  shear  and 
elastic  response.  A  26%  reduction  in  the  nonlinear  material  and  structural  coefficients  were  observed  in 
severely  injured  samples  compared  to  healthy  controls  at  the  location  of  the  contusion  cavity.  The  utility  of  this 
research  is  crucial  for  understanding  the  mechanical  response  of  the  brain  after  TBI.  Knowledge  of  such 


material  properties  may  be  useful  to  uniquely  identify  different  types  of  brain  injury  and  to  better  understand 
the  mechanics  of  repetitive  brain  injuries. 


VI.  LIMITATIONS 

There  are  a  number  of  limitations  that  may  influence  overall  response  of  the  tissue  samples  during 
mechanical  loading.  First,  the  solution  to  the  Boussinesq  problem  assumes  indentation  of  a  semi-infinite  elastic 
half-space.  The  samples  used  in  the  current  study  have  finite  boundaries  and  exhibit  damping  as  well  as  elastic 
properties.  This  solution  was  modified  by  [29]  to  incorporate  the  effect  of  a  rigid  substrate  on  finite  sample 
thickness.  The  corrective  factor  13  was  developed  for  indentation  on  articular  cartilage  and  subsequently 
applied  to  indentation  tests  performed  on  other  soft  biological  tissues  including  brain  [13],  [23].  The  factor  13 
functions  to  reduce  the  applied  load,  measured  within  the  tissue,  by  an  amount  that  depends  on  the  ratio  of  the 
indenter  radius  to  sample  height.  As  the  indenter  radius  increases  or  sample  height  decreases,  higher  forces 
are  transmitted  from  the  substrate  to  the  tissue.  To  avoid  the  effect  of  substrate  on  mechanical  properties,  a 
common  rule  of  thumb  is  to  limit  indenter  penetration  to  depths  of  no  greater  than  10%  of  the  total  sample 
thickness  [40].  Penetration  depths  in  the  current  study  obey  this  rule.  However,  more  recent  work  has 
observed  noticeable  substrate  effects  within  10%  penetration  and  claim  that  the  ratio  of  the  indenter  radius  to 
sample  height  must  be  kept  within  10%  [41].  The  ratio  of  radius  to  sample  height  is  20%  in  the  current  study 
and  13  was  implemented  to  adjust  for  this  limitation.  In  addition  to  substrate  effects,  other  boundary  conditions 
have  been  violated.  In  all  regions  except  E,  the  distance  between  the  indenter  and  edge  of  the  tissue  sample 
was  less  than  the  recommended  distance  of  one  indenter's  width.  This  would  have  the  effect  of  reducing  the 
material  properties  observed  in  regions  A-D  compared  to  E.  Indenter  interference  is  an  additional  concern  in 
this  study.  The  process  by  which  the  coordinates  were  dimensioned  on  the  surface  of  the  tissue  was  imperfect. 
Slight  offsets  in  the  dimensioning  may  have  led  to  overlapping  indentation  tests.  Further,  the  indentations 
made  on  the  surface  of  the  samples  were  within  one  diameter  of  each  other.  However,  adequate  time  was 
allotted  for  tissue  recovery  in  between  subsequent  tests.  Finally,  contusion  cavities  can  be  problematic  for 
material  studies  [23].  They  are  often  near  the  edge  of  the  sample  and  create  discontinuities  in  the  tissue 
surface  making  it  difficult  to  obey  boundary  conditions  under  most  mechanical  tests.  It  is  suggested  that  future 
studies  involving  contusion  cavities  be  made  in  vivo  or  situ  with  the  brain  left  in  the  skull. 

A  number  of  studies  including  the  current  have  used  the  F-test  to  statistically  determine  an  appropriate 
model  for  fitting  experimental  brain  data  [13],  [23],  [34];  however,  there  are  limitations.  The  F-test  assumes  the 
data  to  be  independently  measured  and  normally  distributed.  In  addition,  there  must  be  a  linear  relationship 
between  independent  and  dependent  variables  in  the  models  chosen  to  fit  the  data.  The  force  data  in  the 
current  study  are  not  normally  distributed  and  not  independent  measurements  in  time.  Further,  in  both  the 
linear  viscoelastic  and  QLV  models  the  relationship  between  force  and  time  is  not  linear.  Therefore  it  is  not 
clear  that  hypothesis  testing  using  the  F-test  is  suitable  or  meaningful  to  determine  the  most  appropriate 
model.  In  this  case,  visual  inspection  of  the  model  fit  to  the  data  may  provide  a  better  indicator  of  goodness  of 
fit. 
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IX.  Appendix 


Table  1:  Absolute  Structural  and  Material  Coefficients  of  Jl?tD<y,  and  IZlfcDjT  and  Normalized  Ratios  for  Region  and  Injury 

Treatment. 

Absolute _  _ A _  _ B _  _ C _ 

Unit  Coef.  Sham  Severe  Injury  Sham  Severe  Injury  Sham  Severe  Injury 


mN/mm 

Fi 

140.8 

± 

19.2 

128.2 

± 

33.0 

141.1 

± 

24.1 

155.2 

± 

44.7 

143.4 

± 

30.7 

179.0 

± 

38.7 

mN/mm3 

f3 

111.6 

± 

24.1 

96.6 

± 

19.7 

97.3 

± 

17.4 

115.6 

± 

22.3 

110.8 

± 

23.6 

138.1 

± 

16.4 

kPa 

Po 

8.66 

± 

1.23 

7.89 

± 

2.04 

8.68 

± 

1.51 

9.55 

± 

2.76 

8.85 

± 

2.00 

11.02 

± 

2.40 

kPa/mm2 

P2 

6.88 

± 

1.57 

5.95 

± 

1.23 

5.99 

± 

1.12 

7.11 

± 

1.38 

6.83 

± 

1.53 

8.50 

± 

1.03 

Gi  0.739  ±  0.021  0.763  ±  0.010  0.749  ±  0.021  0.752  ±  0.015  0.746  ±  0.018 

G2  0.122  ±  0.007  0.117  ±  0.007  0.121  ±  0.006  0.123  ±  0.004  0.122  ±  0.006 

G3  0.057  ±  0.005  0.053  ±  0.004  0.055  ±  0.002  0.053  ±  0.006  0.057  ±  0.003 

G4  0.029  ±  0.005  0.026  ±  0.002  0.029  ±  0.008  0.028  ±  0.003  0.027  ±  0.005 

G5  0.023  ±  0.007  0.020  ±  0.005  0.023  ±  0.008  0.022  ±  0.006  0.023  ±  0.004 

Goo  0.029  ±  0.009  0.022  ±  0.005  0.024  ±  0.003  0.022  ±  0.003  0.025  ±  0.004 

_ Absolute  _ D _  _ E _ 

Unit _ Coef.  _ Sham _ Severe  Injury  _ Sham _ Severe  Inj u ry 

mN/mm  Ft  145.7  ±  18.4  163.6  ±  54.1  170.8  ±  24.0  170.7  ±  21.8 

116.2  ±  33.5  106.0  ±  20.6  120.2  ±  32.3  142.8  ± 

8.96  ±  1.15  10.07  ±  3.35  10.49  ±  1.42  10.51  ± 

7.17  ±  2.15  6.52  ±  1.28  7.42  ±  2.07  8.79  ± 

0.746  ±  0.015  0.750  ±  0.009  0.737  ±  0.024  0.756  ± 

0.122  ±  0.006  0.125  ±  0.003  0.122  ±  0.007  0.118  ± 

0.055  ±  0.004  0.054  ±  0.003  0.059  ±  0.004  0.056  ± 

0.028  ±  0.004  0.026  ±  0.002  0.029  ±  0.005  0.024  ± 

0.024  ±  0.007  0.020  ±  0.004  0.025  ±  0.010  0.023  ± 

0.025  ±  0.005  0.025  ±  0.006  0.028  ±  0.006  0.024  ± 

_ A _  _ C _ 

Sham  Severe  Injury  Sham  Severe  Injury 


Ft 

1.017 

± 

0.144 

0.843 

± 

0.145 

0.999 

± 

0.228 

1.187 

± 

0.357 

f3 

1.149 

± 

0.182 

*  0.848 

± 

0.144 

0.992 

± 

0.158 

t  1.361 

± 

0.275 

Po 

1.017 

± 

0.144 

0.843 

± 

0.145 

0.999 

± 

0.228 

1.188 

± 

0.357 

p2 

1.149 

± 

0.182 

**  0.848 

± 

0.144 

0.992 

± 

0.158 

tt  1.362 

± 

0.276 

All  symbols  indicate  a  sta0s0cally  significant  result.  Asterisks  (*,**)  and  daggers  (+,++)  are  comparisons  between  injury  treatment  at  a  particular 
region,  (p  =  resulting  p-value  from  a  student's  t-test  and  13  =  0.025  is  the  Bonferroni  corrected  significance  level).  (*p  =  0.0084,  **p  =  0.0085,  +p  = 
0.0156,  tip  =  0.0156.  All  uncertain0es  are  ±95%CI. 


0.750  ±  0.010 

0.126  ±  0.004 

0.054  ±  0.003 

0.027  ±  0.002 

0.021  ±  0.002 

0.023  ±  0.004 
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terminology:  acoustic  emissions  — >  raw  signal 
Vibro-acoustic  emissions  — >  binned  data  (unnormalized) 


binning  is  important  because  it  identifies  noise  from  actual  signal 

Abstract 

Vibro-Acoustography  (Vibro)  is  an  ultrasound  modality  that  has  proven  to  be  effective  in 
imaging  and  characterizing  tissue  structures  based  on  their  elastic  properties.  While  Vibro  has 
been  implemented  within  a  variety  of  tissue  systems  for  diagnostic  and  therapeutic  effects, 
applications  to  brain  tissue  remain  nascent.  Here  we  investigate  the  ability  of  Vibro  to  detect 
traumatic  brain  injury  (TBI)  in-vivo  via  the  application  of  focused  Vibro-acoustic  emissions  with 
difference  frequencies  between  200-270  kHz.  Findings  from  Vibro  applied  to  rats  (sham  N=6; 
acute  TBI  N=6)  indicate  that  normalized  acoustic  emissions,  captured  near  the  site  of  TBI,  are 
lower  in  magnitude  for  TBI  animals  at  specific  difference  frequencies.  This  ability  of  Vibro  to 
sense  the  acute  changes  in  brain  tissue  caused  by  injury  opens  the  possibility  for  new  TBI 
detection  methods  based  upon  a  Vibro-acoustic  ultrasound  platform. 

Keywords: 

Animal  Models,  Traumatic  Brain  Injury;  Vibro-acoustography,  Rat,  Applied  Radiation  Force 


Introduction 

Traumatic  brain  injury  (TBI)  is  a  leading  cause  of  death  and  disability,  accounting  for 
30.5%  of  all  injury  related  deaths  within  the  U.S.  alone  (Faul  et  al.  2010).  However,  while  TBI 
remains  a  serious  concern  in  domestic  settings,  the  dangers  of  TBI  are  severely  pronounced 
within  combat  zones,  with  over  250,000  cases  of  TBI  reported  by  the  military  between  2002  - 
2012  (Fischer  2013).  The  true  number  of  soldiers  sustaining  TBI  is  believed  to  be  much  higher. 
Military  field  tests  report  missing  40%  of  all  TBI  cases  due  to  high  levels  of  diagnostic 
subjectivity  and  patient  variability  (Miller  et  al.  2010).  While  current  diagnostic  tools  such  as 
computed  tomography  (CT)  and  magnetic  resonance  imaging  (MRI)  are  capable  of  identifying 
TBI,  their  logistics  make  them  ineffective  diagnostic  tools  in  the  field,  where  rapid  triage  can 
significantly  improve  clinical  outcomes.  Here  we  investigate  the  possible  use  of  Vibro- 
Acoustography  (Vibro)  as  a  basis  for  a  portable  TBI-detection  device. 

Vibro  utilizes  confocal  transducers  to  converge  two  ultrasound  beams  of  different  carrier 
frequency  to  a  single  focus,  spatially,  and  generate  vibrations  within  a  tissue  as  a  result  of  the 
applied  radiation  force.  These  tissue  vibrations  will  have  a  frequency  equal  to  the  difference 
between  the  two  carrier  frequencies.  This  emitted  difference  frequency,  acoustic  emission,  can 
be  captured  and  used  to  quantify  the  mechanical  response  of  the  tissue  to  the  applied  radiation 
force  (Fatemi  and  Greenleaf  1998;  Urban  et  al.  2011).  Because  trauma  can  alter  the  mechanical 
properties  of  tissue  (Greenwald  2003;  Boulet  et  al,  201 1  and  2013),  Vibro  induced  vibrations  in 
damaged  tissue  will  generate  different  acoustic  emissions  compared  to  normal  tissue.  In 
essence,  Vibro  exploits  the  differences  in  the  mechanical  properties  between  tissues  to  identify 
and  differentiate  between  normal  and  abnormal  tissue  systems.  Vibro  has  previously  been 
implemented  to  distinguish  between  normal  and  calcified  femoral  arteries  in-vivo  (Pislaru  et  al. 
2008),  detecting  calcifications  within  the  prostate  (Urban  et  al.  2011),  and  the  detection  of  breast 


lesions  in  humans  (Alizad  et  al.  2012).  Vibro  has  also  been  used  to  image  lesions  in  the  liver 
and  to  track  brachytherapy  seeds  used  in  the  treatment  of  prostate  cancer  (Urban  et  al.  2011). 

Here  we  test  the  hypothesis  that  acoustic  emissions  generated  by  transcranially 
delivered  Vibro  are  sensitive  enough  to  detect  the  presence  of  TBI  in  rat  brains  with  focal 
trauma  generated  by  a  controlled  cortical  impact  (CCI). 

MATERIALS  AND  METHODS 
Transducer  description 

A  dual-element  2  MHz  focused  ultrasound  transducer  (Figure  1)(Sonic  Concepts,  Model 
H-148,  Woodinville  WA)  was  used  to  generate  the  acoustic  radiation  force  needed  to  generate 
tissue  vibrations.  This  transducer  has  a  radius  of  curvature  of  63.2  mm,  placing  the  center  of  the 
focus  51 .74  mm  from  the  exit  plane.  The  inner  element  has  an  inner  diameter  of  22.6  mm  and 
outer  diameter  of  48.29  mm.  The  outer  element  has  an  inner  diameter  of  49.05  mm  and  outer 
diameter  of  64.0  mm.  Linear  simulations  of  both  elements  operating  in  unison  (Figure  2)  predict 
a  focus  that  measures  8  mm  axially  and  1  mm  radially  as  measured  at  the  pressure  half¬ 
maximum  value. 

Each  transducer  element  was  driven  independently  but  in  a  synchronized  fashion  by  a 
pair  of  function  generators  (33220A  20MHz  Function/Arbitrary  Waveform  Generator,  Agilent, 
Santa  Clara,  CA)  through  separate  amplifiers  (Model  A150  RF  Power  Amplifiers,  E&l, 

Rochester,  NY).  Driving  signals  were  configured  through  the  use  of  associated  pressure 
outputs,  as  measured  in  degassed  water  with  a  needle  hydrophone  (HNR-100  S/N  1370, 

ONDA,  Sunnyvale,  CA),  to  ensure  that  each  transducer  element  contributed  an  equal  pressure 
to  the  mutual  focus  for  each  of  the  difference  frequencies  considered:  200-270  kHz,  inclusive,  in 
steps  of  10  kHz.  These  selected  difference  frequencies  lay  within  10-13.5%  of  the  transducers 
2  MHz  center  frequency  and  are  within  the  efficiency  band  of  the  transducer.  The  net  pressure 
generated  by  both  transducer  elements  at  a  given  difference  frequency  was  used  to  determine 
that  the  generated  spatial-peak-temporal-average  (SPTA)  intensity  ranged  from  0.1  -  0.7  W/cm2. 


During  experiments  dFU  within  this  intensity  range  was  applied,  at  a  given  difference  frequency, 
for  1000  [is.  Acoustic  emissions  from  the  tissue  were  recorded  for  1200  ps. 

TBI  induction. 

All  animal  procedures  were  approved  by  the  Institutional  Animal  Care  and  Use 
Committee  (IACUC)  of  the  University  of  Washington,  Seattle,  WA.  Male  Sprague  Dawley®  rats 
(Charles  River,  Wilmington,  MA),  weighing  between  351  ±  80  g,  were  split  into  two  experimental 
groups:  a  sham  control  group  (N=6)  and  an  acute  TBI  group  (N=6).  Three  additional  sham  rats 
were  used  to  determine  the  signal-to-noise  ratios  for  various  intensities  of  applied  ultrasound. 
(Three  more  rats  were  used  to  determine  the  difference  between  signals  captured  in  rats  with 
skull  and  those  without.  )| 

Animals  were  anesthetized  with  isoflurane  (5%  for  induction  and  1 .5  -  3%  for 
maintenance,  with  an  oxygen  flow  rate  of  1 .5  L/min)  then  placed  on  a  stereotactic  surgical 
platform  with  metal  bars  that  were  fixed  in  the  rats’  ear  canals  to  stabilize  their  heads  for 
surgery,  injury  induction,  and  data  collection.  The  head  was  then  shaved  and  the  remaining  hair 
removed  with  Nair®  (Church  &  Dwight  Co,  INC.,  Princeton,  NJ,  USA).  Following  epilation, 
subcutaneous  lidocaine  (0.08  ml_,  20  mg/mL)  and  bupivicaine  (0.08  ml_,  0.32  mg/mL)  were 
administered  and  spread  across  the  surgical  site  by  manual  palpation.  A  heating  pad  was 
placed  under  the  body  to  maintain  a  core  body  temperature  near  37°C  throughout  the 
experiment. 

An  incision,  along  the  sagittal  plane,  was  then  made  from  behind  the  ears  to  between 
the  eyes.  Following  this,  an  area  of  skull  approximately  1  cm  in  width  was  fully  exposed  by 
removal  of  the  periosteum.  A  small  bone  drill  was  used  to  create  a  circular  hole  through  the 
skull,  at  a  location  3  mm  anterior  and  2  mm  to  the  left  of  the  lambda,  just  exposing  the 
surface  of  the  brain.  A  rongeur  was  used  to  widen  this  cranial  window  radially  to  a  diameter  of 
approximately  3  mm,  in  order  to  accommodate  the  tip  of  a  controlled  cortical  impact  (CCI) 


Martin  Suarez  4/15/14  9:51  PM 


Comment  [1]: 

add  this  in  ...doesn’t  just  using  the  word 
additional  accomplish  this  meaning?  -LK 


device  (AMS  201,  AmScien  Instruments,  Richmond,  VA,  USA).  This  marked  the  end  of  surgery 
for  the  sham  group. 

For  each  animal  in  the  acute  TBI  group  the  tip  of  the  CCI  device  was  positioned  directly 
over  the  cranial  window  and  lowered  to  the  surface  of  the  brain.  The  device  was  then  fired  to 
induce  a  2  mm  deep  impact  into  the  brain  (tip  speed  of  1.15  m/s),  consistent  with  the  protocol  of 
Lighthall  (1 988)  as  refined  by  Dixon  et  al.  (1 991 ).  This  marked  the  end  of  surgery  for  the  TBI 
group. 

Application  of  diagnostic  focused  ultrasound  (dFU) 

All  Vibro  studies  were  conducted  under  anesthetic  conditions.  The  face  of  the 
transducer  was  covered  with  a  plastic  coupling  cone  filled  with  degassed  water  for  optimal 
ultrasound  transmission  to  the  brain  (Figure  1).  The  front  of  this  coupling  cone  was  sealed  with 
a  sheet  of  acoustically  translucent  latex,  the  rear  was  closed  with  a  focused  polyvinylidene 
fluoride  (PVDF)  hydrophone  (Harisonic®  Y-107-7,  Olympus  Corporation,  Tokyo,  Japan).  The 
focus  of  the  hydrophone  was  selected  such  that  it  is  confocal  with  the  transducer.  The 
hydrophone  was  then  connected  to  an  amplifier  ((Hydrophone  Booster  Amplifier,  Precision 
Acoustics  Ltd.,  Dorchester,  Dorset,  UK)  and  oscilloscope  (LeCroy,  Chestnut  Ridge,  NY,  USA) 
to  allow  visualization  of  acoustic  emissions  from  the  brain.  These  acquired  waveforms  were 
then  linked  to  a  desk-top  computer  (Dell,  Round  Rock,  TX,  USA)  for  MATLAB-based 
(MathWorks,  Natick,  MA,  USA)  data  processing. 

The  transducer  and  hydrophone  were  then  mounted  on  a  three-way  micro-positioning 
system  (Velmex,  INC.,  East  Bloomfield,  NY,  USA)  to  allow  the  ultrasound  focus  to  be  moved  to 
any  anatomical  point  of  interest.  A  laser  guidance  system  was  implemented  in  order  to  identify 
the  ultrasound  focus  during  experiments.  These  lasers  were  affixed  to  the  coupling  cone  such 
that  their  intersection  marked  the  center  of  the  transducers’  focus.  Utilizing  the  lasers,  the 
transducer’s  position  was  adjusted  via  the  micro-positioner  such  that  the  ultrasound  focus 
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occurred  5  or  8  mm  below  the  surface  of  the  skull  (as  desired).  The  focus  could  then  be  moved 
to  each  of  the  two  target  locations  in  our  study,  relative  to  the  site  of  injury,  defined  by  their 
coronal  position  on  the  skull  and  the  depth  of  dFU  (Figure  3).  Acoustic  emissions  from  the  brain 
across  the  range  of  predetermined  difference  frequencies  (200  -  270  kHz,  in  steps  of  10  kHz) 
were  collected  twice  for  each  difference  frequency.  Collected  waveforms  were  exported  to 
MATLAB®  for  analysis.  Within  MATLAB  the  two  waveforms  from  each  difference  frequency 
were  then  averaged  to  yield  a  final  mean  acoustic  emissions  signal.  These  final  mean  acoustic 
signals  were  collected  three  times  for  all  difference  frequencies  and  from  each  target  location. 

Immediately  following  data  acquisition  rats  were  euthanized  with  390  mg/ml 
Beuthanasia®-D  via  intraperitoneal  injection. 

Data  analysis  of  acoustic  emissions 

Raw  acoustic  emissions  from  rat  brain  were  captured  via  the  PVDF  hydrophone  across 
a  frequency  range  of  10  kHz  -  15  MHz.  During  data  analysis  only  emissions  within  the  range  of 
10  kHz  -  2.5  MHz  were  considered  (Figure  4).  Following  acquisition,  acoustic  emissions  were 
decomposed  via  Fast  Fourier  Transform  (FFT)  to  examine  their  frequency  components. 
Frequency  bins  centered  at  each  difference  frequency  with  a  width  of  7  kHz  where  then 
created  to  consolidate  portions  of  the  acoustic  emission  to  be  considered  Vibro-acoustic 
signal.  Additional  frequency  bins  centered  half-way  between  each  applied  difference 
frequency  with  a  width  of  3  kHz  were  also  created  to  consolidate  background  noise 
within  the  range  of  the  applied  difference  frequencies.  This  narrow  frequency  range  for 
noise  was  chosen  so  that  it  would  not  be  skewed  by  the  large  spectral  magnitudes  found 
at  the  carrier  frequencies,  (frequency-frequency  plot,  Figure  5  — >  reassign  location  of 
this  picture. ) 

Preliminary  findings  demonstrated  that  Vibro-acoustic  signals  from  tissue  varied 
significantly  with  the  spectral  magnitude  at  the  carrier  frequencies.  Acoustic  emissions  also 


varied  with  background  noise  in  the  range  of  our  applied  difference  frequencies.  To  account  for 
these  affects  two  normalizations  were  developed: 


N.= -  N,  =  , - — - 

Cf  +  C2  y/Cf  +  C2*RMS(jj) 

A  is  the  maximum  spectral  magnitude  from  within  the  frequency  bin  of  an  acoustic 
emission  generated  by  a  given  applied  difference  frequency.  Ci  and  C2  are  the  spectral 
magnitudes  of  the  two  carrier  frequencies  at  this  same  difference  frequency.  r| is  a  numerical 
set  containing  the  root-mean-square  (RMS)  values  determined  for  each  of  the  7  noise 
bins.  Both  N1  and  N2  yield  non-dimensional  measures  of  Vibro-acoustic  signal  magnitude, 
whereas  unnormalized  signal  was  represented  by  the  spectral  voltage  of  the  acquired  signal  for 
each  applied  difference  frequency. 

N1  strictly  seeks  to  compensate  for  fluctuations  in  the  spectral  magnitudes  of  the  carrier 
frequencies.  Because  induced  tissue  vibrations  are  not  expected  in  the  range  of  the  carrier 
frequencies  changes  in  spectral  magnitude  at  the  carriers,  captured  by  the  hydrophone,  may  be 
attributed  to  an  unexpected  increase  in  ultrasonic  intensity.  Dividing  acquired  emissions  by  the 
carriers  diminishes  the  influence  of  variation  in  the  carriers  on  the  signal  generated  by  vibrating 
tissue.  N2  seeks  to  address  fluctuations  at  the  carriers  as  well  as  noise  within  the  range  of  the 
applied  difference  frequencies.  This  noise  may  arise  from  variations  in  rat  skull  anatomy  as  well 
as  from  bubbles  present  within  the  ultrasound  gel  or  the  coupling  cone.  In  these  cases  any 
captured  emissions  would  be  a  sum  of  scatter  from  skull  or  bubbles  as  well  as  actual  Vibro- 
acoustic  signal.  Dividing  acquired  emissions  by  noise  and  the  carriers  together  allows  a  better 
representation  of  the  signal  generated  by  vibrating  tissue  alone. 


Statistical  Analysis 


The  magnitudes  of  Vibro-acoustic  signals  acquired  from  sham  and  TBI  rats  were 
compared  via  the  Mann-Whitney  U  test  (Sigma  Plot,  Systat  Software  Inc.,  San  Jose,  CA, 
USA),  for  each  location  and  difference  frequency.  In  all  comparisons  a  p-value  <  0.05 
was  taken  to  denote  statistical  significance  between  the  signals  of  TBI  and  sham  rats. 


Data  analysis  of  signal-to-noise  ratio  (SNR) 

Acoustic  emissions  from  three  sham  rats  were  collected  across  a  range  of  applied 
intensities:  0.4,  0.5,  0.6,  and  0.72  W/cm2  SPTA.  This  was  done  to  determine  if  the  application  of 
transcranial  Vibro  could  generate  measurable  Vibro-acoustic  signals  using  intensity  levels  below 
the  FDA  limit  of  0.72  W/cm2  SPTA.  For  calculation  of  the  SNR  the  signal  component  was 

taken  to  be  the  RMS  value  of  the  spectral  magnitudes  within  the  frequency  bin  of  the 
applied  difference  frequency.  The  noise  component  was  taken  to  be  the  RMS  value  of  r| 
(as  previously  defined). 


RESULTS 

Signal-To-Noise  Ratio  (SNR) 

At  an  applied  intensity  of  0.4  W/cm2  SPTA  SNRs  at  the  FR8  location  were  determined  to 
fall  between  3.2  ±  0.9.  SNRs  from  the  FL8  location  fell  between  6.25  ±  2.0. 

Vibro-acoustic  signals  for  sham  versus  TBI  rats 

Analysis  of  un-normalized  Vibro-acoustic  signals  from  sham  versus  TBI  rats  indicated 
that  statistically  significant  differences  were  present  at  several  locations  and  applied  difference 
frequencies  (Table  1),  especially  at  the  FL5  and  FL8  locations  (see  Figure  BLAH  for  an 
example  of  these  results  at  FL8).  In  all  cases  unnormalized  Vibro-acoustic  signals  for  TBI  rats 
were  higher  than  those  for  sham  rats,  across  all  applied  difference  frequencies.  Additionally, 
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in  all  un-normalized  signals  the  spectral  magnitudes  measured  at  the  carrier  frequencies 
were  higher  in  TBI  versus  sham  rats.  However,  once  normalized,  Vibro-acoustic  signals 
from  sham  rats  appeared  higher  than  those  for  TBI  rats  across  all  applied  difference 
frequencies  (see  Figures  BLAH  for  examples  at  FL8). 


Discussion 

In  order  to  verify  that  our  captured  acoustic  emissions  represented  more  than 
simple  backscatter  from  the  skull  we  applied  Vibro  to  rats  (N=3)  with  and  without  skull  at 
the  same  locations  and  depths  as  in  the  primary  study,  each  rat  acted  as  a  self  control. 
Signals  acquired  with  skull  were  compared  to  those  without  skull  via  the  Mann-Whitney 
U  test  (Sigma  Plot,  Systat  Software  Inc.,  San  Jose,  CA,  USA),  once  again  taking  p-values 
<  0.05  to  be  statistically  significant.  Comparisons  indicated  that  the  acoustic  emissions 
acquired  from  rats  with  and  without  skull  for  any  location  and  depth  were  not  statistically 
significant.  Additionally,  Vibro  applied  with  an  intensity  of  0.4  W/cm2  SPTA,  which  marks 
the  lower  end  of  the  applied  intensity  range,  was  able  to  generate  acoustic  emissions  in  which 
signal  components  were  discernible  from  noise  (SNR  >  1).  From  these  analyses  it  follows 
that  measurable  Vibro-acoustic  signals  with  SPTA  intensities  within  FDA  guidelines  can 
be  produced  from  the  application  of  trans-cranial  Vibro.  Furthermore,  comparison  of 
these  Vibro-acoustic  signals  from  sham  versus  TBI  rats  indicates  that  Vibro  is  sensitive 
to  the  changes  in  tissue  elasticity  that  occur  following  trauma.  Given  these  findings 
Vibro  presents  a  lucrative  platform  for  new  TBI  detection  technologies. 


Future  Direction  — >  remove  citation  by  Mao  and  Wagner  2011 

With  the  basis  set  by  these  experiments  work  will  continue  in  order  to  further 
refine  the  application  of  Vibro  for  the  detection  of  TBI.  In  future  studies  additional 
measures  will  be  considered  in  order  to  eliminate  fluctuations  in  transducer  output  and 
amplifier  gain.  While  the  normalizations  (N1  and  N2)  implemented  here  should  minimize 
the  influence  of  these  fluctuations,  further  reduction  may  improve  the  ability  to  detect 
Vibro-acoustic  signals  from  vibrating  tissue.  Additionally,  future  work  would  consider 
larger  animals  that  have  sustained  closed-head  TBI,  such  as  from  a  focused  blast,  this 
would  eliminate  any  influence  the  cranial  window  has  on  emitted  acoustic  emissions. 
Larger  animals  would  also  ensure  that  the  ultrasound  focus  lays  entirely  within  tissue 
and  is  not  effected  by  the  skull  which  is  stiffer  and  less  likely  to  respond  to  Vibro. 

While  Vibro  presents  a  novel  method  by  which  to  interrogate  tissue  for  the 
presence  of  trauma  its  focal  nature  entails  that  scanning  large  regions  would  be 
inefficient.  It  appears  that  Vibro  would  be  most  successful  if  paired  with  a  secondary  TBI 
imaging  modality  with  less  resolution  and  specificity  but  a  broader  field  of  view.  This 
would  combine  the  high  precision  of  Vibro  with  the  capability  to  examine  large  regions 
quickly. 

Conclusion 

Vibro-acoustography  (Vibro)  applied  trans-cranially  is  sensitive  to  the  changes  in 
tissue  elasticity  that  follow  trauma.  As  such,  it  represents  a  foundation  for  novel  TBI 
detection  devices  that  can  be  deployed  in  the  field. 
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Figure  Captions 

Figure  1A:  H-148  HIFU  transducer  (black),  confocal  PVDF  hydrophone  (orange),  and  plastic 
coupling  cone  with  aiming  lasers.  To  couple  the  transducer  to  the  rat’s  skull,  the  cone  is  filled 
with  degassed  water  and  sealed  with  latex.  The  volume  of  water  in  the  cone  can  be  adjusted 
with  a  pair  of  syringes  connected  to  the  cone’s  sides  (not  seen  here). 

Figure  1 B:  Mathematical  simulation  of  the  maximum  pressure  field  generated  by  the  H-1 48 
transducer.  The  transducer  face  is  depicted  on  bottom.  In  the  pressure  field,  red  indicates  high 
pressure  and  blue  low  pressure. 

Figure  2:  Surgical  procedure  schematic.  A  hole  was  drilled  in  the  skull  3  mm  in  diameter  (red 
circle),  located  3  mm  anterior  and  2  mm  to  the  left  of  the  lambda,  which  is  the  visible  rear  suture 
joint  on  the  surface  of  the  skull.  Focused  ultrasound  was  applied  at  the  locations  denoted  by 
the  yellow  stars,  which  are  both  5  mm  anterior  of  the  lambda  and  3  mm  left  and  right, 
respectively.  The  focus  was  aimed  5  and  8  mm  deep  at  each  location. 

Figure  3A:  Sample  signal  from  an  in  vivo  experiment,  Af=220  kHz.  The  top  panel  displays  the 
waveform  collected  by  the  PVDF  hydrophone,  unfiltered.  The  bottom  panel  displays  the  Fourier 
transform  of  the  waveform.  The  two  peaks  occur  at  the  carrier  frequencies  for  the  two 
elements,  in  this  case  2.1 1  and  1 .89  MHz.  The  emissions  at  the  difference  frequency  are 
orders  of  magnitude  below  those  at  the  carriers,  but  in  most  cases  they  stand  out  from  those  at 
neighboring  frequencies  (noise). 


Figure  3B:  Representative  sample  of  the  amplitude  of  the  measured  acoustic  emissions  of  a 
given  data  collection  run,  with  applied  difference  frequency  along  the  horizontal  axis  and 
measured  spectral  frequency  along  the  vertical  axis.  Each  column  represents  the  same  signal 
such  as  that  seen  Figure  3A  but  after  application  of  a  Fast  Fourier  transform.  The  visible 
diagonal  trend  in  the  plot  reaffirms  that  at  each  difference  frequency,  the  acoustic  emissions  at 
that  frequency  stand  out  above  the  noise  floor,  indicating  the  presence  of  vibro-acoustography. 

Figure  4A:  Unnormalized  Vibro-acoustic  emissions  from  the  Front  Left  8mm  location.  Asterisks 
indicate  a  statistically  significant  difference. 

Figure  4B:  Vibro-acoustic  emissions  from  the  FL8  location,  normalized  with  the  N1  formula. 

Figure  4C:  Vibro-acoustic  emissions  from  the  FL8  location,  normalized  with  the  N2  formula. 


TABLE 


Table  1 :  Percent  of  difference  frequencies  at  which  Mann-Whitney  U  testing  indicated  a 
statistically  significant  difference  in  Vibro-acoustic  emissions  between  sham  and  TBI  rats. 
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ABSTRACT 

Annually,  an  estimation  of  1.7  million  people  in  the  U.S  suffers  from  traumatic  brain  injury  (TBI),  in  which 
50,000  people  die  and  another  85,000  people  survive  with  long  -term  disabilities.  CT  and  MRI  are  the 
current  diagnostic  tools  for  TBI,  but  the  prevalent  implementation  of  the  devices  is  limited  in 
undeveloped  areas  and  near  the  battlefield.  Consequently,  a  size  deployable  diagnostic  tool  that  is  able 
to  present  quantitative  evaluation  of  TBI  is  needed.  The  purpose  of  the  project  is  to  optimize  an 
ultrasound  elastography  machine  for  detection  and  quantification  of  TBI  based  on  the  theory  that 
abnormal  brain  elasticity  will  be  observed  in  presence  of  hemorrhage  and  edema.  The  goal  is  achieved  in 
three  steps:  develop  and  test  shear  wave  tracking  algorithm  in  breast  phantom  and  alginate,  test, 
reproduce  and  modify  the  algorithm  on  control  mouse  brains,  and  finally  construct  elastograph  on 
severe  TBI  mouse  and  find  the  most  sensitive  parameter  for  TBI  diagnosis  through  quantitative  statistics. 
The  accuracy  of  elastograph  will  be  ensured  by  SSI  machine  in  breast  phantom  and  alginate.  Successfully 
optimized  ultrasound  elastography  machine  displays  quantitative  elastograph  of  brains,  which  simplifies 
diagnostic  procedures.  Furthermore,  implementation  of  the  device  in  hospitals  and  battlefields  enables 
early  TBI  detection,  which  will  potentially  reduce  TBI-related  deaths,  lower  TBI-related  treatment  costs 
and  greatly  enhance  residents'  life  quality. 


Introduction: 

Statement  of  Problem: 

Traumatic  brain  injury  (TBI)  can  be  classified  into  three  levels:  mild,  moderate  and 
severe.  Top  three  causes  of  TBI  are  car  accidents,  firearms  and  falls.12  Different  levels  of 
TBI  will  result  in  various  extents  of  brain  hemorrhage  and  edema.  TBI  symptoms  include 
continuing  and  worsening  headaches,  loss  of  consciousness  or  memory,  and  in  severe 
cases,  dilation  of  pupils.1- 3  Patients  with  mild  TBI  are  able  to  recover  in  a  few  weeks  or 
months.  However,  many  troops  back  from  the  war  suffer  from  moderate  to  severe  brain 
injury,  which  will  never  fully  recover  and  may  need  lifelong  treatment.2 

Operable  scanning  equipment  such  as  Computed  Tomography  (CT)  and  Magnetic 
Resonance  Imaging  (MRI)  has  been  used  for  TBI  diagnosis.3  However,  their  bulky  size  limits 
their  wide  application;  experienced  operators  are  also  required  in  order  to  give  meaningful 
analysis  of  the  image.4  Consequently,  patients  with  brain  injuries  will  have  to  be 
transported  far  away  in  order  to  receive  proper  diagnosis,  which  is  a  waste  of  time  and 
money  and  can  lead  to  further  damage  to  the  brain.4 

Ultrasound  transient  elastography  is  a  prominent  candidate  for  TBI  diagnosis  due  to 
its  deployable  size  and  potential  to  give  a  quantitative  analysis.  According  to  previous 
works  conducted  on  ultrasound  transient  elastography,  it  is  able  to  detect  changes  in  brain 
elasticity  after  schematic  stroke.4- 5  It  is  also  capable  of  measuring  brain  anisotropy  in 
different  regions.6  However,  ultrasound  elastography  is  not  yet  able  to  detect  TBI  due  to 
factors  such  as  low  signal  to  noise  ratio  and  huge  skull  artifacts. 

The  Verasonics  ultrasound  engine  (VUE] 11  is  a  completely  programmable  PC  -based 
ultrasound  device.  By  coupling  it  with  a  conventional  diagnostic  probe,  which  pushes  tissue 
non -destructively  via  acoustic  radiation  forces  from  focused  diagnostic  ultrasound  (FDU), 
VUE  can  be  programmed  to  perform  ultrasound  imaging.  By  post -processing  the  data  and 
modifying  ultrasound  parameters,  VUE  has  the  potential  to  display  brain  elastography. 
According  to  the  previously  proved  hypothesis,  brain  elasticity  will  vary  in  presence  of 
edema  and  hemorrhage.  Consequently,  the  ultimate  goal  of  the  project  is  to  design  and 
optimize  VUE,  so  that  the  machine  will  be  able  to  generate  TBI  sensitive  ultrasound 
elastography  that  is  able  to  detect,  classify  and  quantify  TBI  with  minimum  artifacts. 
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Significance: 

According  to  the  statistics  compiled  by  the  Centers  for  Disease  Control  and 
Prevention  (CDC),  an  estimated  1.7  million  people  suffer  from  TBI  each  year  and  the 
condition  has  become  a  leading  factor  30.5%  of  injury  -related  deaths  in  the  United  States.1' 

3  Moreover,  “TBI  and  concussion”  is  reported  to  be  the  most  common  combat -related  injury 
and  account  for  nearly  25%  of  combat  casualties.2  Although  protection  strategies  such  as 
body  armor  and  Kevlar  helmets  are  available,  the  head  remains  insufficiently  protected.2 
These  advances  increase  the  number  of  survivors  from  the  war  but  also  results  in  more 
soldiers  with  permanent  defects  in  their  brain  due  to  TBI.2  Since  traumatic  brain  injury  is 
influencing  more  and  more  people’s  lives,  especially  the  welfare  of  soldiers  coming  back 
from  war,  proper  diagnostic  methods  are  definitely  needed  in  order  to  perform  the  right 
treatment.  Because  current  approaches  to  TBI,  i.e.  MRI,  magnetic  resonance  elastography 
(MRE)  and  CT,  have  many  limitations  in  terms  of  accessibility,  accuracy,  and  diagnostic  time, 
the  need  for  a  highly  accurate  and  accessible  TBI  diagnostic  approach  becomes  more  and 
more  urgent.  Not  only  medically,  the  project  is  also  technically  significant  because 
ultrasound  elastography  is  a  recently  developed  technology  that  has  only  been  working  on 
breasts  in  the  clinic.  Applications  of  the  technology  on  many  other  organs,  including  liver, 
kidney  and  brain  are  under  extensive  research.  Consequently  if  the  design  and  optimization 
of  VUE  for  brain  imaging  is  successful,  it  will  become  a  new  scientific  breakthrough,  which 
proves  the  sensitivity  of  shear  wave  in  tissues  other  than  breasts.  Additionally,  the 
technology  will  be  widely  applied  in  the  medical  field  because  of  its  capability  to  produce 
real-time  brain  elastographies  with  little  or  no  artifact  from  the  skull.  By  referring  to  the 
shear  modulus  maps  constructed  by  VUE,  the  operator  may  be  able  to  diagnose  brain  injury 
according  to  the  highlights  in  the  image.  Classification  with  respect  to  level  of  injury  can  be 
achieved  by  looking  at  the  amount  of  change  in  elasticity.  Due  to  its  portability  and 
reasonable  price,  VUE  will  ultimately  be  implemented  in  battlefields,  hospitals  and  other 
places  where  effective  imaging  machines  are  not  available,  and  the  broad  application  of  the 
device  will  enhance  the  lives  of  millions  of  people. 
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Ethical  and  social  issues 

Animal  models  in  the  experiment  can  cause  ethical  concerns  because  rats  and  mice 
will  undergo  mild  to  moderate  TBI  surgery,  imaging  and  finally  perfusion.  However,  the 
Institutional  Animal  Care  and  Use  Committee  (IACUC)  approved  all  surgery  protocols, 
anesthesia  and  chemicals  used  in  the  experiment.7  IACUC  is  an  entity  that  regulates  proper 
use  of  laboratory  animals  for  research  and  instructional  purposes.  By  strictly  complying 
with  IACUC  laws  on  animal  handling  and  surgery,  ethical  issues  associated  with  animal 
experimenting  can  be  avoided. 

Successful  research  outcomes  provide  an  accessible  and  portable  approach  for  early 
detection  and  diagnosis  of  traumatic  brain  injury.  This  will  largely  benefit  society  by 
reducing  TBI -related  deaths  and  enhancing  life  quality  of  residents.  In  addition,  annual  cost 
of  TBI  treatment  in  the  U.S  has  been  estimated  to  be  48.3  billion  dollars,  and  individuals 
suffering  from  severe  TBI  tend  to  spend  4  million  dollars  over  the  entire  life.8  Consequently 
early  detection  of  the  disease  prevents  deterioration  of  TBI  and  can  potentially  benefit  the 
society  by  reducing  medical  costs. 

Economic  issues: 

The  project  requires  large  economic  support  because  experiments  consume  huge 
number  of  animal  models,  including  rats,  mice  and  potentially,  pigs.  Drugs  such  as 
Lidocaine,  Buprenorphine  and  Bupivacaine  and  chemicals  such  as  paraformaldehyde  (PFA) 
and  isoflurane  are  also  inevitable  as  they  are  required  for  surgical  procedures  according  to 
IACUC.  Instrumentation  used  in  the  lab  is  SSI  and  VUE,  which  costs  millions  of  US  dollars 
and  needs  regular  updates  and  renewals. 

The  military  is  the  sponsor  for  the  project  because  a  large  portion  of  casualty  from 
war  suffers  from  blast  -induced  neurotrauma  (BINT),  which  is  a  typical  type  of  TBI.4  Owing 
to  the  fact  that  there  is  currently  no  portable  equipment  near  the  battle  field  able  to 
perform  TBI  diagnosis,  the  military  is  interested  in  funding  the  project  since  the 
construction  of  a  portable  and  accurate  TBI  diagnostic  machine  will  greatly  enhance  the 
clinical  condition  and  welfare  of  soldiers.  Other  organizations  such  as  National  Institute  of 
Health  (NIH)  may  also  be  willing  to  fund  the  project  because  the  diagnostic  machine  can  be 
implemented  in  hospitals  and  ambulances,  which  will  greatly  benefit  residents.4 
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Technical  Background 

Theory: 

The  project  aims  to  quantify  TBI  using  ultrasound  elastography  based  on  the 
hypothesis  that  tissue  elasticity  value  will  vary  in  the  presence  of  injury  due  to  hemorrhage 
and  edema.  The  hypothesis  has  been  proven  through  the  preliminary  work  in  the  lab  on  a 
total  of  60  rat  models,  and  conclusions  have  been  drawn  that  normal  brain  tissues  have  a 
significantly  higher  shear  modulus  than  damaged  tissues. 
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Figure  1:  1)  FDU  probe  sends  out  "push"  beams  at  focal  spots,  inducing  SW  propagation  in  imaging  plane.  2) 
Ultrafast  imaging  system  image  shear  wave  propagation  plane  after  "push"9 


Ultrasound  elastography  displays  tissue  shear  modulus  via  focused  diagnostic 
ultrasound  (FDU]  beams.  The  FDU  transducer  (fig.l)  connected  to  VUE  is  able  to  send  out 
sequences  of  radiation  “pushes”  to  induce  non-destructive  tissue  displacements.  The  tissue 
displacements  generated  by  the  sequences  at  one  location  will  then  form  a  downward  wave 
front  (fig.2),  which  is  able  to  propagate  transversely  through  the  tissue.  The  wave  caused  by 
the  induced  tissue  displacements  is  called  a  shear  wave  (SW).  After  sending  out  "pushes", 
the  probe  will  start  taking  images  (fig.l  and  fig.2)  at  an  adjustable  frequency  with  the 
highest  value  of  8000Hz,  which  records  videos  of  (SW)  propagation  at  different  time 
frames.  Post  processing  and  filtering  of  the  pictures  will  produce  a  video  clip  (fig.3), 
displaying  the  propagation  of  wave  front  caused  by  the  downward  tissue  displacement. 


1  ms 


6ms 


13ms 


Figure  2:  SW  propagation  three  time  points  after  "push"  in  B-mode.  White  scale  indicated  tissue 
displacement.10 


Figure  3:  Fragment  of  SW  propagation  over  time  in  alginate:  red  indicates  local  maximum  vertical 
displacements.  Notice  that  the  SW  starts  on  the  right  and  propagate  toward  the  left. 


Due  to  the  fact  that  similar  shapes  of  wave  fronts  will  form  along  the  length  of  tissue 
when  SW  moves  through,  the  propagation  can  thus  be  tracked  by  cross  correlating 
displacement  data  between  two  different  time  frames  at  each  depth  of  the  image  (fig.4)  .9' 10 
The  cross  correlation  algorithm  calculates  pixels  moved  in  order  to  most  accurately  overlap 
two  waveforms,  and  the  resulting  number  of  pixels  is  equal  to  the  SW  propagation  distance. 
As  a  consequence,  SW  velocity  can  be  found  with  distance  moved  in  the  known  amount  of 
time  difference.  Referring  to  papers  in  the  related  field,  shear  modulus  (G)  is  related  to  SW 
velocity  in  the  form  G  =  3*p*c!,in  which  p  is  the  density  of  tissue  under  imaging  and  c  is 
the  speed  of  SW.10  By  applying  this  equation,  the  shear  modulus  for  each  grid  can  be 
calculated  and  put  together  to  construct  a  map.  In  the  preliminary  stage  of  the  design 
period,  a  commercialized  Aixplorer  Multiwave  Ultrasound  System  (SSI)  for  breast  imaging 
will  be  used  as  a  standard  test  for  the  accuracy  of  VUE. 


Press  E  te»  escape 

*  10  ^  (Reference  Frame)  Depth  -  27.  Frame  -  86, 


Figure  4:  Cross  correlation  method.  The  first  window  is  the  waveform  at  an  early  time  frame,  and  the  second 
window  is  the  waveform  at  the  latter  time  frame.  The  third  window  shows  the  correlation  function  of  the  two 
waves. 
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Review  of  Literature: 

CT  and  MRI,  shown  in  fig.5  have  been  two  commonly  used  ways  in  diagnosis  of  TBI 
over  the  past  twenty  years4.  The  way  they  work  is  to  provide  visualization  of  brain 
structures,  shown  in  fig.  6. 12 

MRI  has  been  widely  used  clinically  to  diagnose  central  nervous  system  disease, 
cardiovascular  disease,  tumor,  and  head  and  neck  health.13  MRI  is  very  safe  because  it  uses 
non— ionizing  radiation.  The  imaging  modality  involves  the  use  of  paramagnetic  contrast 
agents,  thinner  slices  and  surface  coils,  which  provides  high  spatial  resolution  and 
specificity.13  Compared  to  CT,  MRI  is  able  to  give  more  detailed  information  about  the 
brain:  small  brain  damage,  usually  1— 2mm,  can  be  seen  on  an  MRI.15  One  major  limitation  of 
MRI  is  that  it  cannot  be  operated  on  patients  with  ventilation  equipment  or  pacemaker 
because  of  the  magnetic  field  it  creates.  Additionally,  according  to  Sahraian  and  Eshaghi14, 
MR  imaging  has  no  pathological  specificity  when  used  to  infer  microstructural  information 
(i.e.  any  change  in  brain  tissue  composition  (inflammation,  edema,  or  demyelination)  can 
change  MRI  signal). 

Diffusion  tensor  MR  imaging  (DT— MRI)  specifically  targets  brain  white  matter.  The 
technology  studies  the  diffusion  of  water  molecule  in  tissues  containing  a  large  amount  of 
fibers,  such  as  white  matter.35  Water  tends  to  diffuse  fastest  along  the  direction  of  fibers. 
Consequently,  based  on  how  water  diffuses  in  tissue,  DT-MRI  is  able  to  produce  3D 
representation  of  the  tissue.35  Studies  have  shown  that  this  imaging  modality  is  able  to 
have  a  more  complete  assessment  of  brain  health  and  disease.34  35  However,  many  patients 
find  MRI  intolerable  as  the  process  is  long  and  often  takes  place  in  a  confined  space.12 

CT,  on  the  other  hand,  is  better  at  detecting  the  presence  of  fresh  blood  around  the 
brain  and  is  mostly  used  in  acute  TBI  diagnosis.12  Both  MRI  and  CT  require  highly  trained 
and  experienced  operators  so  as  to  give  a  meaningful  interpretation  of  the  image,  which 
increases  the  chance  of  human  bias  and  error  during  diagnosis.12  Moreover  CT  has  poor 
resolution,  low  post— acute  value  and  can  induce  cell  death,  cancer  or  other  health  problems 
due  to  the  exposure  to  x-rays,  which  is  a  high-level  ionizing  radiation.12 
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Figure  5:  [Left]  Schematic  graph  of  MRI  (left)  and  CT  scanner  (right)  12  [Right]  Brain  image  from  MRI 
(left)  and  CT  scanner  (right)  4 


Magnetic  Resonance  Elastography  (MRE)  is  a  new  imaging  technique  (fig6)  that  has 
recently  been  studied  in  vivo  on  both  small  animals  and  humans.  M RE  is  a  non— invasive 
technique;  it  works  by  providing  visualization  as  well  as  quantitative  viscoelastic 
properties  of  the  tissue  being  imaged,  shown  in  fig.6.9  MRE  determines  the  shear  modulus 
of  tissue  by  obtaining  the  propagation  speed  of  shear  wave  generated  via  an  external 
vibrator  based  on  the  fact  that  shear  wave  travels  at  various  speeds  in  tissues  with 
different  elasticity.9  By  measuring  the  elasticity  change  in  tissue,  which  often  times  is  an 
indication  of  disease,  MRE  will  be  able  to  detect  the  region  of  abnormality. 

The  technology  has  been  applied  to  several  parts  of  human  bodies  such  as  breast, 
skeletal  muscle  and  prostate  to  assess  tissue  properties.6 19  Meanwhile,  extensive  research 
targets  the  use  of  MRE  to  non— invasively  study  the  viscoelastic  properties  of  human  brains 
in  vivo.16’17’ 18  Recent  reports  from  Murphy  et  al.  suggests  a  decrease  in  the  elastic 
properties  of  human  brains  with  neuroinflammatory  and  neurodegenerative  disorders, 
including  Alzheimer’s  disease,  using  MRE.20U  Hamhaber  et  al.  used  MRE  to  measure  SW 
velocity  in  human  brain,  and  found  the  average  velocity  to  be  1.88  m/s.21 

However,  when  MRE  is  used  to  evaluate  local  brain  properties,  large  discrepancies 
appear  in  literature  in  terms  of  the  exact  shear  modulus  value  in  the  specific  regions  of  the 
brain.  Referring  to  the  previous  research  conducted  on  MRE  in  vivo,  the  elasticity  of  brain  is 
significantly  different  between  white  matter  and  grey  matter.  According  to  McCracken16, 
the  brain  white  matter  has  a  shear  modulus  of  10.7±1.4  kPa  and  the  grey  matter  has  a 
shear  modulus  of  5.3±1.3  kPa;  while  Michael  A.  Green19  and  his  colleagues  who  measured 
brain  viscoelasticity  via  3D  reconstruction  method  got  a  shear  modulus  value  of  2.7±0.1 
kPa  for  brain  white  matter  and  3.1±0.1  kPa  for  brain  grey  matter.  Many  similar 
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experiments  were  performed  and  the  value  of  shear  modulus  varied  from  1-12  kPa.16'  18< 19' 
21  The  lack  of  consistency  can  be  attributed  to  variation  in  sample  types,  methods,  and 
imperfection  in  the  technology  such  as  artifacts. 

MRE  needs  a  very  long  acquisition  time,  usually  20  minutes,  in  order  to  obtain  a 
quantitative  evaluation  of  the  tissue  under  imaging.22  Short  scan  times,  which  is  usually  the 
case  in  human  study,  leads  to  acquisition  of  data  with  poor  spatial  resolution.22  Due  to  its 
limited  spatial  resolution,  MRE  is  only  able  to  give  averaged  mechanical  properties  over  the 
whole  brain  or  over  large  ROIs,  which  disables  its  application  to  diagnose  localized 
injuries.22  Moreover,  the  amplitudes  of  high  frequency  vibrations  generated  by  MRE 
attenuates  rapidly  in  the  middle  of  the  brain,  which  significantly  lowers  the  signal  to  noise 
ratio  of  the  deformation  data.22  Therefore,  the  application  of  this  technology  as  a  routine 
diagnostic  tool  for  both  research  and  clinical  use  is  limited. 


Figure  6:  [Left]  MRE  set  up  with  external  vibrator,  transducer  mount,  driving  coils  and  bite  bar.  [Right]  MRE  in 

19 

vivo  image  (a)  T2-weighted  MR  image;  (b)  G'  (kPa)  image;  (c)  G"(kPa)  image,  (d).  overlapped  image  of  b)  and  a). 


Ultrasound  transient  elastography  was  first  developed  in  the  1990's. 25  The 
technology  uses  the  same  concept  as  MRE,  in  that  tissue  viscoelasticity  can  be  evaluated  via 
shear  wave  propagation,  an  example  of  shear  wave  is  shown  in  fig. 7.  Different  from  MRE, 
ultrasound  transient  elastography  utilizes  a  transient  mechanical  excitation  of  tissue 
instead  of  prolonged  vibration,  which  naturally  separates  shear  wave  from  compression 
waves  as  shear  wave  travels  three  times  slower.9' 23  This  technology  effectively  reduces 
wave  overlapping  and  thus  increases  SNR.  In  order  to  more  specifically  study  the  shear 
wave  propagation  in  medium,  the  Institute  Langevin  developed  an  ultrafast  ultrasound 
imaging  acquisition  system  in  1997,  which  allows  storing  of  2D  data  acquired  at  a  frame  rate 
of  5000  to  30,000  images  per  second,  followed  by  conversion  of  which  into  a  map  of 
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shear  modulus.9’ 23  The  combination  of  the  two  ideas  (transient  excitation  and  ultrafast 
imaging)  results  in  today’s  shear— imaging  technique.9 


Figure  7:  Comparison  of  ultrasound  wave  to  shear  wave.  [Left]  Longitudinal  ultrasound  wave  induce  medium 
movement  in  the  same  direction  of  propagation.  [Right]  Shear  wave  results  from  perpendicular  medium 
movements  relative  to  the  direction  of  propagation.9 


Shear  imaging  technique  has  been  widely  applied  to  study  many  organs,  including 
thyroid,  gastro -intestinal  tract,  kidney,  and,  most  successfully,  breast  and  liver.24’  25> 26’ 27 > 32 
BK.  Kim  et  al.  reported  that  shear  imaging  provided  an  alternative  non— invasive  method  to 
access  severity  of  liver  fibrosis.28  Additionally,  the  technology  is  able  to  provide  more 
detailed  prognostic  information  than  traditional  liver  biopsy  by  tracking  dynamic  changes 
in  fibrotic  burdens.28  According  to  an  in  vivo  study  conducted  at  Institute  Curie  on  human 
breasts,  50  patients  with  or  without  breast  lesions  were  involved,  and  shear  imaging  was 
demonstrated  to  clearly  differentiate  between  benign  and  malignant  tissue  (fig.8).29 
Further  study  by  W.A.  Berg  et  al.  in  939  patients  with  breast  cancer  showed  that  shear 
imaging  had  high  specificity  for  breast  lesion  characterization.30 


Figure  8:  Ductal  infiltrating  carcinoma.  [Left]  B-mode  image  shows  only  blurry  edge  of  the  lesion.  [Right] 
Shear  imaging  highlights  the  lesion  with  observably  high  shear  wave  velocity,  which  corresponds  to  high 
tissue  stiffness.29 
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Because  of  shear  imaging’s  high  spatial  resolution  to  provide  local  tissue  stiffness, 
efforts  have  been  made  to  apply  the  technique  on  brains  in  small  animals  such  as  rats  and 
mice.  However,  the  technology  fails  to  measure  brain  viscoelasticity  due  to  huge  amount  of 
artifacts  resulting  from  the  brain’s  tiny  size  plus  the  presence  of  skull,  which  reflects  shear 
waves  from  the  edge  and  causes  them  to  overlap.4-6  31  Consequently,  extensive  research  is 
going  on  regarding  optimization  of  ultrasound  transient  elastography  that  can  overcome 
current  roadblocks,  thus  allowing  accurate  evaluation  of  brain  stiffness. 

Previous  Work: 

In  Mourad’s  lab,  experiments  and  data  acquisition  on  rat  and  mouse  brains  after 
surgeries  were  carried  out  using  SSI  with  linear  SL  15--4  MHz  transducer  in  order  to  prove 
the  hypothesis  on  decreased  tissue  elasticity  due  to  injury.  Because  the  SSI  machine  is 
designed  for  imaging  in  breast  tissue,  its  application  on  brain  results  in  huge  artifacts  due 
to  SW  overlaps  at  the  brain  edge  and  skull  reflections.  An  example  of  SSI  imaging  result  on 
mouse  is  shown  in  fig.9. 


Figure  9:  Shear  modulus  map  of  mouse  brain  with  two  cranial  windows  imaged  by  SSI.  The  arrow  indicates 
site  of  artifacts  caused  by  the  skull,  which  has  much  higher  elasticity  value  than  actual  brain.4 

In  stroke  surgeries,  ultrasound  transient  elastography  was  performed  on  mouse 

brains  24  hours  and  72  hours  after  occlusion  in  the  middle  cerebral  artery.5  Shear  moduli 
were  obtained  for  both  ipsilateral  and  contralateral  hemispheres  and  the  data  were 
compared  to  that  of  the  control  mouse  brains.  According  to  the  result  of  the  experiment, 
shown  in  fig.  10,  stroke  mouse  brains  after  24  hours  displayed  significant  shear  modulus 
decrease  in  the  ipsilateral  hemisphere  and  significant  shear  modulus  increase  in  the 
contralateral  hemisphere,  both  with  p  values  less  than  0.01.  No  difference  was  found 
between  72— hour  brains  (p>0.05)  but  the  shear  modulus  for  24— hour  and  72— hour  mouse 
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brains  varied  greatly  in  both  hemispheres  (p<0.01).  The  experiment  proves  the  ability  of 
this  specific  version  of  ultrasound  elastography  to  detect  subtle  changes  in  brain  elasticity 
in  the  presence  of  artifacts  after  induction  of  ischemic  stroke.5 


Figure  10:  [Left]  Normalized  Ipsilateral  Shear  Modulus  for  control,  24-hour  and  72-hour  mice  (P<0.01) 
[Middle]  Comparison  of  normalized  mean  shear  modulus  of  ipsilateral  and  contralateral  hemisphere  for  all 
three  populations:  Control,  24-hour  and  72-hour  mice.  *P<0.01  **P=0.04  [Right]  Normalized  Contralateral 
Shear  Modulus  for  control,  24-hour  and  72-hour  mice  (P<0.01)  5 

In  TBI  surgeries,  brain  shear  modulus  values  were  evaluated  acute  and  24  -hour  after 
moderate  TBI  simulations  using  the  same  SSI  machine.  A  total  of  60  live  rats  were  used  in 
the  experiment,  and  each  20  of  them  were  24— hour,  acute  and  control  rats  respectively.  TBI 
surgery  was  performed  through  a  1— mm  cortical  impact  on  the  right  hemisphere,  and  skulls 
were  removed  on  both  sides  of  the  brain,  leaving  the  middle  line  intact.  MATLAB  was  used 
to  contour  the  images  for  clear  visualization,  and  AN OVA  was  applied  to  provide  statistical 
analysis,  including  mean,  standard  deviation  and  p  value  of  the  two  hemispheres.  The 
brains  were  perfused  after  imaging  to  obtain  histology  results  of  the  injury.  According  to 
coronal  images,  the  average  rat  brain  shear  modulus  was  different  between  24  ■hour  TBI 
and  control  rats  while  no  significant  difference  was  observed  between  acute  TBI  and 
control  rats.  The  experiment  showed  that  24  -hour  after  TBI,  the  injured  hemisphere  had 
significantly  decreased  stiffness,  which,  according  to  histology  results,  was  due  to  presence 
of  edema  and  hemorrhage.  Fig.  11  illustrates  the  experiment’s  statistical  results. 

Normalized  (Visible  TBI) 


Acute  Acute  24  Hour  24  Hour 

Contra  Ipsi  Contra  Ipsi 
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Figure  11:  Bar  plots  of  normalized  acute,  24-hour  contralateral  and  ipsilateral  brain  stiffness  values,  in 
which  the  24-hour  brain  has  decreasing  stiffness  in  injured  hemisphere  and  increasing  stiffness  in 
contralateral  hemisphere.  (P<0.01) 

The  preliminary  data  prove  that  brain  elasticity  will  change  due  to  edema  and 
hemorrhage  caused  by  TBI  or  stroke,  and  the  change  is  statistically  detectable  even  using 
the  SSI  machine,  which  is  not  optimized  and  has  huge  artifacts.  This  result  supports  the 
hypothesis  for  the  project.  Therefore,  modification  and  optimization  of  VUE  to  enable 
detection  and  classification  of  brain  injuries  becomes  a  promising  goal. 

Outstanding  Technical  Issues: 

Presence  of  artifacts  in  SSI  was  the  most  outstanding  issue  in  the  beginning  of  the 
project.  Because  of  the  unavoidable  artifacts  from  skull,  it  became  impossible  to  determine 
the  real  shear  modulus  of  brain  under  imaging.  Consequently,  the  newly  generated  SW 
tracking  algorithm  could  not  track  the  SW  movie  exported  from  SSI.  As  a  result,  the  design 
of  VUE  was  prohibited  due  to  lack  of  a  standard  comparison. 

Another  technical  problem  that  emerged  at  the  outset  was  the  low  signal  to  noise  ratio 
of  displacement  data  exported  from  VUE.  As  SW  propagated  through  tissue,  the  intensity  of 
which  dissipated,  and  there  were  reflections  coming  from  boundaries  due  to  the  tiny  size  of 
the  brain.  This  made  the  detection  of  wave  very  hard. 

The  brains  were  then  put  into  alginates  to  eliminate  reflection  from  skull.  However, 
since  the  SW  impedance  of  alginate  did  not  match  that  of  the  brain,  although  SW  appeared 
to  be  perfect  in  alginate,  it  became  a  mess  when  hitting  the  brain— alginate  interface,  and 
disappeared  or  broke  apart  when  entering  the  brain.  Because  there  is  no  reliable  literature 
that  indicates  brain  SW  impedance,  many  alginate  molds  with  different  percentages  varied 
by  0.1%  were  tested  in  order  for  SW  to  propagate  smoothly.  This  procedure  significantly 
slowed  down  the  experiment  progress. 

Design  of  Tools,  Device,  and  Experiments: 

Overview: 

The  goal  of  this  project  is  to  build  and  optimize  a  deployable  sized  ultrasound 
elastography  device  VUE  capable  of  highlighting  and  classifying  areas  of  damage  caused  by 
moderate  and  mild  traumatic  brain  injury  (TBI).  The  device  should  be  able  to  differentiate 
injury  types  (i.e.  hemorrhage,  edema  and  traumatic  axonal  injury)  and  classify  damage 
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extents  (mild,  moderate]  from  TBI.  Traditional  approaches  such  as  MRI  and  CT  are  mainly 
limited  by  their  size  and  cost;  newly  developed  MRE  uses  external  vibration  source  to 
induce  tissue  displacement  inside  the  brain,  but  the  technology  requires  very  long 
acquisition  time,  which  limits  its  application  for  routine  diagnosis.  VUE  is  an  alternative 
promising  candidate  to  achieve  the  goal:  the  diagnostic  probe  in  VUE  sends  out  radiation 
“pushes"  to  induce  non-destructive  tissue  displacement,  after  the  pushing  sequence  VUE  is 
able  to  take  images  as  SW  propagates  through  the  tissue  and  convert  the  data  to  shear 
modulus  values,  enabling  real  time  visualization  of  elasticity  images.  Because  VUE  is  fully 
programmable,  modifications  can  be  made  on  pushing  sequence,  duration  and  averaging 
pixels,  which  are  capable  of  avoiding  SW  overlap  due  to  skull  reflection.  Successful 
programming  and  design  will  enable  VUE  to  produce  brain  elasticity  maps  that  quantifies 
and  classifies  various  types  of  TBI  with  high  level  of  accuracy. 

Phase  1:  The  goal  of  this  phase  is  to  construct  stiffness  maps  in  breast  phantom  and 
alginate.  This  is  achieved  by  post-processing  of  ultrasound  data  from  VUE  and  performing 
shear  wave  tracking  in  MATLAB.  The  outcome  of  this  phase  is  the  construction  of  shear 
modulus  images  that  are  able  to  display  accurate  breast  phantom  stiffness  and  detect 
lumps. 

Phase  2:  The  goal  of  this  phase  is  to  construct  stiffness  mapping  in  rats  with  moderate  TBI. 
This  is  achieved  by  post-processing  of  ultrasound  data  from  VUE,  and  performing  SW 
tracking  in  MATLAB.  Additional  work  such  as  rat  surgery,  perfusion,  histology  and  artifacts 
prevention  are  required  in  the  process.  The  outcome  of  this  phase  is  the  construction  of 
artifact-free  shear  modulus  image  that  is  able  to  highlight  moderate  TBI  sites. 

Phase  3:  The  goal  of  this  phase  is  to  apply  the  outcomes  from  phase  2  to  detect  mild  TBI  in 
rat  brains.  This  is  achieved  by  modification  of  previous  MATLAB  codes  in  phase  2  to  fit  for 
mild  injury  detection.  Mild  TBI  surgery  on  rats  is  required  in  the  process.  The  outcome  of 
this  phase  is  the  construction  of  a  shear  modulus  image  that  is  able  to  highlight  as  well  as 
differentiate  mild  TBI  sites  from  moderate  ones. 

Revised  Methods: 

In  phase  1  of  the  project,  SW  tracking  in  breast  phantom  and  alginate  turned  out  to 
be  successful  at  early  stage.  However,  in  phase  2,  due  to  the  small  size  of  rodent’s  brain,  and 
the  unavoidable  skull  interference,  SW  propagation  in  brain  tissue  in  vivo  became  very 
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noisy  and  untraceable.  Consequently,  further  experiments  were  carried  out  in  vitro.  Control 
mouse  brains  were  taken  out  and  submerged  in  different  percentage  of  alginate  molds  to 
optimize  imaging.  SW  tracking  algorithm  was  also  modified  to  optimize  brain  elastograph. 
In  phase  3,  both  24-hour  TBI  and  control  mice  were  studied  to  determine  the  most  TBI- 
sensitive  imaging  parameters,  and  give  a  quantitative  idea  of  the  stiffness  of  mice  brains. 
The  purpose  of  this  significant  revision  of  the  design  process  is  to  create  an  imagine  media 
that  is  big  enough  to  fit  the  transducer  head,  and  which  also  has  the  closest  SW  propagation 
impedance  as  the  brain  itself,  so  that  when  FDU  generates  pushes,  SW  will  propagate 
smoothly  from  the  alginate  media  into  the  brain  tissue,  and  the  tracking  algorithm  will  be 
able  to  capture  any  abnormity  of  the  propagation  caused  by  the  presence  of  injury. 

Materials  and  Methods: 

Animal  models  of  mice  were  used  in  both  phase  2  and  3  of  the  study.  Drugs  used 
associated  with  animal  surgeries  were  B-euthanasia,  Lidocaine,  Buprenorphine,  and 
Bupivacaine.  Other  materials  included  breast  phantom,  alginate  powder,  sodium  chloride, 
and  sodium  dithionite.  The  devices  used  were  SSI,  VUE,  Pneumatic  Impact  Device33,  and  L4- 
7  ultrasound  transducer. 

The  design  strategy  was  composed  of  two  major  parts.  One  part  was  to  optimize  SW 
tracking  algorithm.  To  achieve  this  goal,  two  sets  of  algorithms  with  different  correlation 
methods  were  developed  and  tested  in  alginate;  various  time  frame  lag  for  SW  tracking 
were  chosen  and  compared;  weighted  space  averaging  and  interpolation  of  correlation 
coefficients  were  also  coded  into  the  tracking  algorithm.  The  other  part  was  to  optimize  SW 
propagation  through  media.  To  achieve  this  goal,  various  percentages  of  alginates  were 
made  and  imaged;  different  "push”  duration  and  active  "push”  elements  were  chosen  and 
tested;  averaging  effect  of  time  frames  and  pixels  were  studied;  color  percentage  of  SW 
visualization  was  also  manipulated  to  show  the  best  propagation. 

Data  with  different  combination  of  "push"  duration  and  active  "push”  elements  were 
collected  using  SSI  and  VUE  at  the  same  position  on  breast  phantoms,  alginates  and  brain- 
containing  alginate  molds  (the  set-up  of  the  alginate  molds  is  shown  in  fig.12]. 
Displacement  data  from  both  devices  were  exported,  as  well  as  the  elastograph  from  SSI. 
The  exported  displacement  data  were  post-processed  according  to  the  strategies  above, 
and  then  fed  through  different  SW  tracking  algorithms,  giving  elasticity  maps.  Mean  and 
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standard  deviation  were  calculated  from  each  generated  elastograph,  and  the  results  of 
breast  phantom  and  alginate  were  compared  to  the  standard  image  from  SSI. 

Mouse  brains  were  used  to  test  the  alginate  recipe  and  SW  tracking  algorithm  to  give 
the  best  and  most  consistent  result.  Control  animals  were  sacrificed  by  cervical  dislocation 
(mice)  or  B-euthanasia  (rats).  Their  brains  were  later  put  into  alginate  molds.  TBI  animals 
were  given  severe  TBI  surgery,  in  which  a  small  hole  in  the  animal  skull  will  be  drilled  on 
the  right  hemisphere  to  just  fit  the  cortical  impactor,  then  the  animal  will  be  impacted  by 
1mm  into  the  brain.  The  animal  will  be  sutured  up  10  minutes  after  impact  and  sacrificed 
after  24  hours. 


Costs:  (1  page  excel): 

Table  1:  experimental  costs 


Item 

Cost 

VeraSonics  Ultrasound  Engine 

$58,000 

L7-4  Transducer 

$2400 

180-2 2 Og  Sprague  Dawley  Rats 

$24.25/rat 

Mito  mice 

(Free)  bred  in  lab 

Nitrogen  tank 

$13. 66/tank 

Oxygen  tank 

$11. 51/tank 

Roungeurs 

$335.00  each 

Forceps  #45 

$25.50  each 

Buprenorphine  Hydrochloride  Injection 

$31/lml 

Lidocaine 

$5/bottle 

Bupivacaine 

$5/bottle 

B-euthanasia 

$98/bottle 

Isoflurane 

$10/bottle 

Alginate 

$14.10/lb 

Simethicone 

$3.24/1  floz 

Sodium  Chloride 

$36.10/500  g 

Isopropanol 

$209.85/4L 

Sodium  Dithionite 

$15.03/500  g 
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Design  Process: 

In  part  1  of  the  design  process,  data  were  obtained  from  breast  phantom  and  2.5% 
alginates  because  of  their  homogeneous  medium  property.  The  resulting  images  were 
compared  to  that  of  SSI  in  the  same  ROI  in  terms  of  mean,  standard  deviation,  minimum 
and  maximum  values. 


Figure  12:  Example  of  correlation  overview:  correlation  occurs  between  the  1st  frame  and  the  9th  frame,  the  2nd 
frame  and  the  10th  frame  and  so  forth.  K  is  equal  to  8  in  this  case. 

Correlation  methodology: 

Recall  in  fig.3  that  SW  propagation  data  is  composed  of  111  2D  maps  over  a  total 
time  period  of  20ms.  Each  2D  map  represents  the  vertical  displacements  of  elements  in  the 
ROI  at  one  time  frame,  and  correlation  method  was  used  to  find  SW  velocity  and  thus  tissue 
elasticity.  The  overall  correlation  pattern,  shown  in  fig.  12  as  an  example,  was  designed  to 
correlate  data  in  the  nth  frame  with  that  in  the  (n+k]th  frame,  the  (n+l)th  frame  with  that  in 
the  ((n+l]+k)th  frame,  and  so  forth,  in  which  k  is  the  number  of  frame  lags  used  in  the 
tracking.  The  resulting  SW  velocity  and  shear  modulus  was  assigned  to  all  the  horizontal 
pixels  involved  in  the  k  frames. 

Two  separate  correlation  methods  were  designed  and  tested.  Linear  correlation 
method  correlates  the  entire  line  of  data.  Max  value  correlation  finds  the  start  of  wave 
front,  and  defines  region  of  correlation  to  be  pixels  starting  from  wave  front  to  the  end  of 
the  line,  and  same  region  k  frames  later,  as  shown  in  fig.13.  For  each  set  of  data,  post¬ 
processing  was  performed  using  two  correlation  methods  with  frame  lag  k  chosen  to  be  5, 

10  and  15  respectively. 
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Figure  13:  [Left]  Same  line  of  displacement  at  two  time  frames  separated  by  8  units,  correlation  occurs  over  the 
entire  line  of  data.  [Right]  Same  line  of  displacement  at  two  time  frames  separated  by  8  units,  correlation  occurs  in 
regions  of  red  boxes. 

Interpolation  of  cross  correlation: 

Recall  from  fig.4  that  a  cross  correlation  function  gives  correlation  coefficient  and 
horizontal  movement  of  two  waveforms  in  pixels  in  order  for  them  to  overlap.  This  leads  to 
a  shortage  of  the  algorithm  in  that  SW  movement  less  than  1  pixel  cannot  be  detected. 
Consequently,  interpolation  of  the  cross  correlation  function  was  performed  in  MATLAB 
using  "polyfit”,  which  refits  the  function  with  a  polynomial,  and  increases  the  accuracy  and 
spatial  resolution  of  elastograph  by  detecting  SW  movements  within  1  pixel. 

Weighted  average  based  on  correlation  coefficient: 

Weighted  average  of  SW  velocity  and  thus  shear  modulus  was  performed  based  on 
correlation  coefficients.  Correlation  coefficients  vary  from  0  to  1  for  no  correlation  to 
perfectly  positive  correlation.  These  coefficients  are  used  to  weight  the  shear  modulus 
when  averaging.  As  SW  movement  overlaps  spatially  in  each  tracking  (i.e.  1st  and  9th 
frames,  2nd  and  10th  frames  and  etc.),  averaging  algorithm  in  space  was  used  to  enhance  the 
accuracy  of  results.  Fig.14  shows  how  spatial  averaging  was  performed  in  the  horizontal 
direction  at  a  certain  depth. 


Figure  14:  Weighted  average  example:  spreadsheet  shows  the  obtained  SW  velocity  after  each  tracking  at  a  certain 
geometric  depth.  Weighted  average  is  performed  using  correlation  coefficient  and  vertically,  where  there  are 
overlapping  values  in  the  same  x  pixel,  as  circled  out  in  orange. 


18 


Cross  correlation  coefficient  is  used  for  weighted  average  of  SW  velocity  in  each  line. 
According  to  figure  4,  weighted  average  is  performed  through  each  column  among  the  non 
zero  pixels,  and  finally  form  a  line  of  averaged  SW  velocity.  The  velocities  are  then  used  to 
calculate  shear  modulus,  and  the  same  algorithm  is  applied  to  each  line  of  displacement 
data. 

In  part  2  of  the  design  process,  different  percentage  of  alginate  modes,  ranging  from 
4.4%  to  4.7%,  were  made.  0.9%  saline  was  used  as  solvent  to  protect  brain  damage  from 
osmotic  pressure,  and  the  solution  was  kept  at  4°C  to  slow  down  the  alginate  solidifying 
reaction.  2%  of  isopropanol,  0.3ml/100ml  of  Simethicone  and  0.5g/1.0L  of  Sodium 
Dithionite  were  added  to  increase  the  homogeneity  and  decrease  bubble  formation  in 
alginate.  The  goal  was  to  find  an  alginate  recipe  that  has  same  SW  impedance  as  that  of 
brain,  which  is  equal  to  the  product  of  medium  density  and  SW  velocity.  Two  layers  of 
alginate  were  poured  into  the  container,  in  which  the  bottom  layer  acted  as  a  platform  for 
brain  to  sit  on,  allowing  the  top  layer  to  cover  the  brain.  As  alginate  tends  to  change  its 
elastic  property  over  time,  data  were  collected  between  10  "30  minutes  after  solidification 
of  the  top  layer.  Animals  were  sacrificed  after  the  settlement  of  bottom  alginate  layer,  and 
the  brain  was  kept  in  0.9%  saline  solution  to  reduce  its  change  of  property.  SW  movie 
collection  was  performed  at  four  different  positions  around  the  brain,  shown  in  fig.15.  At 
each  position,  five  various  combinations  of  "push”  parameters  were  chosen  respectively 
for  data  acquisition,  and  the  filtering  parameters  were  altered  to  reduce  noise  in  SW 
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Figure  15:  [Left]  Schematics  of  brain  -alginate  setup  and  imaging  positions  14.  In  position  1  and  2,  SW  initiates  in 
alginate,  and  propagates  through  the  alginate  into  the  brain.  In  position  3  and  4,  SW  initiates  in  brain,  and  propagates 
through  the  brain  into  alginate.  The  pink  star  indicates  site  of  TBI  injury  for  TBI  animals.  [Right]  A  real  example  of 
schematics  in  B-mode  image:  the  dotted  square  indicates  ROI,  in  which  the  data  were  post  -processed  in  MATLAB,  the 
pink  arrow  indicates  position  of  FDU  "push",  the  yellow  arrows  indicates  the  brain,  and  the  horizontal  pink  line  at  the 
top  is  the  actual  active  pushing  element  in  the  transducer. 
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Figure  16:  [Left]  Table  of  five  parameter  combinations:  data  were  collected  accordingly  in  order  to  find  the 
push  type  that  works  best.  [Right]  Data  filtering  options:  the  parameters  can  be  manually  altered  to  smooth 
out  SW  propagation  in  brain. 

The  brain  hemispheres  were  accurately  contoured  (accuracy  up  to  0.5mm)  out 
during  elastograph  analysis  according  to  B— mode  images.  The  hemisphere  difference  and 
image  quality  were  first  determined  qualitatively  based  on  mean  value  and  homogeneity  in 
alginate,  presence  of  artifacts,  and  inter— hemisphere  difference  in  brain.  Then,  parameters 
with  equal  image  quality  underwent  quantitative  analysis  in  the  same  specified  region  for 
mean,  standard  deviation,  and  AN OVA  tests.  The  goal  of  the  quantitative  analysis  is  to  find  a 
parameter  combination  that  is  most  sensitive  to  diagnose  TBI. 
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Data: 

Breast  Phantom  and  Alginate  Elastograph: 


At  the  end  of  summer  2012,  a  preliminary  SW  tracking  algorithm  was  designed. 
There  was  no  averaging  or  interpolation  of  data.  The  algorithm  was  tested  in  breast 


phantom  use  SSI. 


Figure  17:  [Left]  Breast  phantom  elastograph  in  the  scale  of  20  kPa  generated  by  SSI.  [Right]  Same  breast 
phantom  elastograph  in  the  scale  of  20kPa  generated  by  SW  tracking  of  the  SW  data  exported  from  SSI. 


In  fall  2012,  the  SW  tracking  pattern  was  modified  to  largely  increase  spatial 
resolution  and  two  SW  tracking  algorithms  were  developed  [linear  and  max  value 
tracking).  The  effects  of  different  frame  lags  k  [recall  that  it  is  the  parameter  that 
determines  the  frame  separation  of  correlation  data  pairs)  were  also  studied.  Data 
acquisition  and  comparison  were  conducted  in  2%  alginate  in  both  SSI  and  VUE. 
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Figure  19:  Elastograph  generated  by  SW  tracking  of  data  exported  from  SSI  in  fig.18  ROI  using  different 
methods.  [Left]  LINEAR  tracking  +  5  frame  lag  [Middle  Left]  MAXVALUE  tracking  +  5  frame  lag  [Middle  Right] 
LINEAR  tracking  +  10  frame  lag  [Right]  MAXVALUE  tracking  +  10  frame  lag 
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Figure  20:  2%  alginate  elastograph  generated  by  SW  tracking  of  data  exported  from  VUE  in  fig.18  ROI  using 
different  methods.  [Left]  LINEAR  tracking  +  5  frame  lag  [Middle  Left]  MAXVALUE  tracking  +  5  frame  lag 
[Middle  Right]  LINEAR  tracking  +  10  frame  lag  [Right]  MAXVALUE  tracking  +  10  frame  lag 


According  to  the  result,  k  =  10  frame  lag  became  preferred  to  5  because  it  gave  more 
homogeneous  and  consistent  results  in  alginate.  However,  in  order  to  increase  spatial  and 
temporal  resolution  of  tracking,  k  =  8  frame  lag,  which  had  comparable  results  to  that  of  k  = 


10,  was  chosen. 
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Figure  21:  [Top]  SSI  elastograph  ofROI  in  2%  alginate.  [Bottom]  Alginate  elastograph  generated  by  SW 
tracking  of  data  exported  from  SSI  in  ROI  using  [Left]  LINEAR  tracking  +  8  frame  lag  [Right]  MAXVALUE 
tracking  +  8  frame  lag 


Figure  22:  2%  alginate  elastograph  generated  by  SW  tracking  of  data  exported  from  VUE  in  same  ROI  using 
[Left]  LINEAR  tracking  +  8  frame  lag  [Right]  using  MAXVALUE  tracking  +  8  frame  lag 

Statistical  analysis  was  conducted  on  a  total  of  12  elastographs  collected  from  3 

different  2%  alginate  using  VUE,  and  compared  to  SSI  results  to  determine  the  accuracy  of 

the  two  tracking  methods.  Max  value  tracking  method  was  chosen  accordingly. 

Table  2:  Statistical  results  of  mean  shear  modulus  difference  between  VUE  and  SSI  and  standard  deviation  of 
VUE  using  max  value  and  linear  tracking  methods. 


2%  Alginate 

Mean  f%  diffl 

STDEV 

Max 

5.02% 

1.0449 

Linear 

8.1% 

0.9270 

Moving  on,  more  optimizations  were  applied  to  SW  tracking  algorithm.  Weighted 
average  and  interpolation  of  cross  correlation  functions  were  encoded  to  increase  shear 
modulus  accuracy  and  give  rise  to  even  higher  spatial  resolution. 


24 


Figure  23:  2%  alginate  elastograph  generated  using  data  from  VUE.  [Left]  Without  averaging  or  interpolation 
[Right]  With  averaging  and  interpolation 

Rat  Brain  Elastograph: 

After  optimizing  alginate  elastograph,  control  and  TBI  rat  brains  submerged  in  2% 
alginate  were  studied  in  winter  2013.  However,  the  result  made  no  sense:  SW  became  a 
mess  when  hitting  the  brain  and  the  elastograph  displayed  unreasonably  high  or  low 
stiffness  in  and  around  brain. 


Figure  24:  [Right]  Control  rat  brain  submerged  in  2%  alginate  in  VUE's  B-mode  ultrasound  image.  The  box 
indicates  ROI.  Ultrasound  "push"  locates  at  the  red  dotted  line.  [Middle]  Static  SW  movie  from  VUE  when  SW 
encounters  the  brain  -alginate  interface  [Left]  Brain  elastograph  generated  using  data  from  VUE. 


Figure  25:  [Right]  TBI  rat  brain  submerged  in  2%  alginate  in  VUE's  B-mode  ultrasound  image.  The  box 
indicates  ROI.  Ultrasound  "push"  locates  at  the  red  dotted  line.  [Middle]  Static  SW  movie  from  VUE  when  SW 
encounters  the  brain  -alginate  interface  [Left]  Brain  elastograph  generated  using  data  from  VUE. 


Mouse  Brain  Elastograph: 

Because  SW  tracking  algorithm  was  working  very  well  in  alginate,  the  hypothesis  for 
such  messiness  in  brain— alginate  interface  is  that  SW  impedance  of  2%  alginate  and  brain 
does  not  match.  According  to  literature,  MRE  was  used  to  measure  SW  velocity  in  human 
brains.  The  value  was  used  in  this  study  to  represent  that  in  mouse  brains.  Measurements 
of  SW  impedance  were  conducted  in  different  percentages  of  alginates  in  order  to  match 
that  of  brain.  Recall  that  SW  impedance  equals  to  density  *  SW  velocity  in  the  medium. 
Mouse  brains  were  used  in  the  study. 

Table  3:  The  measured  SW  impedance  of  3.3%  to  6%  alginate.  SW  did  not  appear  in  5.5%  and  6%  alginate. 
According  to  the  table,  5.3%  alginate  has  the  closest  SW  impedance  to  mice  brain  assuming  the  literature 
measurement  on  human  brain  is  correct. 


Alginate  (96) 

G  (kPa) 

Density  (kg/mA3) 

Velocity  (m/s) 

Density 

Velocity 

Impedance  ! 

1.9113 
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0.828571163 

3.3 
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AVG  = 

1029.66667 

0.83548324 

860.269248 

1.8357 
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3.302S 
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3.7 

3.569 

1004 

1.13224043 
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1096.70217 

2.958 
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1.030776406 

4.5467 
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4.3034 
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1.24328657 
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1.178333496 
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1  4209789 
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1.53275983 
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9.5839 
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8.3632 
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1.733211755 

AVG  - 

1029.66667 

1.68008211 

1729.92455 

5.8666 

1032 

1.451638986 

! 

S5 

NA 

Brain 

llOO 

1.88 

2068 

i 

1  s 

NA 

Referring  to  table  3,  SW  imaging  was  conducted  in  4.4%,  4.5%,  4.6%  and  4.7% 
control  brain— alginate  molds.  Meanwhile,  5  different  combinations  of  ultrasound  push 
parameters  [recall  from  figure  16)  were  chosen  for  SW  movies  collection.  Figures  below 
show  the  elastograph  of  4.4%  and  4.7%  alginates  with  brain. 


Figure  26: 


u  alginate. 
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Figure  27:  Shear  modulus  maps  using  1-5  push  parameters  of  [Top]  4.4%  alginate  alone  and  [Bottom]  Control 
mouse  brain  in  4.4%  alginate. 


i 
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Figure  29:  Shear  modulus  maps  using  1-5  push  parameters  of  [Top]  4.7%  alginate  alone  and  [Bottom] 
Control  mouse  brain  in  4.7%  alginate 


Regarding  to  the  result  above,  4.7%  alginate  was  chosen  to  be  the  most  optimized 
percentage  for  mouse  brains.  Push  parameter  number  5  was  considered  best  for  4.7% 
alginate  SW  imaging.  Therefore,  one  additional  experiment  with  the  chosen  setup  was 
carried  out  in  alginate  to  test  parameter  5,  and  SW  propagation  movies  in  brain  were 
studied  to  reconfirm  the  results. 
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Figure  31:  SW  propagation  movie  of  control  mice  in  4.7%  alginate 


Starting  spring  2013,  a  total  of  5  control  mice  and  5  24— hour  TBI  mice  were  imaged 
in  vitro  using  the  chosen  alginate.  The  imaging  protocol  was  modified  in  that  each  brain 
will  be  imaged  between  20— 30min  after  exposure  and  at  four  positions,  using  various  push 
parameters.  The  positions  were  indicated  below  in  a  chronic  order.  Expectations  were  that 
SW  imaging  would  be  able  to  highlight  TBI  site  with  lower  shear  modulus.  An  example  of 
each  is  shown  below. 
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Figure  32:  [Top  Right]  B-mode  Image  of  Control  mouse  brain  in  positionl.  [Others]  Shear  modulus  map  of  control 
mouse  brain  in  4.7%  alginate  in  position  1  with  5  parameters.  Boxes  contours  the  brain,  circles  indicates  the  two 
hemisphere. 
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Figure  33:  [Top  Right]  B-mode  Image  of  Control  mouse  brain  in  position2.  [Others]  Shear  modulus  map  of  control 
mouse  brain  in  4.7%  alginate  in  position  2  with  5  parameters 
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Right]  B-mode  Image  of  Control  mouse  brain  in  position3.  [Others]  Shear  modulus  map  of  control 
4.7%  alginate  with  parameter  1,  3,  5  in  Position  3 
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Right]  B-mode  Image  of  Control  mouse  brain  in  position4.  [Others]  Shear  modulus  map  of  control 
4.7%  alginate  with  parameter  1,  3,  5  in  Position  4 
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Figure  36:  [Top  Right]  B-mode  Image  of  TBI  mouse  brain  in  positionl.  [Others]  Shear  modulus  map  of  24hr  TBI 
mouse  brain  in  4.7%  alginate  in  position  1  with  5  parameters.  Boxes  contours  the  brain,  circles  indicates  the  two 
hemisphere. 
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Figure  37:  [Top  Right]  B-mode  Image  of  TBI  mouse  brain  in  position2.  [Others]  Shear  modulus  map  of  24hr  TBI 
mouse  brain  in  4.7%  alginate  in  position  2  with  5  parameters 
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Figure  38:  [Top  Right]  B-mode  Image  of  TBI  mouse  brain  in  position3.  [Others]  Shear  modulus  map  of  24hr  TBI 
mouse  brain  in  4.7%  alginate  with  parameter  1,  3,  5  in  position  3,  purple  lines  indicate  edge  of  normal  hemisphere. 
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Figure  39:  [Top  Right]  B-mode  Image  of  TBI  mouse  brain  in  position4.  [Others]  Shear  modulus  map  of  24hr  TBI 
mouse  brain  in  4.7%  alginate  with  parameter  1,  3,  5  in  position  4,  purple  lines  indicate  edge  of  TBI  hemisphere 


According  to  the  qualitative  analysis  in  5  control  and  TBI  mice,  position  1  and  2,  and 
parameter  4  and  5  turned  out  to  have  equal  brain  image  quality.  Consequently,  quantitative 
analyses  of  these  images  were  carried  out  to  determine  the  sensitivity  of  the  parameters  to 
TBI.  The  circled  regions,  indicating  the  two  hemispheres,  were  evaluated  for  mean  shear 
modulus  and  compared  to  each  other.  Position  1  and  position  2  were  analyzed  separately 
at  first  to  study  the  effects  of  time  variance  (due  to  imaging  order:  10—15  minutes 
difference)  and  position  of  hemisphere  (near  or  far  from  the  ‘push’). 
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Figure  40:  Statistical  analysis  (parameter  4  and  5)  of  shear  -modulus  in  the  two  hemispheres  of  control 
mouse  brains  at  position  1  and  position  2.  Box  plots  show  mean,  distribution  and  p  values  of  shear  modulus 
values  in  left  and  right  hemispheres. 
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Figure  42:  Statistical  analysis  (parameter  4)  of  shear  -modulus  in  the  two  hemispheres  of  TBI  mouse  brains 
[Left]  Position  1  and  position  2  TBI  side  [Right]  Position  1  and  position  2  normal  side 
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Figure  43:  Statistical  analysis  (parameter  5)  of  shear  -modulus  in  the  two  hemispheres  of  TBI  mouse  brains 
Left]  position  1  and  [Right)  position  2 
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Figure  44:  Statistical  analysis  (parameter  5)  of  shear  -modulus  in  the  two  hemispheres  of  TBI  mouse  brains 
[Left]  Position  1  and  position  2  TBI  side  [Right]  Position  1  and  position  2  normal  side _ 


After  confirming  that  time  variance  of  imaging  and  the  position  of  hemisphere 
relative  the  ‘push’  have  no  significant  effect  on  the  stiffness  of  brain  tissue,  further  analysis 


was  conducted  combing  the  data  from  position  1  and  position  2. 


Figure  45:  Statistics  results  (parameter  4  and  5)  of  shear  modulus  in  the  two  hemispheres  of  control  mouse 
brains. 
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Figure  46:  Statistics  results  (parameter  4  and  5)  of  shear  modulus  in  the  two  hemispheres  of  TBI  mouse 
brains. 
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Figure  47:  Statistics  results  (parameter  4  and  5)  of  shear  modulus  in  the  two  hemispheres  of  Control  mouse 
and  the  Normal  hemisphere  of  the  TBI  mouse.  _ 


Design  Decision: 

In  the  computation  part  of  the  project,  I  was  responsible  to  design,  test,  and 
optimize  SW  tracking  algorithm.  I  designed  the  tracking  pattern,  two  separate  tracking 
methods  and  the  approach  to  manipulate  SW  data  in  space  and  time  to  give  the  best  shear 
modulus  image.  I  decided  to  do  weighted  average  and  designed  the  cross  correlation 
interpolation. 

In  the  experiment  part  of  the  project,  I  was  responsible  to  perform  surgeries  to 
rodents,  set  up  imaging  protocols,  and  design  an  alginate  setup  for  brain  to  sit  in. 

I  decided  to  work  on  different  percentage  of  alginates  as  an  approach  to  optimize  SW 
imaging.  I  decided  to  test  alginate  percentage  ranging  from  4.4%  to  5%.  I  designed  a  two- 
layer  alginate  setup  that  could  cover  up  the  brain  with  the  correct  orientation. 
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In  the  analysis  part  of  the  project,  I  was  responsible  to  post  -process  SW  data 
exported  form  both  SSI  and  VUE,  study  the  resulting  shear  modulus  map,  and  gave  further 
modification  or  conclusion.  I  decided  to  compare  mean,  standard  deviation,  max  and  min 
values  in  alginate  stiffness  maps  from  VUE  and  SSI.  I  decided  to  contour  the  brain  and  its 
hemispheres  in  stiffness  maps  for  easy  visualization.  Additionally,  I  decided  to  do 
quantitative  analysis  along  with  qualitative  analysis  in  TBI  and  non— TBI  hemispheres  so  as 
to  determine  the  most  sensitive  parameter  for  diagnosis  of  TBI. 

Analysis  and  Conclusions: 

Breast  phantom  and  Alginate  Studies: 

According  to  fig.  17,  the  shear  moduli  in  the  two  images  in  the  same  color  scale  are 
close  in  the  bottom  half.  However,  the  generated  image  has  very  low  spatial  and  temporal 
resolution,  and  inconsistent  values  on  the  top.  This  result  proves  the  correctness  of  the 
preliminary  SW  tracking  algorithm,  and  also  indicates  future  work  in  the  direction  to 
increase  resolution.  In  fig.  18  and  fig.  19,  temporal  resolution  of  the  generated  shear 
modulus  maps  in  alginate  is  significantly  increased  by  the  modified  SW  tracking  pattern 
[recall  in  section  Correlation  Methodology).  Images  with  frame  lag  equals  to  10  are  more 
homogenous  and  have  mean  G  value  much  closer  to  that  of  SSI  [fig.  18  top)  than  images  with 
frame  lag  equals  to  5.  There  are  still  artifacts  existing  in  the  images  using  SSI  data,  indicated 
by  the  bright  red  pattern,  but  the  MaxValue  Tracking  method  results  in  fewer  artifacts  than 
Linear  Tracking  method  [fig.  1 8  bottom).  Images  generated  from  VUE  have  a  lower  average 
G  value  [fig.  19),  but  there  is  no  artifact  present.  This  is  mostly  because  VUE  only  generates 
one  push  into  the  medium,  while  SSI  generates  4  pushes  at  different  positions  into  the 
medium.  Therefore,  artifacts  from  SSI  are  very  likely  due  to  SW  interference  between 
different  pushes.  In  fig.20  and  fig.21,  images  generated  using  frame  lag  8  have  comparable 
homogeneity  as  frame  lag  10,  but  higher  temporal  resolution,  and  therefore  was  chosen  to 
be  the  standard  lag  used  in  further  studies.  Note  in  fig.20  bottom  that  images  generated 
using  SSI  data  have  same  bright  patterns  as  appeared  in  SSI  images  on  the  top.  VUE  images 
have  average  G  values  much  closer  to  that  from  SSI  [fig.21),  and  the  MaxValue  Tracking 
continues  to  have  better  image  quality  and  mean  shear  modulus.  According  to  the  statistical 
analysis  in  table  1,  MaxValue  method  appears  to  work  3% 
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better  comparing  to  Linear  Tracking  method,  and  the  standard  deviation  of  the  two  tracking 
algorithm  do  not  differ  significantly.  Consequently,  MaxValue  tracking  method  is  applied  in 
later  experiments.  In  fig.22,  both  images  have  very  close  mean  shear  modulus  value. 
However  the  right  image  has  less  stripy  patterns  and  higher  resolution  comparing  to  the  left 
image.  This  implies  that  the  cross  correlation  interpolation  and  weighted  spatial  averaging 
are  working  very  well  in  optimizing  the  elastograph. 

Rat  Brain  Elastograph: 

According  to  fig.23,  a  control  rat  brain  is  sitting  vertically  in  2%  alginate  in  B— mode 
image.  From  the  static  SW  propagation  image  (middle],  it  is  very  obvious  that  when  SW 
reaches  the  brain  -  alginate  interface,  it  breaks  up  and  disappears  in  the  brain.  This  gives 
rise  to  the  messy  shear  modulus  map  (right],  which  has  large  non— traceable  area  on  the  left 
side  and  unreasonable  max  G  value  equals  to  78.6kPa  in  brain.  Nonetheless,  the  shear 
modulus  looks  nice  and  homogenous  in  alginate  on  the  right.  Same  phenomenon  is 
observed  in  fig.24  on  a  24— hour  TBI  brain.  The  shear  modulus  map  on  the  right  contours 
the  brain  edge,  but  fails  to  give  reasonable  stiffness  values  inside  the  brain.  Therefore,  one 
possible  conclusion  is  that  when  SW  is  propagating  from  alginate  to  brain,  the  media 
difference  causes  SW  to  break  apart,  dissipate  and  finally  disappear.  This  discrepancy  can 
be  reduced  if  the  SW  impedance  of  alginate  equals  to  that  of  brain.  Another  possible  reason 
can  be  due  to  the  white  matter  in  rat  brains.  By  looking  at  SW  propagation  movie  in  rats, 

SW  tends  to  break  apart  or  disappear  at  the  white  matter  because  the  structure  is  much 
stiffer  than  grey  matter.  Though  this  phenomena  is  not  observed  in  mice  study  (shown 
below],  it  is  very  likely  that  the  larger  white  matter  structures  in  rats  disturb  SW  more 
significantly  so  that  the  filtering  and  averaging  algorithm,  which  works  on  mice,  are  not 
sufficient  to  solve  the  discrepancy  in  rats. 

Mouse  Brain  Elastograph: 

According  to  table  2,  human  brain  SW  impedance,  referring  to  Weaver  et.al18,  is 
closest  to  that  of  5.3%  alginate.  For  alginates  larger  than  5.3%  (i.e.  5.5%  and  6%],  the 
mechanical  property  of  them  is  too  hard  for  the  ultrasound  pushes  to  induce  detectable 
displacements.  Consequently,  5.3%  alginate  becomes  the  top  limit. 

In  fig.27, 4.4%  alginate  images  are  homogenous  throughout  the  five  parameters, 
however  the  mean  shear  moduli  of  which  are  relatively  high  (~7kPa]  compared  to  greater 
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percentage  of  alginates  (i.e.  4.7%  in  fig.29).  This  is  because  alginate  becomes  harder 
through  time,  and  for  this  particular  set  of  imaging,  the  alginate  was  settled  longer  than  the 
other  ones.  According  to  this  finding,  later  alginate  imaging  were  all  controlled  to  occur  in 
the  same  time  frame  (30—40  min].  Reasonable  shear  modulus  values  in  brain  indicate  that 
SW  successfully  propagates  through  the  brain  from  alginate.  Nonetheless,  brain  images  still 
appear  to  have  unexpected  artifacts,  elastograph  from  parameter  1  is  especially  messy. 

This  result  implies  that  the  SW  impedance  of  alginate  and  brain  gets  closer,  but  there  is  still 
space  to  improve.  In  fig.29, 4.7%  alginate  images  are  very  homogeneous,  especially 
parameter  5.  Brain  elastograph  appear  to  have  tiny  spots  of  artifacts  near  the  top  of  the 
brain.  One  possible  reason  for  the  artifacts  is  SW  being  disturbed  when  propagating  through 
white  matter,  which  is  a  lot  stiffer  than  grey  matter.  Other  than  that,  the  elastograph  shows 
reasonable  brain  shear  modulus,  the  brain  seems  invisible  in  alginate  and  there  is  no 
obvious  inter— hemisphere  difference,  which  is  expected.  Consequently,  conclusions  were 
drawn  from  observation  that  parameter  5  is  the  best  for  4.7%  alginate  and  4.7%  alginate  is 
the  best  for  SW  imaging  in  brain. 

In  fig.30, 4.7%  alginate  is  reproduced  using  push  parameter  5  at  four  different 
positions.  According  to  the  elastograph,  the  images  are  consistently  homogeneous  with  a 
steady  mean  G  value  of  ~4.7kPa.  Same  analysis  and  results  were  observed  in  another  4.7% 
alginate.  In  fig.3 1,  static  SW  propagation  movie  in  4.7%  alginate— brain  medium  are  shown. 
Recall  that  fig.3  shows  the  static  SW  propagation  movie  in  alginate.  Comparing  the  two 
propagations,  and  the  ones  in  rats  (fig.23  and  fig.24  middle],  SW  in  mouse  brains  is  totally 
traceable  and  propagates  more  like  in  alginate  as  SW  does  not  break  up  or  get  messy.  Note 
that  the  bottom  part  of  wave  front  does  tilt  a  little  bit  as  SW  propagates  because  of  the 
brain  in  the  medium. 

Fig.32  shows  shear  modulus  images  of  a  control  mouse  brain  in  position  1  using  5 
parameters.  According  to  observation,  the  two  hemispheres  have  slight  difference  in 
parameter  1,  2,  and  3.  In  parameter  4  and  5,  no  difference  can  be  noticed  between  two 
hemispheres.  In  position  2  (fig.33],  parameter  1  results  in  huge  artifacts  within  the  brain 
but  parameter  4  and  5  show  close  shear  modulus  in  both  hemispheres.  In  position  3  and  4 
(fig.34  and  fig.3  5],  there  are  artifacts  existing  in  all  of  the  images  in  the  brain,  and  there  is 
discontinuity  of  shear  modulus.  This  results  primarily  from  the  messiness  of  SW 
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propagation,  which  is  caused  by  the  ultrasound  pushes  generated  into  the  brain  instead  of 
alginate.  Notice  in  images  of  parameter  3  and  5,  the  elastograph  becomes  smooth  and 
homogeneous  again  out  of  the  brain.  Same  imaging  protocol  and  qualitative  analysis  were 
conducted  in  4  more  control  mouse. 

Fig.36  shows  shear  modulus  images  of  a  TBI  mouse  brain  in  position  1  using  5 
parameters.  According  to  observation,  the  TBI  sides  in  all  parameters  are  shown  to  have 
lower  elasticity  than  the  normal  side.  However,  in  parameter  1  and  2,  the  center  of  brain  is 
also  shown  to  be  low  in  stiffness.  Parameter  3,  4  and  5  turn  out  to  have  more  localized  low 
stiffness  region  corresponding  to  the  TBI  side.  In  position  2  (fig.37),  parameter  1  again 
results  in  huge  artifacts  within  the  brain,  while  parameter  2  to  5  all  show  localized  reduced 
elasticity  at  the  TBI  side.  In  position  3  and  4  (fig.38  and  fig.39),  artifacts  again  come  up  in 
the  elastograph.  However,  if  comparing  the  brain  region  in  the  fig.38  and  fig.39,  there  is  still 
an  observable  difference  in  the  two  hemispheres  regardless  of  artifacts.  Same  imaging 
protocol  and  observational  analysis  were  conducted  in  4  more  TBI  mouse. 

According  to  observation  on  5  control  and  5  TBI  mice,  position  1  and  2  are  chosen 
for  further  analysis  because  position  3  and  4  are  way  too  noisy.  Parameter  1  and  2  tends  to 
result  in  artifacts  in  the  brain.  Parameter  3  has  much  fewer  artifacts  but  often  displays 
inter— hemisphere  difference  in  control  brains  and  unfocused  injury  in  TBI  brains.  Parameter 
4  and  5  become  most  prominent  in  diagnosis  of  TBI  because  they  work  well  on  both  control 
and  TBI  brains.  One  hypothesis  for  this  is  that  parameter  4  and  5  have  more  active  pushing 
elements  and  longer  push  duration,  which  induce  larger  and  stronger  tissue  displacement, 
and  enable  SW  to  propagate  smoothly  without  being  dissipated  or  disturbed.  Because  it  is 
hard  to  tell  from  eyes  which  parameter  (4  or  5)  is  the  most  sensitive  to  TBI,  quantitative 
analysis  is  needed. 

Fig.40  [left]  indicates  that  in  control  mouse  brain,  using  parameter  4,  no  statistically 
significant  difference  is  observed  between  position  1  and  2  in  both  contralateral  and 
ipsilateral  hemispheres  (p  =  0.516).  This  means  that  the  time  and  position  variance  do  not 
alter  the  brain  property  great  enough  for  the  elastograph  to  detect  using  parameter  4. 

Same  result  is  reached  using  parameter  5  (p  =  0.115)  in  fig.40  [right]. 

Fig.  41  compares  the  TBI  and  normal  side  of  TBI  brain  at  position  1  [left:  p  =  0.012] 
and  position  2  [right:  p  =  0.016]  separately  using  parameter  4.  Both  statistics  show 
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significant  decrease  of  shear  modulus  on  the  TBI  side.  Then,  in  fig.42,  the  TBI  and  normal 
side  in  both  positions  are  not  significantly  different  (p  =  0.656  TBI;  p  =  0.735  normal].  The 
result  further  affirms  that  the  time  and  position  variance  do  not  affect  the  elastograph  in 
parameter  4.  Same  statistical  analysis  on  parameter  5  reaches  the  same  conclusion  except 
that  there  is  no  significant  difference  found  in  position  1  of  the  TBI  brain  (p  =  0.065], 
though  there  is  a  great  trend  that  they  will  be  different  according  to  fig.43.  This  may  be  an 
indication  of  insensitivity  of  parameter  5  to  TBI,  but  may  also  be  due  to  the  small  sample 
size  (n  =  5]  involved  in  the  study. 

According  to  the  result  above,  data  from  position  1  and  position  2  of  the  same  brain 
are  combined  and  analyzed  again.  In  fig.45,  both  parameter  4  and  5  indicate  no  inter- 
hemispherical  difference  in  control  brains.  Parameter  4  has  p=0.819,  which  is  more  slightly 
greater  than  p=0.799  in  parameter  5.  In  fig.46,  both  parameter  4  and  5  indicate 
significantly  lower  stiffness  in  TBI  side.  Parameter  4  is  the  more  sensitive  (p<0.001],  then 
parameter  5  (p=0.004].  In  fig.47,  control  hemispheres  using  parameter  4  and  5  are 
compared  to  the  normal  side  of  TBI  brains.  TBI  normal  sides  in  both  parameters  have  a 
higher  mean  shear  modulus  than  the  controls,  indicating  a  trend  of  increasing  stiffness, 
however,  the  statistical  analyses  indicate  no  significant  difference. 

As  a  conclusion  according  to  the  results  above,  the  designed  SW  tracking  algorithm 
successfully  track  and  produce  elastograph  of  alginate  and  brain.  Among  the  different 
methods  developed,  MaxValue  tracking  using  a  frame  lag  of  8  plus  correlation  interpolation 
and  weighted  spatial  average  give  rise  to  the  elastograph  with  high  resolution  and  accurate 
shear  modulus  in  alginate.  Low  percentage  of  alginates  (2%]  did  not  work  with  rat  brains 
due  to  mismatch  of  SW  impedance.  After  trails  of  experiments  and  reproducing  in  mouse 
brains,  4.7%  alginate  matches  best  to  its  SW  impedance  and  allow  smooth  SW  propagation 
into  the  brain.  Push  parameters  5  works  best  in  4.7%  alginates.  In  diagnosis  of  TBI,  push 
parameter  4  appears  to  be  the  most  sensitive.  In  control  brains,  parameter  4  reports  an 
average  shear  modulus  of  6.354kPa  and  6.271kPa  for  left  and  right  hemispheres  with  p  = 
0.819.  In  TBI  brains,  parameter  4  reports  an  average  shear  modulus  of  4.962kPa  and 
6.978kPa  for  TBI  and  normal  hemispheres  with  p  <  0.01.  The  shear  modulus  on  the  normal 
side  of  the  TBI  brain,  though  increases,  is  not  significantly  different  from  that  of  control 
hemispheres.  This  result  is  consistent  with  the  preliminary  study  done  with  SSI  in  that 
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brain  stiffness  all  appear  to  decrease  on  the  TBI  side.  Nonetheless,  this  result  differs  from 
the  preliminary  study  in  that  brain  stiffness  on  the  normal  side  of  TBI  brain  does  not 
appear  to  be  higher.  This  difference  may  be  due  to  the  ex  vivo  imaging  protocol,  (for 
preliminary  study,  rats  are  imaged  in  vivo)  which  is  likely  to  alter  mechanical  properties  of 
brain  due  to  exposure.  Moreover,  the  analysis  of  this  project  is  accurate  to  0.5mm  regarding 
to  the  location  of  TBI.  This  gives  another  hypothesis  of  the  discrepancy  in  shear  modulus  in 
that  preliminary  studies  chose  regions  without  artifacts,  which  could  be  anywhere  in  the 
hemisphere.  The  study  also  gives  reasonable  shear  modulus  values  for  both  control  and  TBI 
brains,  which  has  not  been  done  yet  in  the  existing  literature. 

Future  Work: 

The  first  step  of  future  work  should  focus  on  collecting  more  data  on  mice.  Due  to 
limited  time,  only  5  TBI  and  5  control  mice  were  studied  in  the  project.  More  data  from 
animal  studies  can  create  a  better  normal  distribution  for  ANOVA  tests,  and  thus  give  more 
accurate  average  brain  shear  modulus.  Secondly,  elastograph  should  be  further  optimized 
for  more  focused  diagnosis  of  TBI.  As  due  to  heavy  filtering  and  averaging  of  the 
displacement  data,  the  elastograph  tends  to  spread  low  G  value  to  non  TBI  regions  as  in 
fig.36  parameter  2  and  3.  Making  the  elastograph  more  focusing  will  certainly  increase  the 
contrast  between  TBI  and  normal  hemispheres.  Thirdly,  this  technology  should  be  applied 
onto  rat  brains,  as  they  are  larger  animal  models.  Since  all  the  rodent  studies  will  be  ex  vivo 
because  of  their  limited  skull  size,  the  next  step  of  future  work  should  move  on  to  larger 
animal  models  to  better  represent  human  brains.  The  optimized  SW  tracking  algorithm 
should  also  be  implemented  into  VUE,  so  that  VUE  can  display  elastograph  simultaneously 
with  imaging.  Additionally,  more  sophisticated  TBI  inducing  surgeries  or  procedures  also 
need  to  be  designed,  i.e.  procedures  that  will  cause  massive  and  irregular  destruction  of 
brain  tissue.  Finally,  the  technology  needs  to  be  optimized  so  that  it  is  no  longer  invasive. 
The  technology  can  proceed  on  to  human  clinical  trials  if  all  the  above  works  are  finished 
successfully. 
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